Re^taduced 
luf  the 


ARMED  SERVICES  TECHNICAL  INFffiMffl^  AGENCY 
ARLICTON  HALL  STfflffll 
ARLINGTON  12,  VIRIRNIA 


HOTICE:  When  govemment  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  Sc 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  lAatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  foirmilated,  furnished,  or  in  any  way 
supplied  the  s£ild  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  ri^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


METEOROLOGICAL  rkSPECTS 

of 

RADIO  RADAR  PROPAGATION 


<X.  CO 
CJ>  <c 


^-u,Cu>  -  ^ 

NOX 


A  S  T  i  A 

l^r SEP  27  1360 


FRONTISPIECE  -  The  Radar  Hole. 


METEOROLOGICAL  ASPECTS 


OF 

RADIO -RADAR  PROPAGATION 


appl:*d  ntURCH 

-A-R-#-W-A- 

ORERAYIONAL  WEATHER  ANALYSIS 


U.  S.  NAVY  WEATHER  RESEARCH  FACIUTY 
BUILDING  R-48 
U.  ».  NAVAL  AIR  STATION 
NORFOLK  11,  VIRGINIA 

June  1980 


1 


FOREWORD 


This  report  has  been  prepared  under  Task  31,  “Meteorological 
Aspects  of  Refraction  and  Propagation  of  Radar  Waves,”  which  was 
assigned  to  the  U.  S.  Navy  Weather  Research  Facility  for  the  purpose 
of  developing  new  techniques  for  predicting  atmospheric  propagation 
parameters.  Fleet  and  field  tmits,  ope'  ating  microwave  radio  and  radar 
equipment,  require  atmocrheric  inforiL  ation  in  order  to  evaluate  and 
predict  the  performance  oi  naval  communications  and  detection  systems. 

This  publication  is  intended  as  a  comprehensive  reference  in 
radio-radar  meteorology  for  Naval  Weather  Service  personnel.  It  is 
not  a  procedure.!  manual  and,  more  often  than  not,  suggests  several 
solutions  to  a  given  problem  without  prescribing  the  “best  one”.  Our 
objective  has  been  to  assemble  sufficient  general  information  so  as  to 
facilitate  the  solution  of  new,  as  well  as  routine,  problems  as  they  arise 
in  the  field.  In  an  era  of  rapid  technological  advance  and  of  ever- 
changing  concepts  as  to  the  tactical  employment  of  radio  and  radar  sys¬ 
tems,  new  problems  must  be  dealt  with  on  an  almost  daily  basis. 

The  central  ideas  presented  herein  are  largely  the  result  of  re¬ 
search,  supported  in  part  by  theU.  S.  Navy  Weather  Research  Facility, 
at  Smyth  Research  Associates  in  San  Diego,  California  and  at  the  Boulder 
Laboratories  of  the  National  Bureau  of  Standards  in  Boulder,  Colorado. 
Chapters  5,6,  and  7  were  takenalmost  wholly  from  a  manuscript  written 
under  contract  by  Mr.  L.  J.  Anderson  of  Smyth  Research  Associates. 
Chapters  8  and  9  were  taken  almost  wholly  from  technical  reports  pre¬ 
pared  by  Mr,  B,  R,  Bean  and  associates  at  the  Boulder  Laboratories 
of  the  National  Bureau  of  Standards.  Chapters  1,  2,  3,  and  4  were  as¬ 
sembled  by  the  undersigned  from  current  literature. 

■j.nere&<=r4.r.'.b  program,  which  laid  the  foundations  for  this  report, 
was  organized  under  the  direction  of  Corner.  J.  P,  King,  USN,  with  the 
assistance  of  Mr.  Paul  L.  Hexter.  Excellent  editorial  assistance  '..as 
provided  during  the  assemblv  of  this  publication  by  Mr.  John  M.  Mercer. 


DANIEL  F.  REX 
Commander,  U.  S,  Navy 
Officer  in  Charge 
Na^vy  Weather  Research  Facility 
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INTRODUCTION 


A  specific  responsibility  of  the  meteorological  *'ficcr,  as  stated 
in  U.  S.  Navy  Regulations^,  is  to  provide  "the  informulion  necessary 
for  evaluation  of  the  effects  of  sea  and  atmospheric  conditions  on  radar 
and  sonar  performance,"  In  order  to  provide  this  information  (insofar 
as  it  relates  to  radar  performance)  the  weather  officer  must; 

1.  Understand  the  effects  of  atmos>.  .eric  variables  on  (radio 
and)  radar  performance, 

2,  Know  the  fleet  requirements  for  atmospheric  information 
in  terms  of  current  operational  employment, 

3,  Collect  the  necessary  data,  and  make  those  analyses  per¬ 
tinent  to  a  description  of  existing  atmospheric  propaga¬ 
tion  conditions. 

4.  Describe  the  propagation  conditions  that  will  be  encoun¬ 
tered  at  some  future  time. 

It  should  perhaps  be  emphasized  that  success  can  only  be  achieved  by 
taking  each  of  these  four  actions. 

It  is  the  aim  of  this  manual  to  provide  the  Naval  Weather  Service 
officer  with  information  of  direct  assistance  in  connection  with  items 
1,3,  and  4  above.  A  knowledge  of  fleet  requirements  (item  2)  can  only 
be  obtained  from  experience  and  through  frequent  and  detailed  discus¬ 
sions  with  command  and  operational  personnel.  The  weather  officer 
must  take  the  initiative  in  order  to  fully  inform  himself  concerning  fleet 
(operational;  requirements. 

It  should  bementioaed  here  that  since  the  material  in  this  manual 
was  gathered  from  maxiy  sources,  the  letters  and  sj’mbo]'?  used  in  the 
text  and  equations  are  necessarily  uniformly  defined  throughout  the 
entire  manual;  i.e.,  the  symbols  may  only  be  consistent  in  meaning 
within  the  section  where  they  are  used. 

In  part  I  of  thi.s  manual,  some  of  the  basic  principles  of  radio- 


'li.  3.  Navy  Uegulations,  Chapter  9,  Article  0928. 


xiii 


radar  pi’opagation  are  presented  in  order  to  familiarize  the  reader  with 
the  nature  of  electromagnetic  energy  and  how  it  propagates  through  the 
atmosphere.  No  attempt  has  been  made  to  present  a  complete  treatise 
on  the  subject;  only  that  information  considered  to  be  e>  ential  to  an 
elementary  xmderstanding  of  the  subject  has  been  included,  x-artll  con¬ 
tains  a  rather  general  discussion  of  atmospheric  refraction,  the  methods 
available  for  obtaining  refractive  index  data,  and  the  various  modified 
indices  in  common  use.  In  part  III,  working  methods  for  computing 
radio-radar  performance  are  presented,  tof  jther  with  numerical  ex¬ 
amples  (in  some  cases)  intended  to  clarify  ti.j  procedures  used.  In 
order  to  apply  the  methods  of  part  III  intelligently,  the  reader  should 
become  as  familiar  as  possible  with  the  background  material  presented 
in  parts  I  and  II. 

The  final  section  of  this  manual,  part  IV,  presents  information 
available  from  various  climatological  and  synoptic  studie  3  of  atmos¬ 
pheric  refractivity.  This  material  has  been  included  not  only  because 
of  its  interest  in  connection  with  planning  problems,  but  more  particu¬ 
larly  because  of  its  importance  in  predicting  future  refractive  condi¬ 
tions.  It  seems  unlikely  that  our  ability  to  forecast  pressure,  tempera¬ 
ture,  and  humidity  will  improve  rapidly  enough  within  the  next  decade 
to  permit  one  to  predict  future  propagation  conditions  from  direct  fore¬ 
casts  of  these  parameters,  Refractivity  forecasts  will  more  probably 
be  made,  within  the  foreseeable  future,  from  a  consideration  of  the  ex¬ 
pected  climatological  and  synoptic  influences  on  present  (observed) 
conditions. 

In  concluding  this  introduction  it  may  be  well  to  mention  that  ab¬ 
normal  radio-radar  performance  may  be  the  result  of  improper  or 
abnormal  fxmetioning  of  the  radar  equipment,  or  it  may  be  the  result 
of  atmospheric  effects.  In  the  first  instance  the  electronics  (or  radar) 
officer  must  take  action.  In  the  second  instance,  the  weather  officer 
and  the  electronics  officer  must  take  joint,  coordinated  action  in  order 
to  minimize  detriruental  ''>nd  to  take  advantage  of  favorable)  atmos¬ 
pheric  conditions.  Work  with  your  electronics  -  CIC  -  radar  teammatesi 
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1.  ELECTROMAGNETIC  PROPAGATION 


1.1  Characteristics  of  Electromagnetic  Waves 

The  wave  nature  of  visible  light  was  first  demonstrated  by  Thomas 
Young  in  1801.  James  Clerk  Maxwell,  after  a  study  of  Faraday’s  scien¬ 
tific  papers,  published  several  mathematic.  '  treatises  on  the  electro¬ 
magnetic  theory  of  light  in  the  late  1860's.  It  was  not,  however,  until 
1873  that  this  theory  appea  'ed  in  a  fully  developed  form  in  his  "Elec¬ 
tricity  and  Magnetism’’.  Maxwell  reduced  all  electric  and  magnetic 
phenomena  to  stresses  and  motions  of  a  material  medium  and  proved 
that  electromagnetic  waves  in  a  homogeneous,  transmitting  medium 
tend  to  propagate  in  straight  lines  with  a  velocity  equal  to 


where  c  is  the  velocity  of  light,  ^  is  the  magnetic  permeability,  and  k 
is  the  dielectric  constant  of  the  medium. 


Since  Maxwell’s  time  it  has  been  found  that  there  is  an  entire 
spectrum  of  electromagnetic  waves,  covering  at  least  20  decades  of 
wavelength,  over  which  the  physical  characteristics  of  the  waves  are 
identical,  except  for  wavelength  itself.  A  representation  of  a  simple 
wave  motion  is  shown  in  figure  1.1.  In  this  figure  a  periodic  variable, 
A,  is  plotted  graphically  as  a  function  of  time,  For  the  case  of  an 
electromagnetic  wave  the  quantity  A  may  be  considered  to  be  either  the 
electric  field  strength  or  the  magnetic  field  strength  at  some  fixed  point 
in  space.  If  we  define  one  wavelength,  X,  as  the  distance  from  crest  to 
crest  (or  trough  to  trough)  and  the  frequency,  £,  as  the  number  of  cycles 
completed  per  unit  time;  i,e„  the  number  of  wavelengths  which  pass  a 
fixed  point  in  uni+  time,  the  velocity  of  propagation,  v,  is  equal  to 


V  -  Xf  , 


(1.2) 


or,  from  equation  (1.1) 


c 


Xf  =  .  (1.3) 

^  c 

Ir  hoiiLOg-iiecus  m-.  diuru  the  quantity is  a  constant  (i.e.,  the  prop¬ 
agation  velocity  is  constant)  so  that  'either  wavelength  or  frequency, 
which  are  inversely  related,  may  be  used  to  uniquely  describe  any  par- 
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Figure  1.1.  A  Simple  Sine  Wave. 


ticular  electromagnetic  wave. 

Electromagnetic  energy  propagates  best  in  a  vacuum,  where  both 
H  and  k  equal  unity;  however,  most  gases  and  some  solids  are  trans¬ 
parent  to  it  aswvll.  The  velocity  of  propagation  in  a  vacuum  is  identical 
over  the  entire  wavelength  spectrum  and  is  equal  to  the  velocity  of  light 
(3.0  X  10®  cm. /sec.).  One  of  the  basic  axioms  of  physics  is  that  this 
is  the  maximum  veloci.^r  atiainable  oy  energy  oc  matter-. 

In  some  of  its  properties  electromagnetic  energy  behaves  as  though 
it  consisted  of  energy  packets  (quanta)  instead  of  waves.  This  has  given 
rise  to  the  branch  of  physics  known  as  “quantum  mechanics’’,  which 
has  done  much  to  clarify  many  of  the  puzzling  phenomena  which  have 
been  observed,  particularly  at  extremely  short  wavelengths.  For  the 
purpoi  cs  c.f  ihi?  manual  how-e/er,  one  may  regard  electromagnetic 
energy  as  being  pu'ely  a  wave  phenomenon.  In  the  radio-radar  spec¬ 
trum,  tlie  wave  nature  of  the  energy  is  sufficiently  general  to  describe 
all  the  phenomena  with  which  we  are  concerned. 


-  5  - 


1.2  The  Electromagnetic  Spectrum 


Figure  1.2  is  a  representation  of  the  entire  electi  nagnetic  spec¬ 
trum  as  it  is  now  recognized.  It  is  divided  into  numerous  wavelength 
bands  as  shown.  The  boundaries  are  somewhat  arbixrary  because  over¬ 
lapping  occurs,  particularly  at  shorter  wavelengths,  but  the  figure  illus¬ 
trates  the  variety  of  phenomena  to  which  the  name  “electromagnetic 
energy”  is  properly  applied.  The  wavele:  T+h  and  frequency  scales  are 
logarithmic.  The  vertical  lines  are  divided  into  increments,  each  of 
which  represents  a  tenfold  change  in  wavelength.  The  units  commonly 
used  to  designate  wavelength  and  frequency  are  shown  just  to  the  left 
and  to  the  right  of  the  vertical  lines. 

The  narrow  band  just  below  the  1  micron  (ju)  line  represents  the 
visible  spectrum  which  covers  the  range  of  sensitivity  of  the  human 
eye.  Late  in  the  last  century  the  remainder  of  the  vast  spectrum  (non¬ 
visual)  shown  in  figure  1.2  was  discovered,  resulting  in  a  flood  of  new 
knowledge  and  techniques  whose  usefulness  is  now  gradually  being  ex¬ 
ploited. 

The  portion  of  the  spectrum  discussed  in  this  manual  is  that  shown 
above  the  far  infrared  band  in  figure  1,2,  particularly  the  region  above 
the  1- centimeter  wavelength  line.  This  is  the  radio- radar  spectrum 
in  which  waves  are  generated  by  coherent  electronic  transmitters  of 
various  kinds. 

1.3  Wave  Fronts  and  Rays 

In  dealing  with  the  behavior  of  electromagnetic  energy  it  is  often 
convenient  xo  employ  the  concepts  of  wave  fronts  and  rays.  The  Dutch 
scientist  Huyghens  in  1690  stated  the  principle  which  leads  to  a  simple 
geometrical  construction  of  the  wave  front,  and  which  allows  us  to  deal 
with  a  wave  moti«:>:'.  a*  .  t^raight- ’ine  prop'-gation. 

In  figure  1.3  the  point  P  represents  apoint  source  of  energy  from 
which  electromagnetic  waves  r.  e  proceeding  through  a  homogeneous 
medium  with  equal  speeds  in  all  directions.  At  the  end  of  a  certain 
time  the  disturbance  will  have  arrived  at  all  points  which  lie  on  tht 
spherical  surface  Sj.  Since  the  points  are  all  equidistant  from  P,  the 
locus  ei  such  points  is  _  ..x.here,  and  this  sphere  represents  the  initial 
wave  front.  According  to  Huyghens’  principle  every  point  on  the  wave 
front  becomes  immediately  a  new  source  loum  which  secondary  waves 
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Figure  1.3,  Huyghens’  Construction  of  the  Wave  Front. 

or  wavelets  spread  out  in  all  directions.  By  construction  then,  let  us 
draw  circles  of  equal  radius  from  each  point  on  the  surface  S2  repre¬ 
senting  the  wavelets.  These  innumerable  wavelets,  starting  together 
from  all  points  influenced  by  the  principal  wave,  overlap  and  interfere 
with  each  other  so  that  their  resultant  effect  is  produced  only  at  the 
points  in  the  resultant  surface  which  at  any  instant  envelop  all  the  secon¬ 
dary  wave  fronts.  Thus,  the  new  principal  wave  front  will  be  S2,  which 
is  the  enveloping  surface  or  envelope  of  the  wavelets.  The  tangent  sur¬ 
face,  constructed  from  the  wavelets  according  to  Huyghens’  principle, 
thus  gives  the  position  ot  the  wave  front  av  a  later  time.  The  result  is 
the  same  as  if  the  old  wave  front  had  expanded  into  a  new  •■vave  front 
S2,  the  disturbance  progressing  in  straight  lines  in  all  directions. 

It  can  be  seen  that  in  a  homogeneous  medium  (i.e.,  propagation 
velocity  constant  at  all  points  in  the  medium)  the  enveloping  surfac  0 
or  ,vAVw  front  etn&natmft  J’rom  a  point  source  will  be  a  sphere  concen¬ 
tric  with  all  other  earlier' and  later  wave  fronts,  and  that  straight  lines 
or  rays,  radiating  out  from  the  center,  will  be  the  paths  of  the  distur¬ 
bance.  The  ray  (see  fig.  1.3),  which  is  the  path  of  the  disturbance,  is 
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thus  normal  to  the  wave  front,  which  is  the  envelope  of  the  disturbance. 
1.4  Propagation  Characteristics  in  Free  Space 

It  is  convenient  to  begin  our  study  of  electromagnetic  propagation 
characteristics  with  a  consider!  tion  of  the  behavior  of  electromagnetic 
waves  in  free  space;  i.e.,  in  a  vacuum  and  without  the  presence  of  dis¬ 
crete  bodies. 


i.4.1  Free  Space  Fields 


A  fundamental  factor  governing  the  behavior  of  electromagnetic 
waves  is  the  inverse  distance  spreading  of  energy.  If  a  transmitting 
source  emits  a  given  amount  of  energy  per  second,  £.^,  when  one  is  at 
dist?aice  r  from  the  transmitter  this  energy  is  spread  out  over  the  sur¬ 
face  of  a  sphere  of  radius  £  and  area  47rr2  (see  fig.  1.3).  Thus,  the 
power  received  ^  at  distance  r  from  a  point  source  P  is 


(1.4) 


watts  per  unit  area.  The  power  decreases  as  the  square  of  the  distance 
from  the  source  (or  transmitter),  and  the  electromagnetic  field  strength, 
which  is  proportional  to  the  square  root  of  power,  decreases  as  the 
first  power  of  the  distance.  This  decrease  is  known  as  the  ‘  ‘free  space 
attentuation”  of  the  field,  since  it  occurs  even  in  a  vacuum,  and  without 
the  earth  present.  It  is  a  direct  consequence  of  the  fundamental  law  of 
the  conservation  of  energy,  and  applies  both  in  free  space  and  in  the 
terrestrial  environment. 


1.4.2  Interfereiice 

Whenever  two  (or  more)  wave  trains  travelling  over  different 
paths  impinge  (interser*.)  a;  inr-i'ace,  tlit  phenomenon  known  as 

interference  occurs.  If  the  two  waves  arrive  at  the  pom+  in  question 
in  phase  they  reinforce  one  another  and  produce  a  resultant  field  strength 
greater  than  either  of  the  two  compo  .ent  waves  taken  alone.  If  the  two 
waves  arrive  at  the  point  in  question  out  of  phase  they  partially  cancel 
one  another  and  produce  a  weaker  resultant  field  strength.  This  phe¬ 
nomenon  is  fundam.ental  to  all  wave  moxion  and  occurs  both  in  free  space 
and  in  >'hetcrrost’-iai  eu  riror,ii.e.it.  Figure  1.4  illustrates  the  phenom¬ 
enon  of  interference  for  two  simple  sine  waves  of  the  same  frequency 
(or  wavelength),,  The  figure  shows  the  time  variation  of  the  electro- 


Figure  1.4.  The  Interference  of  Two  Wave  Trains. 
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magnetic  field  strength,  E^,,  at  point  in  space,  which  results  from  the 
interference  of  two  waves  of  equal  frequency  and  with  fi.  H  strengths 
E^  and  The  component  waves  E^  and  Eg  are  shown  in  ^  hase,  90’’ 
out  of  phase,  and  180“  out  of  phase. 
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2.  RADIO  RADAR  PROPAGATION 


2.1  The  Radio-Radar  Spectrum 

That  portion  of  the  electromagnetic  spectrum  utilized  for  com¬ 
munication  extends  from  millimeter  wavelengths  to  those  thonsJinds  of 
meters  long.  This  broad  zone  covers  ai-  ’•he  frequencies  of  radio  and 
radar  transmission  (see  fig.  2,1).  Because  the  effect  of  tropospheric 
variables  is  most  pronounced  on  the  shorter  wavelengths,  our  primary 
interest  in  this  manual  will  be  focused  on  the  short-wave  end  of  the 
radio-radar  spectrum  from  about  1  centimeter  to  10  meters.  At  longer 
wavelengths,  the  ionosphere  becomes  the  controlling  influence,  rather 
than  the  comparatively  thin  troposphere.  At  wavelengths  shorter  than 
1  centimeter,  on  the  other  hand,  the  atmosphere  becomes  relatively 
opaque  for  at  least  two  decades  of  wavelength,  and  long  distance  prop¬ 
agation  is  impractical. 

2.2  Atmospheric  Effects 

When  Marconi  demonstrated  in  1901  that  radio  signals  could  be 
transmitted  across  the  Atlantic  Ocean,  investigations  were  soon  initiated 
to  account  for  the  presence  of  signals  beyond  the  horizon.  Four  prin¬ 
cipal  mechanisms  were  proposed;  diffraction,  atmospheric  refraction, 
atmospheric  scattering,  and  reflection  from  an  elevated  layer  of  ionized 
gases.  The  reflecting  layer  hypothesis,  proposed  by  Ksnnelly  and  Hea-vi- 
side  in  1902  and  proven  by  Watson  in  1919,  is  generally  accepted  as  the 
primary  mechanism  which  permits  extended  range  transmission  at 
wavelengths  greater  than  about  10  meters.  Although  communication  at 
these  longer  wavelengths  is  subject  to  various  anomalies  (fading,  holes, 
etc.)  wMch  may  become  serious  problems  in  practice,  the  principles 
of  ionospheric  transmission  lie  outside  the  scope  of  this  manual. 

In  the  1930’s  new  developments  In  eiecuonic  oscillcilors  made 
shorter  wavelengths  useable  for  practical  communication.  It  was  soon 
observed  that  signals  shorter  than  about  10  meters  were  often  not  re¬ 
ceived  at  ranges  greatly  exceeding  the  horizon  distance.  The  shorter 
the  wavelength  the  more  transparent  the  Kennelly-Hea'vriside  layer  he- 
cemes  to  incident  rad':.tion,  and  the  ionospheric  reflection  mechanism 
jr.coiL.es  reifitively  ummportant  in  explaining  transmission  beyond  the 
horizon.  Once  again  the  attention  of  researchers  was  focused  on  the 
phenomena  of  diffraction,  refraction,  ann  scattering  as  the  mechanisms 


Figure  2.1.  The  Radio-Radar  Spectrum. 
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which  might  explain  extended- range,  short-wave  (<10  m.)  communica¬ 
tion. 


With  the  development  of  radar  in  the  1940's  (on  even  shorter  wave¬ 
lengths),  it  was  foxmd  that  the  shorter  the  wavelength  in  use  the  more 
frequent  were  occurrences  of  propagation  well  beyond  the  radio  horizon. 
Surface  objects  v/ere  often  observed  at  distances  many  times  greater 
than  the  visual  range,  and  these  occurrenv,  '  were  found  to  be  associated 
with  imusual  distributions  of  temperature  and  relative  humidity  in  the 
layers  near  the  earth's  surface.  This  phenomenon  is  known  as  "duct¬ 
ing”  or  "trapping”.  It  is,  in  effect,  an  extreme  case  of  superrefraction, 
or  atmospheric  refraction  greater  than  normal. 

The  phenomenon  of  superrefraction  was  observed  even  prior  to 
1940  for  radio  propagation  jji  the  range  of  1.6  to  5  meters;  comprehen¬ 
sive  studies  of  superrefraction  have  been  principally  stimulated,  how¬ 
ever,  by  the  development  of  radar  during  and  after  World  War  II,  The 
7-meter  and  1.5-meter  radar  sets  established  at  coastal  sites  along 
the  Mediterranean  reported  superrefraction  effects  much  more  intense 
than  those  reported  on  the  same  wavelengths  around  the  British  Isles. 
The  entire  coast  of  Sicily,  for  example,  could  frequently  be  seen  with 
1,5-meter  radars  on  Malta,  and  from  time  to  time  echoes  from  Greece 
and  Sardinia  at  ranges  on  the  order  of  650  kilometers  were  obtained. 
Ships  were  sometimes  plotted  to  a  range  of  300  kilometers  in  locations 
where  the  geometrical  horizon  was  less  than  30  kilometers  away.  Per¬ 
haps  the  most  distant  targets  ever  observed  were  seen  on  a  1.5-meter 
radar  at  Bombay,  India[13],  whichreceu  ;d  echoes  from  points  in  Arabia 
at  a  distance  of  2^00  kilometers.  Although  meteorological  effects  usually 
act  to  extend  surface  radar  coverage,  it  should  be  mentioned  here  that 
subrefraction  also  does  occur.  On  several  occasions,  radar  equipment 
on  Fisher’s  Island,  New  York,  was  imable  for  many  hours  to  pick  up 
Block  Island  despite  the  fact  that  it  is  only  35  kilometers  away  and  ■'/isi- 
ble  to  the  naked  f- ve. 

In  the  years  following  World  War  11  with  the  continuously  increas¬ 
ing  application  of  airborne  rudar,  an  even  more  disturbing  phenom- 
enonhas  been  frequently  observed,  the  "radar  hole”.  This  effect,  pro¬ 
duced  by  the  presence  of  superrefractive  layers  aloft,  may  at  ximr.s 
drastically  reduce  the  effectiveness  of  airborne  search  radar  (see  Fron- 
•.’spiacfti, 

2.3  Propagation  Characteristics  in  the  lerrestrial  Envirornnent 


In  order  to  find  an  explanation  for  the  anomalies  described  in  the 
foregoing  paragraphs  we  shall  consider  in  some  detail  tho'''rjhenomena 
which  control  propagation  in  the  troposphere;  (1)  the  abo./rption  of 
radiant  energy,  (2)  the  interference  of  waves,  resulting  in  regions  where 
two  waves  enhance  and  in  adjacent  regions  where  they  cancel  each  other, 
(3)  the  diffraction  of  waves  around  the  earth  and  other  objects,  and  (4) 
the  variations  in  propagation  velocity  caused  by  changes  in  pressure, 
temperature,  and  humidity;  which  give  rise  to  ^>efractive,  reflective, 
and  scattering  effects, 

2.3.1  Absorption 

Whenever  electromagnetic  radiation  propagates  through  a  real 
medium  (be  it  solid,  liquid,  or  gaseous),  a  fraction  of  the  incident  energy 
is  aosorbed  in  passing  through  each  increment  of  the  medium.  The 
energy  absorbed  is  imparted  to  the  atomic  or  molecular  structure  of 
the  medium  and,  in  most  cases,  is  eventually  given  up  as  heat.  The 
effect  may  be  thought  of  as  a  kind  of  frictional  loss. 

Any  particular  substance  does  not  absorb  radiation  of  different 
wavelengths  equally,  nor  do  different  substances  absorb  radiation  of 
one  particular  wavelength  equally  —  as  anyone  familiar  with  the  spec¬ 
trograph  will  know.  The  particular  absorption  characteristics  of  a 
medium  (e,g,,  air),  therefore,  depend  upon  its  composition,  upon  the 
atomic  and  molecular  structure  of  each  constituent,  and  upon  the  phy¬ 
sical  characteristics  of  the  incident  radiation. 

Fortunately,  in  practice  the  absorption  mechanism  does  not  often 
produce  serious  performance  anomalies  in  the  operation  of  radio  or 
radar  equipment.  Absorption  may,  however,  become  an  important  de¬ 
sign  consideration  in  selecting  the  wavelength  to  be  used  in  a  radar 
system;  e.g,,  for  aii'  zczrr^  radar  one  would  avoid  selecting  a  wave¬ 
length  close  to  one  of  the  water  vapor  absorption  bands,  since  water 
vapor  is  present  in  the  natural  atmosphere. 

2.3.2  Interference 

In  the  radio- radar  spectrum  interference  produces  very  pro¬ 
nounced  f^’ictuations  in  signa’  < '^rength  (see  fig,  1.4).  These  are  not 
only  pro  luceU  by  elec  tromagnelic  wave  trains  from  two  or  more  sources 
(transmitters)  arriving  at  a  particular  point  in  ^pace  (receiver),  but 
also  by  wave  trains  from  a  single  source  arriving  at  a  particular  point. 
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having  travelled  over  different  transmission  paths.  For  examp].e,  two 
different  ray  paths  from  source  to  receiver  might  be  the  direct,  line- 
of- sight  ray  and  a  reflected  ray  from  the  transmitt.  antenna  to  the 
earth’s  surface  and  thence  to  the  receiving  antenna  (see  fig,  2,3), 

Interference  bet w eeit  signals  of  different  wavelengths  from  two  or 
more  sources  is  eliminated  by  the '  ‘tuning  section”  of  the  receiver  which 
is  nothing  more  than  a  sharply  tuned  byp,.  "5  filter.  The  filter  rejects 
all  signals  except  those  having  frequencies  within  a  very  narrov^  band 
centered  on  the  frequency  ;o  which  the  receiver  is  “tuned”.  Any  elec¬ 
tromagnetic  wave  (signal)  having  a  frequency  within  the  bypass  band 
will  be  "passed”  on  through  the  amplification  and  detection  stages  of 
the  receiver;  waves  having  other  frequencies  are  rejected. 

Signals  of  the  same  (or  nearly  the  same)  frequency,  whether  from 
several  sources  or  from  a  single  source,  produce  interference  which 
cannot  be  eliminated  by  a  filter  (tuner).  The  case  of  two  or  more  sources 
is  of  practical  importance  in  connection  with  electronic  counter  meas¬ 
ures  (ECM)  but  will  not  be  discussed  further  here.  The  case  of  inter¬ 
ference  produced  by  surface  reflections  from  a  single  source  will  be 
developed  more  fully  in  the  following  paragraphs, 

A  spherical  earth  enveloped  by  a  homogeneous  atmosphere  is 
shown  in  figure  2.2,  In  this  figure,  as  in  all  others  depicting  curved 
earth  geometry,  the  height  scale  has  been  greatly  exaggerated  com¬ 
pared  to  the  horizon'  i  distance  scale.  Assume  that  a  radio  transmitter 
is  located  at  ^  at  height  H  above  the  earth.  Since  the  atmosphere  is 
homogeneous,  energy  will  leave  the  antenna  in  straight  rays.  If  one  now 
draws  a  straight  ray  from  R  tangent  to  the  earth’s  surface,  this  line  de¬ 
fines  the  horizon  distance  D  and  the  horizon  height  h  at  point  P  which 
is  at  a  distance  d  beyond  the  tangent  point.  Referring  to  figure  2,2, 

a’-  H  D-  -  (a  +  K)"  ;  (2.1) 


expanding,  we  have 

=  2  a  H  +  .  (2.2) 

Since  H  is  much  smaller  than  2  a,  the  second  order  term  in  the  above 
f  "v-.t -  on  ir  .j.y  be  neglected.  Therefore, 


D  =  \/2  a  H  , 


(2.3) 
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Figure  2.2.  Ray  Geometry  for  Spherical  E=>rth — Homogeneous  Atmosphere, 
and  similarly, 

d^sfrTh  .  (2.4) 

The  hor'-zon  line  beyond  the  tangent  point  divides  radio  space  into 
two  regions.  The  region  above  the  line  is  always  illuminated  directly 
by  the  transmitter  and  is  called  the  interference  region;  the  region  be¬ 
low  the  line  is  not  dir***-  :  illuminated  and  is  generally  called  the  dif¬ 
fraction  region  (shaded  in  fig.  2,2)c  The  illuminated  region  is  referred 
to  as  the  interference  region  because  the  radio  field  strength  at  a  given 
point  above  the  horizon  line  is  determined  by  the  interference  between 
two  rays,  as  described  in  section  1.4.2.  One  is  the  direct  ray  from  the 
transmitter  to  the  point;  the  other  is  a  ray  reflected  from  the  earth' s 
surface  to  the  point.  Since  the  reflected  ray  path  is  somewhat  longer, 
the  ray  -u  ive?  at  the  poi  iomev/hat  later  than  the  direct  ray.  If  it 
arrives  half  a  wavelength  behind  (180®  out  of  phase)  the  two  waves  will 
tend  to  cancel  each  other,  and  if  they  are  of  jal  strength  the  resultant 
field  at  the  point  will  be  zero.  On  the  other  hand,  if  the  point  is  higher. 
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so  that  the  reflected  wave  arrives  a  full  wavelength  behind  the  direct 
wave,  the  two  will  be  in  phase,  and  the  field  strength  will  beinlensified. 
This  alternate  cancellation  and  enhancement  in  the  fieiv.  strength  as  the 
path-length  differences  are  increased,  gives  rise  to  a  vertical  inter¬ 
ference  pattern,  consist! ig  of  a  series  of  narrow  fingers  or  "lobes” 
of  strong  field  strength  with  weaker  fields  in  between  (see  fig.  6,1). 


In  figui-e  2.3,  showing  a  spherical  earth  surrounded  by  a  homo¬ 
geneous  atmosphere,  the  path  difference,  6 ,  between  the  direct  and  re¬ 
flected  waves  may  be  shown  to  be 


(2.5) 


where  Hjj  and  Hrj.  are  the  heights  of  the  transmitter  at  point  R  and  of 
the  target  at  point  ^  above  the  reflecting  plane  and  d  is  the  direct  path 
distance.  The  reflection  plane  is  tangent  to  the  earth  at  the  point  of 
reflection,  such  that  the  angle  of  incidence  equals  the  angle  of  reflec¬ 
tion  0 . 


figure  2.3.  Path-Difference  Geometry- -Spherical  Earth. 
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On  a  spherical  earth,  arid  Krp  are  not  the  actual  heights  above 
the  surface.  Referring  to  figure  2.3,  it  may  be  shown  that 


and 


(2.6) 


(2.7) 


where  hj.  and  are  the  heights  of  the  transmitter  and  target  above  the 
earth’s  surface  and  dj^  is  the  distance  from  the  reflection  point  to  the 
projection  of  the  target  on  the  reflection  plane.  Since  the  angle  of  in¬ 
cidence  is  equal  to  the  angle  of  reflection,  we  also  have 


d 


1 


d  -  dj 


(2.8) 


One  can  eliminate  Ht.,  and  di  from  the  above  four  expressions,  and 
solve  for  6  in  terms  of  h^.,  1^,  and  d.  The  resulting  expression,  how¬ 
ever,  is  to'o  complicated  to  warrant  further  discussion  in  this  manual. 
Domb  and  Pryce  [12]  have  developed  a  convenient  means  of  plotting 
interference  patterns,  to  which  the  reader  is  referred  for  practical 
assistance. 


The  picture  just  described  is  further  complicated  in  that  the  re¬ 
flected  ray  unde  •>'goes  asudden  phase  shift  of  one-half  wavelength  (180°) 
on  reflection;  minimum  fields,  therefore,  occur  when  the  path  difference 
is  an  integral  number  of  wavelengths  (1,  2,  3  ...  etc),  and  maxima  are 
centered  at  those  poir>t<.  • '>v!ng  1/2,  11/2,  2  1/2  ...  etc.  wavelength 
path  differences.  Furthermore,  the  strength  of  tns  reflected  wave  is 
reduced  by  scattering  if  the  reflecting  surface  is  rough  and  by  diver¬ 
gence  if  it  is  smooth,  so  that  the  direct  and  reflected  waves  are  not  of 
equal  strength.  These  factors  modify,  but  do  not  invalidate,  Ihe  general 
characteristics  of  the  interference  pattern  as  described  above. 

In  rro'-tice,  interfere-v; .  patterns  for  a  given  antenna  are  com¬ 
puted  for  various  angles  of  elevation  of  the  antenna  (tilt);  these  com¬ 
puted  values  cire  checked  and  corrected  with  .!  oservational  data  \/hen- 
ever  possible.  The  weather  officer  should  consult  the  electronics  offi- 
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cerfor  further  information  .^oncermag  the  characteristics  of  a  particu¬ 
lar  antenna  array.  The  phenomenon  of  interference  is  not  affecteci  per 
se  by  changing  meteorological  conditions,  except  that,  we  shall  see, 
refractive  effects  may  produce  path  length  changes. 

2.3  c3  Diffraction 

The  only  way  in  which  energy  can  pea  .trate  the  region  below  the 
horizon  line,  in  a  homogeneous  atmosphere,  is  by  the  process  of  dif¬ 
fraction  (see  fig.  2.2),  This  process  res\ilts  from  the  wave  nature  of 
the  transmitted  energy  and  is  similar  to  simple  diffraction  over  a  knife 
edge.  The  mathematics  of  diffraction  over  a  spherical  obstacle  is  quite 
formidable  and  will  not  be  attempted  here.  It  turns  out,  however,  that 
the  depth  of  penetration  of  energy  below  the  horizon  line  is  proportional 
to  the  wavelength,  and,  at  a  constant  height  above  the  surface,  it  falls 
off  exponentially  with  distance  beyond  the  horizon  point.  Domb  and 
Pryce  [12]  have  developed  relatively  simple  methods  for  computing 
diffraction  fields,  and  the  reader  is  referred  to  their  paper  for  a  dis¬ 
cussion  of  practical  computational  techniques. 

As  was  the  case  with  interference,  the  mechanism  of  diffraction 
is  not  primarily  weather-dependent,  and  as  such  its  effect  upon  radio¬ 
radar  performance  is  more  properly  a  problem  for  the  electronics 
engineer. 

2.3.4  Reflection  and  Refraction 

We  now  complicate  the  simple  case  of  figure  2.2  and  2.3  by  sur¬ 
rounding  spherical  earth  with  a  nonhomogeneous  atmosphere  (such 
as  the  earth's  real  atmosphere).  Inhomogeneties  produce  changes  in 
the  speed  of  propagation  of  electromagnetic  energy  through  different 
portions  of  the  an'ior; -'cric  envelope,  giving  rise  to  the  phenomena  of 
reflection,  refraction,  and  scattering. 


The  index  of  refraction,  n,  of  a  medium  is  defined  as 

■>  =  7  .  (2.0' 

tvne.  .-  c,  is  the  prcpagatio*' velocity  of  electromagnetic  energy  in  a  vac¬ 
uum,  and  V  is  the  propagation  velocity  in  the  medium  in  question.  The 
index  of  refraction  of  a  medium  is  therefor  i  a  measure  of  the  propaga¬ 
tion  speed  (relative  to  the  velocity  of  light}  of  an  electromagnetic  wave 
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through  the  medium. 

(a)  Reflection  and  Refraction  at  a  Discontinuity 

Whenever  an  electroms  gnetic  wave  strikes  a  surface  of  discon¬ 
tinuity  in  n  (i.e.,  a  surface  across  which  there  is  an  abrupt  change  in 
propagation  speed)  the  incident  wave,  in  gen  ral,  is  partially  reflected 
and  partially  refracted  (see  fig.  2.4).  The  lav  of  refraction  at  a  dis¬ 
continuity  may  be  derived  from  a  consideration  of  propagation  speeds 
in  the  two  media.  In  figure  2.4  we  assume  electromagnetic  energy 
strikes  the  boundary  between  two  media  with  indices  of  refraction  of 
nj  and  n^,  where  n^  is  associated  with  the  initial  medium  and  ng  is 
greater  than  n^;  i.e.,  the  propagation  speed  is  greater  in  the  initial  than 
in  the  final  medium.  In  a  certain  length  of  time,  the  wave  front  tra¬ 
vels  from  position  1-2  to  position  4-3,  and 

S2_3  =  V|  t  ;  S2_4  =  V2  t  ,  (2.10-2.11) 

where  is  the  distance  traveled  and  v^  is  the  speed  of  propagation 
in  medium  1,  and  82.4  and  V2  are  the  corresponding  distance  and  speed 
in  medium  2. 


Also,  we  have 


sin  ©2  =  cos  /l32  = 


’2-3 

’1-3 


(2.12) 


and 


Therefore, 


sin  Oj  =  cos  A13  =  - 


1-4 


’1-3 


(2.13) 


sin  0, 


sin  8- 


’2-0 

n-4 


Ill 

^2* 


(2.14) 


and,  r  =  we  may  wr:?-  The  law  of  refraction 


sin  0. 

_ ^ 

sin  0,, 


n. 


(2.15) 
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where  ^  is  the  angle  of  incidence  and  is  the  angle  of  refraction. 

It  may  be  well  to  note  that  in  figure  2.4  the  dir'^^+ion  of  propaga¬ 
tion  may  be  reversed  without  invalidating  the  conclusions  just  reached 
so  long  asv/e  interchange  the  designation  of  incident  and  refracted  wave 
fronts  and  the  angles  of  incidence  and  refraction.  The  incident  and  re¬ 
fracted  rays  lie  in  the  same  plane,  and  the  incident  rays  are  always 
bent  toward  the  medium  of  higher  index  c.  "efraction  (i.e.,  lower  pro¬ 
pagation  speed)  in  passing  through  a  surface  of  discontinuity. 

With  the  aid  of  Huyghens’  principle  one  may  derive  the  law  of  re¬ 
flection  which  states  that  the  incident  and  reflected  rays  lie  in  the  same 
plane  and  the  angle  of  incidence  is  equal  to  the  angle  of  reflection;  i.e., 
0.  =  0J,  in  figure  2.4.  It  will  be  noted  that  the  incident  wave  front  is 
inverted  by  reflection  (the  incident  wave  is  shifted  in  phase  180°),  as 
has  been  mentioned  previously;  in  refraction  there  is  no  corresponding 
phase  shift. 

Whenever  energy  is  passing  through  a  discontinuity  between  a 


Figure  2,4.  Refraction  and  Reflection  at  a  Surface  of  Discontinuity. 


-  22  - 


medium  of  higher  index  of  refraction  to  a  medium  of  lower,  the  phe¬ 
nomenon  of  total  reflection  will  occur  when  the  angle  of  incidence  ex¬ 
ceeds  a  certain  critical  value  given  by 

8-  =  sin~^  ,  (2.16) 

1  ^2 

where  n2  >  n^^  (see  fig.  2,4).  At  angles  of  int  '••'nee  gree.ter  than  this 
critical  value  the  refracted  ray  is  suppressed  (does  not  emerge  above 
the  surface  of  the  discontinuity)  and  the  incident  ray  is  totally  reflected. 

The  relative  amounts  of  energy  carried  in  the  reflected  and  in  the 
refracted  waves  are  dependent  upon  the  transmission  characteristics 
of  the  two  media  on  either  side  of  the  discontinuity,  upon  the  physical 
characteristics  of  the  surface  of  discontinuity,  upon  the  angle  of  inci¬ 
dence  of  the  initial  wave  train,  and  upon  the  wavelength  of  the  incident 
wave.  For  example,  a  sheet  of  ordinary  screen  wire  is  a  good  reflector 
of  10-centimeter  radar  waves  but  a  poor  reflector  of  electromagnetre 
energy  in  the  visual  spectrum;  whereas,  a  silvered  mirror  is  a  good 
reflector  of  visible  radiation  but  a  poor  reflector  of  gamma  rays. 

(b)  Refraction  in  an  Environment  of  Continuously  Varying  Inde 
of  Refraction 

If  an  electromagnetic  wave  traverses  a  region  in  which  the  prop¬ 
agation  velocity  (and  therefore  also  the  index  of  refraction)  varies  from 
point  to  point  but  in  a  continuous  fashion;  i.e.,  no  discontinuities  in  n 
present,  there  will  be  no  reflection  phenomena,  but  refraction  of  the 
wave  will  occur. 

In  figure  2.5  the  dashed  vertical  line  at  the  left  represents  a  wave 
front  at  time  t  =  0.  The  wave  is  traveling  horizontally  toward  the  right. 
Assume  that  the  refracTiv..  ’  ;dcx  is  be'^tom  of  the  wave  front 

and  is  n^  -  An  at  the  top.  As  the  wave  proceeds  along,  the  +on  part  v.  ill 
travel  faster  than  the  bottom  part,  be^cause  the  refractive  index  at  the 
top  is  less  than  at  the  bottom.  Afte"  time  At,  the  top  end  of  the  wave 
will  have  traveled 


+  Ad 


cAt 

'^top^^  ~  Oq  -  An  ' 


(2.17) 


where  c  is  the  propagation  velocity  in  vacuum.  The  bottom  end  will 
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Figure  2,5.  Refraction  in  anEnvironment  of  Continuously  Varytng  Index  of  Refraction. 


have  travelled 


cAt 

“  ^bottom  ~  Hq 


(2.18) 


and, thus 


Ad  =  cAt 


1 _ 

-  An 


-  An 


(2.19) 


Since  {n^  -  An)  is  very  close  to  unity  we  may  disregard  it  and  write 


Ad  =  dpAn  . 


(2.20) 


Since  the  lop  of  the  wave  has  travelled  A^  farther  than  the  bottom,  the 
wave  front  is  now  tipped  downward  by  an  angle  a  =  Ad/ Ah,  where 
is  the  distance  increment  between  the  top  and  bottom  of  the  wave  front. 
Substituting  Ad  =  o>.s  the  eycpi'ession 


a 


An 

Ah 


(2.2il 


for  the  angle  through  which  the  ray  is  bent. 


But  is  simply  the  curvature,  of  the  ray  and  since  thebend- 
i.-g  ui-op';i iionai  co  .i.j.ance,  the  ray  path  is  an  arc  of  a  circle  of 
radius  R  =  ,  so  long  asthe  n-gradient,  ■^,  remains  constant.  Thus, 

one  can  say  that. 
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C 


R 


An 

Ah 


(2.22) 


This  is  a  very  important  expression  in  radio  wave  jcefraction; 
the  curvature  of  a  ray  is  equal  to  the  gradient  of  the  refractive  index 
normal  to  the  ray  path.  Staten  another  way,  the  angle  through  which  a 
ray  will  be  bent  in  passing  through  a  refractive  index  gradient  is  equal 
numerically  to  the  gradient  normal  to  the  path,  '  'mes  the  distance  tra¬ 
velled. 


.  In  the  natural  atmosphere  the  index  of  refraction  becomes  smaller 
as  one  proceeds  to  higher  and  higher  altitudes,  but  abrupt  changes  (dis¬ 
continuities)  in  n  are  not  commonly  observed.  For  this  reason  refrac¬ 
tive  effects  in  the  earth’s  atmosphere  are  generally  more  important  in 
the  study  of  anomalous  radio- radar  performance  than  are  reflective 
effects.  As  has  already  been  implied,  however,  the  surface  of  the  earth 
acts  as  a  good  reflector  of  electromagnetic  waves  in  the  radio- radar 
spectrum. 

(c)  Snell’s  Law 

One  is  frequently  interested  in  tracing  the  trajectory  of  a  ray  as 
it  travels  through  a  medium  such  as  the  atmosphere  in  which  the  re¬ 
fractive  index  changes  with  height.  Theory  indicates  that  the  ray  will 
follow  a  trajectory  having  the  shortest  optical  path.  This  principle  is 
expressed  in  a  simple  mathematical  form  known  as  Snell’s  Law,  It  was 
originally  applied  to  ray  tracing  through  lenses  and  prisms,  but  applies 
equally  well  in  the  radio- radar  case.  If  one  assumes  that  the  refrac¬ 
tive  index  in  the  atmosphere  changes  only  with  height  above  the  spher¬ 
ical  earth,  and  not  with  horizontal  distance,  then  Snell’s  Law  may  be 
expressed  as  ' 


'  r.  (a  -  n)  cos  p  ~  a  cos  ,  (2.23) 

where  n  is  the  refractive  index  at  a  given  height  h,  n^  is  the  refractive 
index  at  the  earth’s  surface,  a  is  th..-  earth’s  radius,  is  the  ray  incli¬ 
nation  angle  at  height  h,  and  is  the  ray  inclination  at  the  earth’s 
surface.  Figure  2.6  shows  a  typical  ray  leaving  the  earth  at  an  angle 
Gf~  and  passing  upwards  through  the  atmosphere.  By  the  time  it  has 
arriv.'d  ai  h'^igh'.  h,  it  ha  ^  been  Dent  downward  by  a  total  angle  t,  which 
can  be  calculated  -'f  one  knows  the  way  in  which  n  changes  with  height. 
This  figure  will  be  referred  to  later  when  the  method  for  computing 
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total  bending,  t,  and  the  associated  radar  elevation  angle  error,  is 
developed. 

2.3.5  Scattering 

The  presei-Mie  thermodynamically  dissimilar  parcels  of  air  in 
the  atmosphere  is  responsible  for  a  propagacion  mechanism  which  has 
not  been  considered  previously  and  which  is  known  as  “scattering” 
It  has  been  found  experimentally  that  when  one  proceeds  away  from  a 
transmitter  with  a  receiver  near  the  surface,  the  fields  beyond  the  hor¬ 
izon  decrease  exponentially  at  a  rate  which  agrees  with  diffraction 
theory  for  quite  some  distance  (see  fig.  7.1),  As  one  proceeds  furtluT 
!UVc„'  om  the  transmit*-  r,  however,  the  attenuation  rate  suddenly  de¬ 
creases  tthe  fields  decrease  with  distance  at  a  much  slower  rate).  In 
other  words,  beyond  a  certain  critical  di  tance,  the  receiver  fields  are 
much  stronger  than  can  be  accounted  for  by  diffraction  theory.  The 
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mechanism  by  which  these  fields  are  produced  is  called  scattering. 
Each  parcel  above  the  horizon,  which  is  illuminated  by  the  transmitter, 
scatters  a  small  part  of  the  incident  energy  in  directio*.  other  than 
along  the  straight  line  path  from  the  transmitter  through  the  parcel. 
Most  of  the  cnarg^^.  of  course,  continues  straight  through.  If  a  particu¬ 
lar  parcel  is  above  the  horizon  for  both  the  transmitter  and  the  re¬ 
ceiver,  part  of  the  scattered  energy  will  be  directed  toward  the  re¬ 
ceiver.  The  scattered  field  at  any  point  may  Le  regarded  as  the  sum 
of  the  energies  scattered  by  a  multitude  of  such  parcels  ar'd  directed 
toward  the  point  in  question.  Theories  have  been  developed  for  com¬ 
puting  the  fields  to  be  expected  by  such  a  mechanism;  by  posxulating 
reasonable  An  values  and  a  spectrum  of  parcel  sizes,  one  can  compute 
fields  which  agree  well  with  those  observed.  One  of  the  current  prob¬ 
lems  in  tropospheric  propagation  research,  centers  around  proving  or 
disproving  this  theory  of  scattering.  Refractometer  data  show  that  at 
the  heights  required  for  the  parcels  to  be  visible  from  both  terminals 
over  long  paths,  the  actual  n  fluctuations  are  smaller  than  those  re¬ 
quired  by  the  theory  to  explain  the  observed  fields.  This  has  led  to  the 
postulation  of  a  multiscattering  mechanism,  in  which  the  energy  is  scat¬ 
tered  from  parcel  to  parcel  several  times  before  reaching  the  receiver. 
In  such  a  model,  the  only  parcel- height  requirement  is  that  a  particular 
parcel  be  in  the  line  of  sight  of  the  previous  and  the  succeeding  parcel, 
and  not  necessarily  be  in  the  line  of  sight  of  both  terminals.  This  per¬ 
mits  the  scattering  to  occur  lower  in  the  atmosphere  where  the  ^  val¬ 
ues  are  much  greater.  The  only  difficulty  with  this  proposal  is  that 
an  analytical  treatment  of  multiscattering  becomes  very  complex  and 
has  in  fact  not  yet  been  worked  out.  A  further  complication  occurs  when 
one  considers  that,  because  of  atmospheric  stability,  the  parcels  are 
probably  flattened  like  pancakes,  rather  than  being  globular.  A  great 
deal  of  experimental  evidence  exists,  however,  which  seems  to  support 
the  multiscattering  hypothesis  with  flattened  parcels. 

In  any  cas»  scatter  fields  do  exist;  they  are  stronger  at  night  than 
in  midafternoon,  and  are  stronger  in  summer  than  in  winter.  Even  at 
their  weakest,  they  are  much  strong»>r  than  diffraction  fields  over  paths 
greater  than  a  few  hundred  mil**''  '=  length. 

2.4  Summary 

.n  this  chapter  we  have  described  briefly  some  of  the  operational 
anomalies  actually  observed  in  the  performanc  ?  of  radio- radar  equip¬ 
ment  and  have  discussed  the  physicr^  phenomena  which  are  important 
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in  producing  the  anomalous  performance  observed.  Of  these  phencmena, 
absorption,  interference,  :'.:id  diffraction  are  not  significantly  al'fected 
by  varying  atmospheric  parameters  and  do  not,  therefore,  produce  crit¬ 
ical  weather  problems  in  the  field. 

On  the  other  hand,  refraction,  reflection,  and  scattering  are  most 
commonly  (except  in  the  case  of  reflection)  produced  by  inhomogeneities 
in  the  atmosphere  which  in  turn  are  attri’’.utable  to  variations  in  mete¬ 
orological  parameters.  Of  these,  by  far  thi  most  important  in  explain¬ 
ing  or  predicting  the  anomalous  performance  of  naval  radio- radar  equip¬ 
ment  is  refraction.  It  is  for  this  reason  that  the  remaining  three  parts 
of  this  manual  are  devote  J  principally  to  a  detailed  discussion  of  atmos¬ 
pheric  refraction. 


PART  II 


THEORY  OF  ATMOSPHERIC  REFRACTION 
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3.  THE  ATMOSPHERIC  REFRACTIVE  INDEX 

3,1  Definition 


It  was  mentioned  previously  that  the  propagation  velocity  of  all 
electromagnetic  energy  in  a  vacuum  is  a  onstant  (c  =  3,0  x  10^  cm./ 
sec.).  In  atmospheric  gases  this  velocity  is  1  educed  slightly.  The  ratio 
of  the  velocity  of  propagation  of  electromagnetic  energy  in  a  vacuum 
to  that  in  a  particular  medium,  is  defined  as  the  refractive  index  of  the 
medium  as  was  described  in  section  2,3.4,  This  quantity,  designated 
as  n,  generally  varies  with  wavelength  and,  of  course,  from  medium  to 
medium.  In  tropospheric  air  (within  the  wavelength  band  from  1  cm. 
to  10  m.)  the  refractive  index  is,  for  ail  practical  purposes,  independent 
of  wavelength  and  has  a  value  slightly  greater  than  unity.  At  sea  level, 
values  range  from  approximately  1.000250  to  1.000450,  with  values  near 
1.000350  frequently  found  over  an  ocean  surface. 

To  facilitate  numerical  computation  it  is  more  convenient  to  de¬ 
fine  a  derived  index,  called  the  “refractivity",  as  follows 

N  =  (n  -  1)  X  10®  .  (3.1) 

The  atmospheric  refractivity  ranges  from  250  to  450  N-units  at  sea 
level.  It  is  also  possible  to  express  N  as  a  function  of  total  atmospheric 
pressure,  temperature,  and  water  vapor  pressure,  as  follows 


where  P  and  e  are  in  millibars  and  ^  is  in  degrees  absolute.  The  con¬ 
stants,  A  and  B,  ’-?en  determined  experimentally  by  a  number  of 
investigators;  A  is  the  dielectric  constant  for  dry  air  and  ^  is  the  water 
vapor  dipole  moment  constant.  Smith  and  Weintraub  [22J  have  evaluate^' 
the  results  obtained  by  numero’^s  investigators  and  have  recommended 
the  value  of  77.6  for  A  and  the  value  of  4,810  for  B.  Equation  3.2  is 
frequently  written  in  the  following  form: 

N=  -.  V.6(-|-)+3.73  X  10®(-^)  ,  (3.3) 

where  the  first  and  second  terms  are  called  the  “dry”  and  the  “wet” 
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terms,  respectively.  Since  both  P  and  £  normally  decrease  with  height 
above  sea  level,  N  also  decreases  with  height.  The  exact  way  in  which 
N  varies  with  height  is  of  course  determined  by  meteor -’'^gical  pro¬ 
cesses  and  has  a  profound  effect  upon  radio-radar  propagation. 

In  the  event  that  available  moisture  data  are  expressed  as  rela¬ 
tive  humidities,  equation  3.3  may  be  written  follows: 

N  =  77.6  (y)  ^  3.73  x  10^  (-^^)  ,  (3.4) 

where  eg  is  the  saturation  vapor  pressure  in  millibars  and  RH  is  the 
relative  humidity  in  percent. 

It  may  be  well  to  point  out  that  both  n  and  N  are  space  and  time 
functions  whi.ch  in  general  vary  from  place  to  place  in  the  atmosphere 
and  which  in  general  vary  with  time  at  a  particular  point.  The  pre¬ 
ceding  formulae  permit  one  to  determine  the  refracstivity  for  one  par¬ 
ticular  place  and  time;  namely,  the  place  and  time  represented  by  the 
values  of  pressure,  temperature,  and  humidity  used.  In  general  the  path 
which  a  radio-radar  ray  will  follow  is  primarily  determined  by  the  gra¬ 
dient  of  N  rather  than  by  N  itself  as  was  discussed  in  section  2.3.4(b). 
N-gradients  may  be  determined  from  a  series  of  discrete  values  of  N 
in  the  region  of  interest,  or  the  gradients  may  be  determined  by  direct 
measurement  with  suitable  equipment. 

3.2  Methods  of  Measurement  and  Computation 


Direct  measurements  of  the  propagation  velocity  of  electromag¬ 
netic  energy  in  a  particular  medium  can  be  made  and,  when  compared 
with  the  speed  of  light,  give  a  direct  measure  of  the  refractive  index 
or  refractivity.  Because  of  the  complexity  of  the  equipment  needed, 
this  approach  is  impractical  for  operational  use. 

3.2.1  Refractometer 

A  more  practical  way  to  measure  the  atmospheric  refractive  index 
was  developed  in  the  1940’s  by  Birnbaum  [7]  at  the  National  Bureau  of 
Stand.:aJ'’  "ud,  somewhat  in  a  slightly  different  form  by  Crain 

and  lean  jSj  at  the  Umversity  of  Texas.  The  refractometer,  as  the 
instrument  is  caLed,  essentially  measures  changes  in  the  dielectric 
constant  of  the  air  inside  a  perforated  cylindrical  cavity.  This  is  done 
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by  feeding  electromagnetic  srergy  of  about  a  3 -centimeter  \vavelength 
into  a  cavity  which  has  been  machined  to  resonate  over  a  very  narrow 
band  of  wavelengths.  Changes  in  the  dielectric  constant  of  the  air  in¬ 
side  the  cavity  cause  the  resonant  frequency  of  the  nty  to  change 
proportionally.  By  using  two  cavities  (one  being  hermetically  sealed 
Sind  containing  either  dry  a.r  or  vacuum  and  the  other  being  perforated 
to  allow  passage  of  outside  air)  one  can  detect  very  small  changes  in 
the  resonant  frequency  of  the  perforated  '"dxy.  Since  the  refractive 
index  in  this  case  is  simply  the  square  root  of  the  dielectric  constant, 
these  changes  can  be  recorded  as  changes  in  N.  The  sensing  lag  of  the 
device  is  limited  only  by  the  time  needed  to  introduce  new  air  into  the 
cavity.  Since  the  volume  of  air  contained  in  the  cavity  is  very  small, 
the  speed  of  response  may  be  reduced  to  a  fraction  of  a  second  by  care¬ 
ful  design  of  the  ventilating  ducts  or  ports. 

In  spite  of  various  simplifications  which  have  been  made  in  re¬ 
cent  years,  the  refractometer  is  a  relatively  complex  and  heavy  instru¬ 
ment,  which  is  not  yet  suitable  for  expendable  balloon  borne  operations. 
Its  main  use  has  been  aboard  aircraft,  and  many  hundreds  of  vertical 
and  horizontal  profiles  of  refractive  index  have  been  recorded,  particu¬ 
larly  in  the  continental  United  States  and  Alaska,  Very  recently,  re¬ 
fractometer  s  have  been  installed  in  AEW-type  aircraft  and  in  certain 
carrier- based  types.  In  the  future  refractometer  data  may  be  expected 
to  become  available  to  most  large  operational  forces.  It  should  always 
be  used  in  preference  to  derived  N-data,  whenever  there  is  a  choice, 
because  of  its  greater  precision  and  because  it  reproduces  the  "fine” 
structure  of  atmospheric  refractivity. 

3.2,2  Temperature,  Pressure,  and  Humidity  Method 

Since  N  is  a  function  of  P,  T,  and  e,  it  follows  that  whenever 
measurements  of  these  quantities  are  available  they  also  serve  to  de¬ 
fine  N,  At  the  Si, t fare  ihe  CC44.\  Jtional  Lcrometer,  thei’mometer,  and 
hygrometer  may  be  used  to  determine  N;  aloft  radiosond**  provides 
the  necessary  measurements.  From  values  of  T,  and  e,  N  may  Le 
computed  using  equations  3.3  or  3.4.  The  particular  features  of  the 
radiosonde  which  make  it  useful  for  propagation  purposes  are  that  it 
measures  the  required  quantities  as  functions  of  height,  and  that  recoi 
data  are  available  which  have  been  obtained  over  periods  of  many  years 
!•  id  trsTU  ‘■.andreds  sutions  having  a  wide  geographic  distribution. 

All  meteorologists  recognize  that  the  radiosonde  has  certain  de- 
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ficiencies  which  make  it  somewhat  less  than  ideal,  particularly  in  the 
measurement  of  humidity.  The  high  lag  coefficient  of  the  humidity 
sensing  element,  particularly  at  low  temperatures,  causes  the  resulting 
N  (namely,  the  height  profile)  to  deviate  from  its  true  shap.  The  effect 
of  these  deviations  on  radio  propagation  perhaps  is  not  as  serious  as 
it  may  be  for  other  purposes.  The  reason  for  this  is  that,  as  will  be 
shov/n  later,  anomalous  propagation  effects  are  principally  conditioned 
by  the  gross  features  of  the  N-profile;  the  fine  structure  is  not  of  first 
order  importance.  This  does  not  mean  that  for  propagation  work,  im¬ 
provements  in  the  radiosonde  are  unnecessary;  rather,  they  are- greatly 
desired  (see  section  3.3.3).  It  does  mean,  however,  that  record  radio¬ 
sonde  data  can  be  used  to  good  advantage  in  spite  of  its  shortcomings. 

Comparisons  between  N-values  obtained  from  refractometer 
measurements  and  those  computed  from  radiosonde  data  are  in  good 
agreement  as  to  the  gross  structure  present,  but  of  course  the  fine  de¬ 
tails  are  largely  missed  in  the  computed  radiosonde  values.  Figure 
3.1  is  an  example  of  such  a  comparison.  Where  appreciable  deviations 


Figure  3.1.  Comparison  of  Brjiosonde  and  Refractometer  N-Profiles. 
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occur,  it  has  been  found  th;\t  thesffv;ere  caused  by  pronounced  humidity 
fluctuations,  as  would  be  expected. 

In  propagation  problems,  one  is  continually  faced  v»*ih  the  problem 
of  reducing  P“»  T-,  and  e-data  to  N- values  at  numerous  heights.  Over 
the  years,  various  methods  have  been  devised  for  computing  N,  These 
take  the  form  of  tables,  slide  rules,  anp’.ogue  computers,  and  nomo¬ 
grams, 

(a)  Refractive  Index  Nomogram 

Perhaps  the  most  convenient  of  these,  for  use  when  great  pre- 
cisionis  not  required,  is  the  nomogram  shown  infigure  3.2.  Pressure, 
temperature,  and  humidity  are  shown  in  the  units  commonly  used  in 
radiosonde  data  tabulations;  i.e.,  pressure  in  millibars,  temperature 
in  “C.,  and  dewpoint  in  °C.  It  is  particularly  convenient  in  that  the  fig 
ure  and  a  straight  edge  are  all  that  is  required.  Assuming  no  error  in 
entry  data,  the  accuracy  of  the  N- values  so  derived  is  ±  1  N-unit,  which 
is  compatible  with  the  accuracy  of  radiosonde  data.  In  addition  to  fig¬ 
ure  3.2,  a  larger,  full-size  version  of  the  refractive  index  nomogram 
is  provided  as  a  loose-leaf  inclusion  to  this  manual  and  is  marked  as 
chart  I. 

The  required  simple  instructions  for  use  of  the  nomogram  appear 
thereon.  For  example,  suppose  the  value  of  pressure  is  1,000  milli¬ 
bars,  the  temperature  is  20'’C.,  and  the  dew  point  is  4‘’C.  First,  find 
the  intersection  of  the  1,000  millibar  pressure  line  with  the  20®C.  tem¬ 
perature  line.  Then  connect  this  point  with  the  4®C.  division  on  the  dew¬ 
point  scair  and  read  off  an  N- value  of  300  where  the  straight  edge  inter¬ 
sects  the  refractive  index  (refraciiAaty)  scale, 

(b)  Refra^ii'v..-  Diagram 

Another  method  for  determining  r  efractivity  irom  pressure,  ten 
perature,  and  humidity  data  is  illustrated  in  figure  3.3.  This  diagram 
has  been  designed  for  use  in  aircraft,  in  that  the  arguaments  used  are 
pressure  altitude  in  feet,  dry  bulb  temperature  in  °C.,  and  wet  bulb  tem¬ 
perature  in  °C.  Values  of  N  may  be  read  directly  from  the  diagram 
viii.  oprec’sicii  of  ±1  unit,  assuming  the  entry  data  is  without  error. 
In  addition  to  figure  3.3,  a  larger,  full-size  version  of  the  refractive 
index  diagram  is  promded  as  a  loose-le.  f  inclusion  to  this  manual  and 
is  marked  as  chart  11. 
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Figure  3.2.  Refractive  Index  Nomogrami.  (Chart  I) 
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Figure  3.3.  Refractive  Index  Diagram.  (Chart  ID 
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tise  of  the  diagram  is  explained  by  the  inset  in  the  upper  right  - 
hand  corner  of  figure  3.3  and  ir.  the  simple  instructions  which  appeair 
on  the  back  of  chart  II.  For  example,  suppose  the  pressure  altitude  is 
13,000  feet,  the  dry  bulb  temperature  is  10°C.  and  the*\  ^  bulb  tem¬ 
perature  is  -2®C.  First  locate  the  10“C.  and  -2°C,  points  along  the  hori¬ 
zontal  13,000  foot  pressure  altitude  line.  Then  follow  the  straight, 
sloping  dry  bulb  temperaijre  curve  up  and  to  the  left  from  the  10®C« 
point  and  the  curved,  concave- downward  wet  bulb  temperature  curve 
up  and  to  the  left  from  the  -2‘*C.  point,  until  ii.ese  two  temperature 
curves  intersect.  At  the  point  of  intersection  read  off  an  N-value  of 
175  from  the  curved,  convex- downward  refractivity  curves. 

(c)  Refractive  Index  Overlay  (Arowagram) 

Whenever  pressure,  temperature,  and  moisture  data  are  already 
displayed  in  the  form  of  a  plotted  sounding  on  the  Arowagram  (the  ther¬ 
modynamic  chart  in  common  use  in  theU.  S.  Navy),  it  is  most  convenient 
to  obtain  N-values  for  various  levels  of  the  sounding  by  using  the  re¬ 
fractive  index  overlay  (see  fig.  3.4).  This  overlay  has  been  especially 
prepared  for  use  with  the  large-size  Arowagram;  it  cannot  be  used  with 
other  forms  of  thermod3rnamic  charts. 

Use  of  the  overlay  is  explained  in  the  following  instructions  and 
recommended  computational  procedures: 

1.  Plot  temperature  (T)  and  dew  point  (T^j)  on  Arowagram 

and  place  overlay  in  register. 

a.  Carefully  place  a  dot  on  the  overlay  in  register 

with  each  plotted  point  (T,  T^j)  with  a  grease  pen¬ 
cil.  Connect  each  pair  of  temperature-dew-point 
dots  at  the  same  pressure  level  with  a  line  so  as 
to  identify  of  values. 

b.  If  the  aew-point  data  is  missing  extend  c.  line  +he 
left  of  the  plotted  T- value  for  about  10°C. 

c.  If  “motorboating”  is  reported  for  the  dew  point; 
determine  the  saturation  mixing  ratio  at  tempera¬ 
ture  T  and  multiply  by  15  percent  to  obtain  a  “mo¬ 
torboating”  mixif^i-  ratio.  Plot  the  dew  point  for 
thia  moiorboating  mixing  ratio  at  the  same  pres¬ 
sure  level.  Connect  the  points  with :  line  as  before. 
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2.  Remove  overlay  and  place  on  a  white  surface.  Deter¬ 
mine  N  by  following  the  projection  scheme  given  in  the 
inset  at  the  upper  right-hand  coi-nerof  figure'^  4,  mark¬ 
ing  each  point  with  a  grease  pencil  of  a  differt-l  color 
than  the  one  used  for  the  T-  and  Tj-dots.  Connect  each 
dot  as  it  is  determined  to  the  previous  dot  with  a  straight 
line.  This  is  to  avoid  inverting  th'-  sequence  of  the  dots. 

3.  Read  the  value  of  N  from  the  curves  and  place  this  value 
alongside  each  dot  as  it  is  read.  Interpolate  to  the  near¬ 
est  unit  value  of  N. 

4.  Relabel  the  temperature  lines  in  the  lower  left-hand  cor¬ 
ner  of  theArowagram  as  follows:  -80®C.  as  100;  -60°C. 
as  200;  -40‘'C.  as  300;  -20“C.  as  400,  and  plot  the  N- 
values  as  determined  in  (3)  above  at  the  correct  pres¬ 
sure  level  on  theArowagram.  Connect  these  points  with 
straight  line  segments.  This  is  the  N-profile  as  a  fimc- 
tion  of  pressure.  This  plot  maybe  directly  (and  easily) 
compared  with  the  temperature  and  dew-point  curves 
which  are  already  plotted  further  to  the  right.  The  N- 
profile  may  thus  be  related  to  the  existing  synoptic  sit¬ 
uation. 

The  refractive  index  overlay  was  completed  prior  to  the  survey 
of  the  constants  A  and  JB,  which  was  mentioned  in  section  3.1,  and  the 
values  A  =  74.4  and  B  =  4,973  were  used  in  its  construction  rather  than 
the  currently  accepted  correct  values  (77,6  and  4,810).  As  a  result  N- 
values  as  determined  from  figure  3.4  are  consistantly  too  low  (5  to  15 
N-units)  and  should  not  be  used  together  with  N-values  determined  by 
any  other  means  in  computing  N-gradienls  in  a  layer.  So  long  as  both 
N-values  (top  and  bettom  of  the  layer)  are  determined  from  figure  3,4, 
the  error  introduced  in  the  computed  value  of  the  N-gradient  isentirel;/ 
negligible.  As  will  be  seen  in  later  paragraphs,  it  is  tne  N-gradient 
(not  the  discrete  values  of  N)  which  affects  the  performance  of  radio¬ 
radar  equipment.  For  this  reason  the  refractive  index  overlay  may  be 
used  for  most  practical  computations  in  spite  of  the  inherent  errors 
contained  in  individual  N-values  as  determined  from  it. 

In  atidUxOn  lo  figure  3.4,  a  larger,  full-size  overlay  (for  use  with 
the  large  Arowagram)  has  been  printed  on  '  ‘'ansparent  paper  and  is 
provided  as  a  loose-leaf  inclusion  to  this  manual  (chart  III).  Although 
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chart  III  can  be  used  directly  in  its  present  form,  it  is  suggested  that 
working  copies  be  made  on  transparent  plastic  stock  using  the  Ozalid 
(or  other)  duplicating  process.  Routine  use  of  cha  '  m  in  its  present 
form  will  soon  result  in  its  destruction, 

3.2.3  Radar  Method 

Since  variations  in  the  vertical,  refractive  index  structure  in  the 
troposphere  are  the  principal  cause  of  anomalous  radar  propagation 
one  can  use  the  measured  characteristics  of  radar  field  strength  versus 
height  and  range  to  deduce  an  "effective”  refractive  indexprofile.  The 
details  of  an  operational  method  for  doing  this  are  described  in  section 
6.4.  There  are  certain  advantages  and  some  disadvantages  to  using 
this  approach.  The  principal  advantage  is  that  the  profile  deduced  is  an 
effective  profile  averaged  over  many  miles,  rather  than  one  represen¬ 
tative  of  a  single  set  of  points  as  a  radiosonde  ascent  would  be.  Fur¬ 
thermore,  one  would  expect  such  a  profile  to  have  a  better  correlation 
with  propagation  effects,  since  it  was  derived  from  actual  propagation 
measurements.  The  disadvantages  of  this  approach  are  that  it  requires 
radar  equipment  and  aircraft  for  making  the  measurements,  and  that 
one  can  deduce  only  ceridn  gross  features  of  the  refractive  index  pro¬ 
file,  rather  than  its  detailed  structure.  There  are  operational  situations, 
however,  when  current  radiosonde  data  are  not  available,  in  which  the 
radar  method  can  be  very  useful, 

3.3  Refractivity  Gradients 

Most  tropospheric  propagation  anomalies  are  caused  by  the  bend¬ 
ing  (refraction)  of  radio-radar  waves  and  such  bending  is  in  turn  caused 
by  gradients  in  the  refractive  index,  rather  than  by  the  absolute  value 
of  N  itself.  In  fact,  in  an  atmosphere  of  constant  N  no  refraction  occurs 
irrespective  of  .ht  -  u  of  N.  Since  gradients  in  pressure,  tempera¬ 
ture,  and  humidity  exist  in  the  atmosphere,  ^f^fractivity  gradients  must 
also  occur. 

By  far  the  strongest  and  most  persistent  gradients  in  refractivity 
are  observed  in  the  vertical  plane.  The  largest  and  most  persistent 
contribution  to  the  decrease  in  refractivity  v/ith  height  is  the  pressure 
gradient  wmeh  ?mc:.u.ts  to  about  35  millibars  in  the  first  1,000  feet 
above  sea  level.  For  a  similar  pressure  change  to  take  place  in  the 
horizontal  plane,  even  in  association  w.  :h  a  deep  low  pressure  center, 
would  require  a  distance  of  many  hundreds  of  miles.  Similarly,  vertical 
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gradients  of  temperature  and  humidity  are  usually  much  stronger  than 
those  in  the  horizontal  plane.  For  these  reasons,  in  the  study  of  prop¬ 
agation  anomalies  we  are  primarily  concerned  with  grac’^-'nts  in  the 
vertical  plane.  As  will  be  apparent  to  the  meteorological  officer,  these 
vertical  gradients  are  strongly  affected  by  meteorological  processes 
such  as  advection,  subsidence,  turbulence,  radiation,  and  evaporation. 

3,3,1  Standard  Refractivity  Profile 

If  one  averages  out  in  time  the  local  variations  in  refractivity  at 
various  heights,  there  results  what  is  defined  as  the  "standard  refrac¬ 
tivity  profile",  such  as  that  shown  in  figure  3.5.  In  the  troposphere, 
the  vertical  decrease  in  N  is  exponential,  reaching  a  certain  constant 
value  of  N  at  the  tropopause,  regardless  of  the  surface  value.  In  the 
stratosphere  a  single  exnonential  curve  is  sufficient  to  describe  the 
average  profile.  The  mathematical  form  of  each  of  the  two  segments 
of  this  profile  may  be  written  as 

—ch 

Njj  =  Ng  exp.  ,  (3.5) 

where  ^  is  the  refractivity  at  height  h,  Ng  is  the  refractivity  at  the 
earth's  surface,  h  is  height  in  thousands  of  feet,  and  c  is  a  constant 
which  depends  upon  Ng.  If  one  assumes  a  constant  value  of  N  at  the 
tropopause  of  61.0  units  at  a  height  of  42,500  feet,  the  constant  c,  in 
the  troposphere,  may  be  expressed  as  the  following  function  of 

In  Ng 

c  =  -  0.0968  .  (3.6) 

Values  of  £  corrr  spondingto  several  values  of  Ng  are  given  in  table  3.1. 

TABLE  3.1 


Ns 

c 

N- Units 

Per  Thousand  Feet 

465 

0.0478 

400 

0.0442 

350 

0.0411 

318 

0.0388 

'Ut; 

0.0375 

273 

0.0351 

250 

0.0331 

Values  of  the  Exponential  Coefficient  (c)  for  Various  Values  of  Surface  Refractivity  (Ng). 
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Figure  3.5.  Standard  Refractivity  Profile. 
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The  stratospheric  segment  of  the  standard  refractndty  profile  may  be 
expressed  as 


Njj  =  465  exp. 


-0.0478  h 


(3,7) 


where  h  is  the  height  in  thousands  of  feet- 


This  two- segment  form  of  the  standard  re^  activity  profile  was 
derived  from  data  presented  by  Schulkin  [20]  and  has  been  found  to  fit 
observations  quite  adequately.  It  should  perhaps  be  noted  that  the  stan¬ 
dard  refractivity  profile  here  defined  is  not  necessarily  identical  with 
either  the  U.  S,  or  the  ICAO  Standard  Atmospheres^  By  selecting 
equal  to  the  value  computed  from  the  U,  S.  Standard  Atmosphere  sur¬ 
face  values  (U.  S.  Standard  Dry,  Ng  =  273  and  U,  S.  Standard  60%  RH, 
Ng  =  318),  an  N-profile  is  obtained  which  agrees,  in  the  troposphere, 
with  that  obtained  by  calculating  N-values  from  U,  S.  Standard  Atmos¬ 
phere  pressure,  temperature,  and  humidity  data  (using  equation  3.3  or 
3.4). 


One  may  construct  the  standard  refractivity  profile  described 
above,  using  semilogarithmic  plotting  paper  as  follows: 

1.  Layout  height  on  the  linear  scale  and  refractivity  on  the 
logarithmic  scale. 

2.  Plot  the  tropopause  point  at  N  =  61.0,  h  =  42,500  feet. 

3.  Plot  a  point  at  N  =  465,  h  -  0  feet. 

4.  Lay  a  straight  edge  on  these  two  points  and  draw  in  the 
stratospheric  segment  above  h  =  42,500  feet. 

5.  Plot  a  point  at  the  desired  Ng  at  h  =  0  feet. 

6.  Lay  a  straight  edge  on  this  point  and  the  tropopause  point 
and  draw  in  the  desired  tropospheric  segment  below  h  = 
42,500  feet. 


P  shou.’...»  1.:  rr'ilized  that  the  r'lr 'cttire depicted  in  figure  3.5  is  a  long¬ 
term  a\  iragfc,  and  that  significant  departures,  particularly  at  low  alti¬ 
tudes,  are  to  be  expected.  Furthermore  the  he  'jht  of  the  tropopause 
varies  seasonally  and  with  latitude,  but  these  departures  will  not  cause 
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serious  deviations  from  thr  N-profile  as  described  above.  The  b€:havior 
of  radio  propagation  in  a  standard  atmosphere  will  be  discussed  in  de¬ 
tail  in  the  next  chapter. 


The  gradient  of  refractivity,  N,  in  a  standard  atmosphere,  maybe 
obtained  from  the  general  equation  of  the  standard  refractivity  profile 
by  differentiation  (of  equation  3.5)  with  r  jspect  to  h  as  follows: 


dh 


or 


(3.8) 


dN 

dh 


-cN  exp. 


-ch 


dN 

where  is  given  in  N-units  per  thousand  feet,  Ng  is  the  surface  re¬ 
fractivity,  c  is  the  coefficient  given  in  table  3.1,  and  h  is  the  height  ex¬ 
pressed  in  thousands  of  feet. 


We  see  that  the  N-gradient  is  negative  in  the  standard  atmosphere; 
i.e,,  the  value  of  N  decreases  with  height,  and  that  the  magnitude  of  the 
gradient  also  decreases  with  height.  Values  of  N-gradient  for  several 
altitudes  in  the  standard  troposphere  are  given  in  table  3.2;  surface 
values  for  pressure  and  temperature  were  taken  from  U.  S,  Standard 
Atmosphere  with  60  percent  RH;  i.e„  Ng  =  318.  It  will  be  noted  that 
the  N-gradient  in  the  standard  atmosphere  falls  from  an  average  value 
of  about  minus  12  units  per  thousand  feet  in  the  lowest  layers  to  an 
average  value  of  about  minus  6  units  per  thousand  feet  in  the  midtropo¬ 
sphere. 


TABLE  3,2 


“  1 

h 

— 

uN 

Feet 

N-Units /Thousand  Feet 

Surface 

-12.3 

1,000 

-11.9 

-  8.4 

20,000 

-  5.6 

30,000 

-  3.8 

Values  of  N-Gradient  in  che  Standard  Atir.ospherd  (Ng  =  318)  for  Several  Alti*udes  in 
the  Troposphere. 
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3.3,2  Actual  Atmosphere 

In  order  to  evaluate  N-gradients  in  the  actual  a."  "osphere  it  is 
recommended  that  a  plot  of  N  versus  height  be  constructed.  The  height 
scale  may  be  constructed  in  erms  of  geometric  height,  pressure  alti¬ 
tude,  pressure,  etc.,  depending  upon  the  application  intended.  For  most 
purposes  a  simple  linear  plot  of  N  versus  ge>  '»'etric  height  or  pressure 
altitude  will  suffice.  It  is  recommended  the  ordinary  10x10  per  inch 
graph  paper  be  ruled  off  with  height  (2,000  feet  to  the  inch)  as  the  ordi¬ 
nate  and  with  refractivity  (50  N-units  to  the  inch)  as  the  abscissa.  Such 
anN-profile  is  shown  diagramatically  in  figure  3.6.  In  order  to  rapidly 
evaluate  the  N-gradient  in  the  various  layers  cf  the  profile  an  N-gra- 
dient  overlay  has  been  prepared  and  is  shown  in  figure  3.7.  This  over¬ 
lay  has  been  designed  for  use  with  N-profiles  constructed  as  shown  in 
figure  3.6;  it  may  be  hSed  on  any  N-profile  chart  having  linear  h-  and 
N-scaleswith  aratio-^-  equal  to 40.  Inaddition  tofigure3.7,  an  over¬ 
lay  has  been  printed  on  transparent  paper  and  is  provided  as  a  loose- 
leaf  inclusion  to  this  manual  (chart  IV),  As  in  the  case  of  chart  III  it 
is  suggested  that  working  copies  be  made  on  transparent  plastic  stock. 

The  overlay  is  used  by  matching  one  of  the  sloping  lines  on  the 
overlay  with  the  slope  of  a  section  of  the  plotted  N-profile  (the  N-  and 
h-reference  lines  on  the  overlay  must  be  parallel  with  the  vertical  and 
horizontal  lines  on  the  graph  paper),  'i he  value  of  the  selected  sloping 
line  is  the  estimated  N-gradient  for  the  section  of  the  profile  concerned. 
This  process  can  be  repeated  until  the  gradients  of  all  layers  are  de¬ 
termined. 

In  ordfci  to  classify  the  refractive  properties  of  layers  conven¬ 
iently,  four  zones  are  defined  as  follows; 

1.  Subrefra-.-tivc  ^- >  0);  ^  increases  with  height.  Rays 
curve  upward  with  reference  to  a  Gt;‘?.’2ht  line  (opposite 
in  sense  to  the  curvature  of  the  earth’s  surface).  Radio¬ 
radar  ranges  are  signi^’lcantly  reduced;  occurrence  quite 
rare, 

2.  Normal  (0  >'2^  >-  [  ooo^*  ^  decreases  with  height. 

Rc'Vo  curve  dov.'i./<-'drd  (in  the  same  sense  as  the  curva¬ 
ture  of  the  earth’s  surface)  but  not  as  sharply  as  in  zone 
(3).  Radio-radar  performance  is  f.onerally  undisturbed. 
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FjgLi-e  3.6.  Analyzed  N-Pr-Ji'iie. 


NT  OVERLAY 


Figv.  rn  3.7.  N-Gradienf  Overlay  (Linear  Plot).  (Chart  IV) 
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3.  Superrefractive  i^qOO  ^  "  'i,00b)>’  N  decreases; 

with  height.  Rays  curve  downward  (in  same  sense  as  the 
curvature  of  the  earth’s  surface)  more  shar^Wthan  zone 
(2)  but  not  as  much  as  the  curvature  of  the  earth’s  sur¬ 
face.  Radio-radar  ranges  are  significantly  extended; 
occurrence  quite  frequent. 

48  AM 

4.  Trapping(-  1  000  N  Rays 

curve  downward  more  sharply  than,  and  in  the  same 
sense  as,  the  curvature  of  the  earth’s  surface.  Radio¬ 
radar  performance  is  greatly  disturbed;  ranges  extended 
greatly,  appearance  of  radar  holes,  etc.;  occurrence  not 
frequent. 

In  practice  construct  the  N-profile  from  refractometer  values  or 
from  radiosonde  data  and,  using  the  overlay  (fig.  3.7),  evaluate  the  N- 
gradient  in  each  layer.  Label  each  layer  with  the  gradient  value  and 
designate  the  type  of  refractive  zone  by  color  shading  in  accordance 
with  the  color  code  given  on  the  overlay  (see  fig,  3,6  as  an  example). 
All  superrefractive  and  trapping  layers  may  then  easily  be  identified 
from  an  inspection  of  the  completed  N-profile.  Each  of  these  layers 
should  be  examined  for  possible  influences  on  the  particular  radio- radar 
system  in  operation,  using  the  methods  described  in  part  III  of  this 
manual.  For  computational  purposes  use  the  actual  N-values  obtained 
for  the  base  and  top  of  each  layer,  not  the  N-gradient  as  read  from  the 
overlay  (fig.  3.7). 

3.3.3  Limitations  in  Using  Radiosonde  Data 

In  using  radiosonde  data,  as  distinguished  from  refractometer 
data,  the  meteorological  officer  must  keep  certain  shortcomings  of  the 
system  in  mind,  -f.cr'.’ly,  are: 

1.  The  large  lag  (long  response  time)  of  the  humidity  ele¬ 
ment  which  has  a  tendency  to  reduce  the  values  of  all 
recorded  gradient?,  to  increase  the  apparent  thickness 
of  all  superrefractive  and  trapping  layers,  and  to  raise 
the  apparent  height  of  the  top  and  bottom  of  all  layers. 

This  lag  Ume  principally  a  function  of  temperature; 
decreasing  with  increasing  temperature. 

2.  The  sequential  method  of  transmitting  temperature. 
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humidity,  and  refer«n-jes  v/hich  make  it  entirely  possi¬ 
ble  to  miss  a  superrefractive  or  trapping  layer  of  any¬ 
thing  less  than  about  300  feet  in  thicloiess.  The  response 
of  the  temperature  element,  while  the  humidity  ^  ^ment 
is  “off  the  air”,  can  at  best  only  given  an  indirect  clue 
as  to  the  presence  ox  a  refractive  layer,  since  N  varies 
only  slowly  with  temperature. 

3.  The  poor  low-temperature  response  of  the  humidity  ele¬ 
ment.  At  temperatures  below  -40®C.  the  response  of 
this  element  virtually  ceases  which,  along  with  the  se¬ 
quential  system,  makes  evaluation  of  refractive  layers 
in  the  middle  and  upper  troposphere  nearly  impossible 
with  the  present  instrument. 

4.  Incases  whei  aN-profiles  are  plotted  from  coded  radio¬ 
sonde  reports,  even  when  both  standard  and  significant 
levels  are  available,  errors  may  be  introduced  due  to 
the  application  of  the  Circular  "P”  instructions.  The 
criteria  defining  a  significant  level  from  the  humidity 
trace  operate  in  such  a  manner  as  to  often  mask  out  the 
all-important,  rapid  decreases  in  humidity.  When  local 
soundings  are  being  evaluated,  this  deficiency  may  be 
overcome  by  instructing  all  raob  observers  to  be  alert 
for  such  transients  and  to  evaluate  humidity  significant 
levels  critically. 

3.3.4  Radiosonde  N-Gradient  Overlay  (Arowagram) 

In  order  to  minimize  the  effects  of  lag  in  the  computation  of  N- 
gradients,  an  overlay  for  use  with  the  Arowagram  (fig.  3.8)  has  been 
prepared  by  Comdr,  T C.  Clarke,  USNR,  and  associates  at  U.  S.  Fleet 
Weather  Facility,  Argentia,  Newfoundland  [24].  Ihis  overlay  has  been 
designed  so  as  to  permit  the  analysis  of  N-profiles  plotted  on  the  large- 
size  Arowagram  as  described  in  section  3.2.2(c). 

The  principal  parameter  in  radio- radar  meteorology  which  de¬ 
termines  ray  geometry  (as  has  been  shown)  is  the  gradient  dN/dh  or 
AN !  4b  where  h  is  the  true  o-  vecmetric  height  increment.  The  eval¬ 
uation  -f  AN/di"-  from  N-.xrofiles  plotted  as  shown  in  section  3.2.2(c); 
i.e.,  refractivity  v(  rsus  pressure;  requires  re'erence  to  a  pressure- 
height  curve  for  the  particular  sounding  in  question.  It  has  been  shown 


Figure  3.8.  Radiosonde  N-Gradient  Overlay  (AROWAGRAJyl).  (Chart  V) 
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[24]  that  within  the  range  of  current  radiosonde  inaccuracies,  the  gra¬ 
dient  AN/ Ah  may  be  replaced  by  AN/AZp  {Zp  =  pressure  altitude)  with 
sufficient  accuracy  for  computational  purposes. 

We  have  by  definition 

Ah  =  AZp  +  AD  .  (3.9) 

then 


(3.10) 


where  D  is  the  altimete"  correction  and  S  is  the  specific  temperature 
anomaly  after  Bellamy  [5],  Possible  values  of  S  range  from  about  plus 
100  feet  per  1,000  feet  to  minus  150  feet  per  1,000  feet,  with  plus  or 
minus  50  feet  per  1,000  feet  being  probable  extremes  in  the  lower  tro¬ 
posphere.  Hence,  the  approximation  equals may  be  accepted 
without  introducing  errors  greater  than  about  5  percent. 


If  the  radiosonde  actually  reported  true  values  of  temperature 
and  humidity  as  functions  of  pressure,  the  identification  of  trapping  lay¬ 
ers  would  require  simply  the  construction  of  an  overlay  with  isolines 
drawn  for  a  gradient  of  minus  48N-imits  per  1,000  feet  (Zp)  at  suitable 
intervals;  AN/AZp  equal  to  minus  48  N-units  per  1,000  feet  being  the 
critical  gradient  for  trapping  to  occur.  A  visual  comparison  between 
the  slope  of  these  isolines  and  the  slope  of  various  sections  of  the  N- 
profile  as  plotted  in  section  3.2.2(c)  would  delineate  trapping  layers. 
Because  of  the  effects  of  radiosonde  humidity-and-temperature- element 
lag,  as  described  in  section  3.3.3,  many  real  trapping  layers  are  not 
detected  by  this  method. 

Consider  a  radiosonde  ascending  at  a  fixed  rate  ■  =  R  (a  con¬ 
stant),  then 

AN  an/ At  1  /  AN\ 

-  = -  ,  (3JI) 

AZp  AZp /At  R  At  / 

which  shows  that  *he  gradient  of  N  is  a  direct  function  of  AN /At.  From 
equation  3.2  we  have 


and  by  differentiation  with  respect  to  time 


dN  _  A  /dP\  ^  2ABeN.  ^  ^  / de\  ..  , 

dt  TVdty  \t2  ■  t3  /  dt  T2\dt/  * 

Clarke  [24]  has  skew'll  that  fora  typical  ascensional  rate(R  =  300 
m./min.  =  16.4  ft. /sec.)  the  contribution  of  the  first  term  on  the  right- 
hand  side  of  equation  3.12  ranges  from  about  minus  0,16  N-units  per 
second  at  the  surface  to  minus  0.10  N-units  per  second  at  500  millibars. 
Similarly,  the  contribution  of  the  second  term  in  equation  3.12  is  less 
than  one  N-unit  per  second  in  magnitude,  while  under  typical  tropo¬ 
spheric  conditions  the  third  term  may  have  values- as  large  as  minus 
10  N-units  per  second.  Thus,  the  contribution  of  the  first  and  second 
terms  may  be  neglected  in  comparison  with  the  third,  and  we  have 


AN 

At 


(3.13) 


Since  most  real  trapping  layers  are  quite  shallow  (200  to  400  ft.)  we 
may  assume  that  no  significant  temperature  changes  will  occur  within 
the  layer.  Therefore,  from  equation  3.11  we  have 


where  the  “constant”  C  = 
change  of  N  is  principally 
pressure. 


AN 

AZ 


^  C 


(— ) 


(3.14) 


P 
AB 

•  Thus,  we  have  shov/n  that  the  rate  of 
a  function  of  the  rate  of  change  of  the  vapor 


The  humid’iy  ..  ’'=Trent  in  the  present  radiosonde  is  a  relatively 
inaccurate  and  slow-responding  sensor.  Its  lag  in  a  changing  environ¬ 
ment  is  a  function  of  the  temperature,  the  vapor  pressure,  and  the  d.- 
rcction  of  change  of  the  relative  humidity.  Of  these  factors,  tempera¬ 
ture  is  the  most  important,  and  it  may  be  taken  as  the  only  factor  in¬ 
fluencing  the  humidity  sensor's  lag,  without  introducing  appreciallo 
error. 


Using  the  approximations  described  in  the  proceeding  paragraph 
and  using  an  error  function  of  the  form 
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(3.1.5) 


where  c  is  the  error  in  N  introduce(i  by  lag  and  J  is  the  ’•'g  coefficient 
in  seconds,  Clarke  [24]  has  derived  an  expression  for  the  apparent  N- 
gradient  in  a  layer  (as  measu-ed  by  the  radiosonde)  required  for  trap¬ 
ping.  This  expression  is 


trapping 


^ apparent 
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=  -  48  (l  -  e^)  .  (3.16) 
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where  R  istheascentional  rate  in  feet  per  second.  Ah  is  the  layer  thick¬ 
ness  in  feet,  and  ^  is  the  humidity  sensor  lag  coefficient  in  seconds. 

In  order  to  compute  the  apparent  critical  N- gradient  from  equa¬ 
tion  3.1 6,  ^  must  be  specified  and  this  requires  a  knowledge  of  the  tem¬ 
perature  structure.  The  ICAO  Standard  Atmosphere  (S  =  O)  can  be  used 
as  an  initial  reference,  but  since  the  variation  of  ^  with  is  large, 
auxiliary  standard  atmospheres  for  S  =  ±50  and  S  =  +  100  must  also  be 
en  lOyed  to  provide  the  range  of  temperature  structures  found  in  the 
real  atmosphere.  Typical  values  of  the  apparent,  critical  N-gradients, 
given  the  N-units  per  1,000  ieet  are; 

S  =  -50  S  =  0  S  =  ±50  S  =  +100 

For  Zp  =  0  -15  -?.-6  -36  -42 

For  Zp  =  10,000  feet  -5  -12  -21  -31 

If  one  remembers  that  the  true,  critical  N-gradient  is  minus  48  N-units 
per  1,000  feet,  the  large  error  introduced  by  the  lag  of  the  humidity  cle¬ 
ment  is  at  once  apparent. 

li.crder  iocv>nstruct  the  radiosonde  N -gradient  overlay  (fig.  3,8), 
the  apparent,  cri.ical  N-g.radient  was  compu  cd  at  selected  pressure 
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levels  and  was  plotted  on  anArowagram  at  the  correct  N  and  prtsssure 
values.  Since  the  N-values  at  a  given  temperature  T  are  a  function  of 
the  humidity,  and  since  the  ICAO  Standard  Atmosphere  makes  no  pro¬ 
vision  for  moisture;  three  values  of  relative  humidity  were  used;  0  per¬ 
cent,  60  percent,  and  100 percent  in  order  to  cover  tht  ’nge  of  all  pos¬ 
sible  situations.  In  orde  that  heights  remain  unchanged,  the  -sdrtual 
temperature  must  be  used.  Smooth  equally- spaced  curves  were  drawn 
parallel  to  the  computed  gradients.  Th*-  se  curves  define  three  areas 
on  the  N  versus  P  plot  on  the  Arowagram  corresponding  to  the  three 
atmospheres;  S  =  +  50,  S  =  0,  and  S  =  -  50.  It  was  found  convenient  to 
extend  each  area  to  the  right  so  as  to  include  the  adjacent  portion  of 
the  next  area;  e.g.,  S  =  -i-  50  was  extended  to  the  right  so  as  to  niulude 
apparent  gradients  for  a  S  =  +  100  atmosphere. 

Use  of  the  radiosonde  N-gradient  overlay  is  explained  in  the  fol¬ 
lowing  instructions  and  recommended  computational  procedures: 

1.  Carefully  place  the  overlay  in  register  over  the  Arowa¬ 
gram.  The  S  =  0  curves  are  aligned  initially. 

2.  Examine  the  temperature  sounding  (T  vs,  P)  carefully 
to  determine  which  layers  of  the  sounding  fall  within  the 
S  =  0  parallelogram.  Place  grease  pencil  tab  marks  on 
the  N-profile  so  as  to  properly  delineate  this  (these) 
layer(s).  Repeat  this  procedure  for  those  layers  which 
fall  into  the  S  =  +  50  and  S  =  -  50  areas, 

3.  Layers  within  the  S  =  0  parallelogram. 

a.  Carefully  compare  the  gradient  (slope)  of  each 
straight  line  segment  of  the  N-profile  with  the  slope 
of  the  apparent  trapping  gradient  shown  on  the  over¬ 
lay.  If  the  slope  of  the  N-profile  is  greater,  then 
that  layer  is  marked  as  being  a  trapping  layer, 

b.  If  i>  segment  of  iht  i'l-profii^,  is  exactly  parallel  to 
a  trapping  gradient  line  on.  tbo  c  vcrlay,  mark  it  as 
being  a  trapping  layer. 

c.  If  segments  of  the  N-proiile  have  slopes  close  but 
not  equal  to  the  apparent  trapping  gradient,  these 
layers  may  be  trapping  layers.  Examine  the  sound¬ 
ing  to  sCv  there  are  meteorological  indications 
which  suggest  that  a  trapping  layer  should  be  pres¬ 
ent.  If  so,  evaluate  the  layei  (layers)  as  a  trapping 
layer(s). 
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4.  Layers  S  =  +  50  and  S  =  -  50. 

a.  Invert  the  overlay.  Place  the  appropriate  curve 
in  register  and  proceed  to  evaluate  as  in  (3)  above. 

Each  trapping  layer  identified  should  be  examined  for  possible  influ¬ 
ences  on  the  particular  radio- radar  system  in  operation,  using  the 
methods  described  in  part  III  of  this  manual.  For  computational  pur¬ 
poses  use  the  actual  N-value  observed  at  the  la^  :r  base  and  the  actual 
critical  N-gradient  ( =  -  48  N-units/1,000  ft.)  to  determine  the  most 
probable  N-value  at  the  layer  ;op.  Heights  of  the  base  and  top  (and 
layer  thickness)  may  be  determined  for  all  layers  of  interest  from  the 
pressure  height  curve. 

In  addition  to  figure  3.8,  a  larger,  full-size  overlay  {for  use  with 
the  large- size  Arowagram)  has  been  printed  on  transparent  paper  and 
is  provided  as  a  loose-leaf  inclusion  to  this  manual  (chart  V).  As  in 
the  case  of  charts  III  and  IV  it  is  recommended  that  working  copies  of 
chart  V  be  made  on  transparent  plastic  stock. 
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4.  THE  STANDARD  ATMOSPHERE 


in  the  preceeding  chapter  we  have  defined  a  “standard  refrac- 
tivity  profile”  in  a  "standard  atmosphere”  and  have  disc... ssed  the  com¬ 
putation  of  N- gradient  va-ues  in  both  the  standard  and  in  the  actual 
atmosphere.  In  this  chapter  we  shall  discuss  the  effect  of  refractive 
gradients  in  the  standard  atmosphere  upc.  ->ropagation  and  shall  des¬ 
cribe  the  various  "modified”  refractive  indices  in  current  use. 

4,1  Refraction  in  a  Standard  Atmosphere 


We  have  seen  that  N  decreases  exponentially  with  height  in  the 
standard  atmosphere  and  that,  as  a  consequence,  electromagnetic  rayt. 
are  bent  downward  by  refraction  in  such  an  atmosphere.  Figure  4.1 
shows  a  radio  frequency  source  of  electromagnetic  radiation  above  the 
surface  of  the  earth.  The  lit  zone  (or  interference  zone)  forms  the  locus 
of  aU  direct  paths,  while  the  crosshatched  area  below  the  grazing  path 
is  the  diffraction  region.  Note  that  the  radio  horizon,  C,  lies  beyond 
the  geometric  horizon.  A,  and  the  optical  horizon,  B,  In  a  standard 
atmosphere  the  curvature  of  a  radio  ray  is  about  one-fourth  the  curva¬ 
ture  of  the  earth  so  that  the  radio  horizon  extends  beyond  the  geometric 
horizon  by  about  15  percent.  Tables  have  been  calculated  for  radio  hori¬ 
zon  distaii  :e  versus  transmitter  height  in  a  standard  atmosphere  for 
several  values  of  surface  refractivity,  Ng.  The  results  of  such  a  cal¬ 
culation  are  plotted  in  graphical  form  in  figure  4.2. 

4,2  The  Four-Third’s  Earth  Radius  Concept 


In  the  first  few  thousand  feet  above  the  surface,  the  exponential 
decrease  of  N  with  height  (in  the  standard  atmosphere)  may  be  closely 
approximated  by  a  linear  decrease.  The  upper  sketch  in  figure  4,3 
shows  the  earth’s  E.u.ria  -  and  a  langcnl  ray  having  a  curvature  equal 
to  one-fourth  the  earth’s  curvature.  At  a  givq;'.  distance  t.eyond  the 
tangent  point,  a  straight  tangent  line  will  be  above  the  surface,  xH 
we  neglect  second  order  tern.s  (see  fig,  2,2).*’  At  the  same  distance 
away,  the  straight  line  will  be  above  the  curved  ray.  Assuming 

is  constant  and  equal  to-^,  it  follows  that  the  ray  will  be 
'»bove  the  earth  at  the  distance,  d.  The  same  height-distance  relatior- 
lip  li.'twten  rays’ll!  sjuiface  maybe  obtained  if  one  makes  the  ray 
straight  and  modifies  the  earth’s  radius  to  a  new  value  Bq.  such  that 
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Figure  4.3.  Comparative  Ray  Geometry. 
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Thus,  by  making  the  “effective"  earth’s  rc  lius  four-thirus  as  great  as 
the  true  radius,  the  situation  has  been  restOL  ed  to  that  of  the  homoge¬ 
neous  atmospheric  case  (straight  ray  and  curved  earth)--see  middle 
sketch  in  figure  4.3.  This  method,  ad-vanced  by  Schelling,  E arrows,  and 
Ferrel  [19],  assumes  an  earth  suitably  larger  than  the  actual  earth  and 
allows  radio  rays  to  be  drawn  as  straight  lines  over  this  earth  rather 
than  as  curved  rays  over  an  earth  with  true  radius  a. 


In  general,  the  effective  earth’s  radius  a^  =  ka,  where  k  is  the 
effective  earth’s  radius  factor  and  a  is  the  true  radius  of  the  earth. 
The  value  of  k  is  usually  assumed  to  be  four-thirds  in  the  United  States. 
It  is  desirable,  however,  in  some  geograpMc  areas  to  use  a  somewhat 
different  value  of  k,  for  example  in  desert  areas  it  is  more  nearly  six- 
fifths.  One  may  write  in  general 


k  = 


1  +  a 


An 

Ah 


(4.2) 


where  is  the  assumed  linear  gradient  of  the  index  of  refraction 
(usually  negative),  applying  to  a  particular  area  or  season  and  a  is  the 
sea  level  radius  of  the  earth. 


This  method  of  accoimting  for  atmospheric  refraction  permits  a 
tremendous  simplification  in  the  computation  of  radio  field  strengths 
even  though  the  distribution  of  refractive  index  (Linear  with  height)  im¬ 
plied  by  this  method  *?clistic  only  in  th^  lowest  few  thousand  feet  of 
the  atmosphere. 

4.3  The  Flat  Earth  Concept 


For  some  purposes,  it  is  desirable  to  portray  the  earth  as  beinj,' 
flat  with  upward  curving  rays.  This  is  done  by  subtracting  the  earth’s 
cnrvo.  uxo  frnm  both  raj  r  earth.  This  results  in  an  upward  ray  cur¬ 
vature  equal  to  the  downward  curvature  given  to  the  eai  th  in  the  situa¬ 
tion  described  in  section  4.2  (see  bottoms  etch  in  fig.  4.3).  This  rep- 


-  62  - 


resentation  is  not  now  widely  used  because  of  its  unrealistic  portrayal 
of  the  earth-ray  relationsiiip.  However,  there  are  instances,  such  as 
ray  plotting  by  analogue  computers,  where  it  is  advantageous  to  use  the 
“flat  earth”  type  of  presentation, 

4,4  Modified  Refractive  Indices 


The  B_“,  M-,  A-,  and  K-modifications  of  i.  efractivity  are  useful, 
particularly  in  ray  tracing  problems  and  in  synoptic  studies  of  refrac- 
tivity.  Each  results  in  simplifications  in  ray  geometry  in  certain  situa¬ 
tions,  not  only  in  plotting  raytrajectories  but  also  in  the  mathematical 
relationships  describing  various  propagation  mechanisn^s  and  in  the 
study  of  climatic  and  synoptic  relationships.  There  is  much  to  be  said, 
however,  for  retaining  the  unmodified  r  efractivity  N  itself.  It  is  a  basic 
physical  quantity,  and  as  such  is  more  easily  understood  than  are  B, 
M,  A,  or  K.  For  the  purposes  of  uniformity  and  initial  simplicity,  N- 
units  are  therefore  the  principle  refractive  index  units  employed  in  this 
manual  (except  in  part  IV).  For  reference  purposes, however,  a  des¬ 
cription  of  each  of  the  modified  indices  is  given  in  the  following  four 
subparagraphs, 

4.4.1  B-Units 


It  was  shown  in  section  4,2  above,  that  one  can  compensate  for  a 
linear  refractive  index  gradient  merely  by  changing  the  effective  radius 
of  the  earth.  The  most  convenient  way  to  represent  departures  of  a 
given  refractivity-profile  from  the  standard  atmosphere  (at  least  in  the 
lowest  layers  where  the  departure  of  the  N-gradient  from  the  assumed 
constant  value  's  negligible)  is  to  add  12  N-units  per  1,000  feet  to  all 
N-values,  This  modified  index,  B,  is  defined  as 

^  =  Nj,  0.012  h  (4.3) 


where  ^  is  the  value  of  refractivity  at  any  height  h  in  feet.  A  B-pro- 
filc  will  have  zero  gradient  in  the  lower  layers  of  a  standard  atmos¬ 
phere.  The  B-modification,  therefore,  is  a  logical  consequence  of  the 
four-third  earth  radius  concept  described  in  section  4.2  because  the 
B-gradient  is  zero  for  the  same  situation  where  the  ray  curvature  is 
zero  A  .-'oii.dar<^  otrnvu'phei '  a-profile  is  compared  with  standard  N-, 
A-,  a»id  M- profiles  in  figures  4.4  and  4.5.  Note  that  in  the  lowest  few 
thousand  feet_B  has  a  zero  gradient  but  at  higi  .r  levels  (where  the  de¬ 
parture  from  an  assumed  constant  surface  value  of  minus  12  N- 
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REFRACTIVITY  N 

Figure  4.4,  N-,  B-,  and  M-Profiles  for  a  Standard  Atmosphere. 


Figure  4.5.  N-,  A-,  B-,  and  M-Profiles  for  a  Standard  Atmosphere. 
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units  per  1,000  feet  has  become  appreciable)  B  shows  a  small  positive 
gradient. 


4.4.2  M-Units 

Another  form  of  modified  refractive  index  has  been  widely  used, 
and  is  associated  with  the  flat  earth  concept  of  . '  :;tion4,3.  In  this  case, 
one  adds  48N-units  per  1,000  feet  to  all  N- values.  This  modified  index, 
M,  is  defined  as 


M  =  +  0,048 h  ,  (4.4) 

where  is  the  value  of  refractivity  at  any  height  h  in  feet.  When  the 
M-gradient  is  zero  the  ray  curvature  is  zero  in  the  flat  earth  case. 
This  is  another  v/ay  of  saying  that  when  the  N- gradient  is  minus  48  units 
per  1,000  feet,  the  ray  has  the  same  curvature  as  the  earth.  Thus,  if 
the  earth’s  surface  is  represented  as  a  straight  line,  the  ray  will  also 
be  straight.  Figures  4,4  and  4.5  show  a  standard  atmosphere  profile, 
expressed  in  terms  of  N-,  A-,  B-,  and  M-units  for  comparison  pur¬ 
poses. 


4.4,3  A-Units 

As  we  have  seen  the  B-unit  modified  index  overccrrects  for  the 
decrease  in  N  with  height  in  the  standard  atmosphere  as  shown  in  figure 
4.5  by  the  increase  of  ^  with  height.  This  overcorrection  results  from 
the  assumption  of  a  consteint  N-gradient  (i.e,,  linear  N-profile)  with  a 
value  representative  of  the  lowest  few  thousand  feet  of  the  standard 
atmosphere.  A  modified  refractive  index  (potential  refractivity)  ex¬ 
pressed  in  A-units  has  been  defined  by  Bean,  Riggs,  and  Horn  [4]  in 
order  tc  eliminate  +bis  overcorrecticn.  The  expression  for  A,  which 
is  written  in  terms  cv  an  expontiuidi  height  lu.iction,  is  given  by 

A  =  Nh  +  [l  -  axp.='»]  , 


or 

A  =  +  313  [l  -  exp.“^- (4.5) 

where  is  the  'efractivity  at  height,  h,  in  kilometers.  The  decay  co¬ 
efficient  c  =  -  0.144  was  determined  so  as  to  correspond  best  with  the 
average  station  value  of  Ng  =  313  for  the  United  States.  Figure  4.5 
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shows  a  comparison  between  the  staiidard  atmosphere  profile  expressed 
in  terms  of  N-,  A-,  B-,  and  M-units^ 

To  facilitate  the  transformation  fromN-units  to  '  •units,  a  poten¬ 
tial  refractivity  chart  was  prepared  and  is  given  in  figure  4,6,  This 
chart  eliminates  the  necessity  ©fusing  exponential  tables  in  each  indi¬ 
vidual  calculation  of  A  and  thus  lends  considerable  ease  to  the  prepara¬ 
tion  of  charts  of  the  new  parameter. 

4.4.4  K-Units 

In  order  to  remove  the  height-dependent  variations  of  refractivity 
in  homogeneous  air  (the  problem  for  which  Bean  proposed  the  A-unit 
modification)  and  to  develop  a  modified  index  from  fundamental  meteor¬ 
ological  relationships  (rather  than  from  essentially  empirical  consid¬ 
erations  as  Bean  has  done),  a  modified  index  (potential  refractive  index) 
expressed  inK-imits  has  been  defined  byJehn  [16].  The  expression  for 
K  is  given  by 


where  is  the  conventional  "dry"  term  and  the  conventional 
"wet"  term  in  the  refractivity  equation  and  P  is  the  pressure  in  milli¬ 
bars  at  the  point  in  question.  Tables  or  graphs  of  the  two  multipliers 
could  be  prepared  to  facilitate  the  conversion  from  N-units  to  K-units, 
but  as  yet  these  are  not  available. 

It  should  be  noted  that  K- values  are  referred  to  the  1,000  millibar 
pressure  surface,  and  as  a  result  the  station-to- station  variability  of 
refractivity  occas’rtf-.  ’  by  differences  in  station  elevation  are  strongly 
suppressed.  In  synoptic  radio- meteorology  this  becomes  sin  important 
factor. 


A  =  N  +  NsD-e‘^^3 

WHERE  Ns=3l3 

c  =-0.144 
h  =  height  in  Km 
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5.  THE  ACTUAL  ATMOSPHERE 


Thus  far,  the  discussion  has  been  concerned  witi,  '’e  propagation 
of  radio  waves  under  two  artificial  situations;  an  earth  enveloped  by  a 
homogeneous  atmosphere,  and  one  with  an  atmosphere  in  which  refrac- 
tivity  decreases  smoothly  (linearly  or  expo'  entially)  with  height.  As  a 
final  step,  it  is  now  in  order  to  consider  the  actual  atmosphere  and  its 
effects  on  propagation. 

There  are  numerous  ways  in  which  the  troposphere  departs  from 
the  simple  linear  or  exponential  gradient  model.  In  a  well  mixed  tro¬ 
posphere,  for  example,  there  are  a  myriad  of  small  random  variations 
in  N,  both  in  vertical  and  horizontal  planes.  Even  at  a  fixed  point,  N 
is  found  to  vary  in  a  random  way  about  some  mean  value.  These  fluc¬ 
tuations  are  caused  by  parcels  of  air  with  different  temperatures  and/ 
or  vapor  pressures  than  their  surroundings.  Most  parcels  are  caused 
by  unequal  heating,  during  the  daytime,  of  air  in  contact  with  the  ground. 
The  warmer  parcels  rise  convectively  and  drift  with  the  wind  until  they 
are  broken  up  and  eventually  lose  their  identity  through  turbulent  mixing 
with  the  surrounding  air.  Similarly,  air  overlying  lakes,  rivers,  and 
oceans  becomes  more  moist  than  its  surroundings  and  rises  because 
of  its  lower  density.  Refractometer  measurements*  show  that  the  root 
mean  square  value  of  the  N  fluctuations  decreases  exponentially  with 
height  in  a  well  mixed  atmosphere.  This  is  a  priori  evidence  that  the 
source  of  the  small-scale  variability  of  N  is  at  the  surface, 

5.1  Anomalous  Propagation 

In  some  cases  the  actual  refractive  index  profile  in  the  atmosphere 
departs  from  the  linear  and  exponential  models  enough  to  cause  radio 
and  radar  fields  to  ’  '•h  etronper  well  boyond  the  horizon  than  would 

be  expected  from  standard  refraction  and  diffractioii.  Thece  situations 
are  generally  caused  by  the  presence  of  steep,  negative,  vertical  gra¬ 
dients  of  refractivity  in  layers,  which  in  turn  cause  the  radio  waves  to 
be  bent  downward  sufficiently  to  continue  some  distance  on  around  the 
figure  of  the  earth.  A  discussion  follows  of  the  development  and  of  the 
characteristics  of  these  so-called  "trapping”  or  "ducting”  layers. 

5.1.1  Surface  Layers 

The  actual  troposphere  departs  most  .frequently  from  standard  in 
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the  lowest  fewhimdred  feet.  These  surface  layers,  in  which  the  N- gra¬ 
dient  is  much  steeper  than  minus  12  units  per  1,000  feet,  are  caused 
by  nocturnal  radiational  cooling,  by  evaporation  or  by  advection.  Fre¬ 
quently  the  N-gradient,  particularly  just  above  the  ground,  ■  consider¬ 
ably  stronger  than  minus  48  uni:s  per  1,000  feet,  such  that  a  ray  actually 
has  a  downward  curvature  greater  than  the  earth’s  curvature.  If  the 
layer  is  thick  enough  compared  to  the  propag  tion  wavelength,  energy 
is  guided  around  the  figure  of  the  earth  in  thela>er,  causing  fields  be¬ 
yond  the  horizon  to  be  very  much  stronger  than  normal.  Because  of  the 
wave  nature  of  radio  energy,  it  turns  out  that  the  layer  thickness  must 
exceed  a  certain  critical  value,  which  is  a  function  of  wavelength  and 
of  the  shape  of  the  N-profile,  before  guiding  or  “ducting”  of  the  wave 
can  occur. 

When  the  layer  becomes  particularly  intense,  the  horizontally 
directed  energy  from  a  transmitter  located  within  the  layer  will  be 
completely  contained  in  the  layer;  in  this  instance  the  fields  will  de¬ 
crease  as  the  square  root  of  distance  rather  than  as  the  first  power  of 
distance.  At  some  distance  beyond  the  horizon,  these  fields  would  then 
theoretically  not  only  be  stronger  than  the  diffracted  fields,  but  would 
even  be  stronger  than  the  free  space  fields.  Such  intense  layers  have 
not  actually  been  found  in  the  atmosphere.  The  strongest  fields  observed 
in  surface  ducts  have  been  about  equal  to  the  free  space  fields.  When 
the  fields  approach  this  level,  the  energy  is  said  to  be  “trapped”.  The 
refractivity  profile  requirements  for  trapping  radio  energy  of  various 
wavelengths  are  described  in  more  detail  in  section  6.1. 

(a)  Radiation  Ducts 

Nocturnal  radiational  cooling  of  the  grovmd  surface  may  produce 
“radiation”  ducts.  They  are  particularly  intense  over  dry  or  desert 
areas,  where  there  is  v/rter  vripnr  in  the  air  to  intercept  and  re¬ 
turnoutgoing,  longwave  radiation  back  to  the  ground.  Calm  nights  tend 
to  intensify  duct  formation  because  there  is  little  mechanical  turbulence 
to  mix  the  cool  lower  layers  with  *he  warmer  air  above.  The  typical 
N-profile  is  exponential  in  shape,  merging  v/ith  the  standard  gradient 
at  heights  of  a  few  hundred  feet.  Measurements  in  southwestern  Arizona 
by  Day  and  Trolese  [10]  showed  that  usually  the  profile  is  standard  at 
sunscl,  ir.ie  suxface  N  '  -.iue  increasing  as  the  night  progresses. 

At  siuirise,  the  surface  N-’/alue  has  typically  increased  by  about  15 
units  v^hile  the  values  above  500  feet  haveincre.  sed  by  only  2  or  3  units. 
Usually  *he  lop  of  the  duct  is  arbitrarily  taker  to  be  the  height  where 
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the  gradient  is  minus  48  units  per  1,000  feet.  A  typical  value  of  this 
height  at  sunrise  in  Arizona  is  about  50  feet.  S\ich  a  duct  would  be  ex¬ 
pected  to  trap  wavelengths  shorter  than  about  30  centimeters,  and  ex¬ 
perimental  results  [10]  show  this  to  have  been  the  case. 

(b)  Evaporation  DuctJ 

Over  the  ocean,  surface  refracting  l../  rs  are  caused  by  a  dif¬ 
ferent  process  than  that  described  in  (a)  above.  Diurnal  temperature 
changes  over  large  bodies  of  water  arc  usually  quite  small.  This  is 
because  solar  energy  penetrates  well  beneath  the  sea  surface.  There¬ 
fore,  a  relatively  large  amount  of  heating  is  required  to  change  the  sea 
surface  temperature,  Radiational  cooling,  which  does  occur,  makes 
the  surface  water  heavier  than  the  under Ijdng  water,  so  it  sinks,  bring¬ 
ing  up  warmer  water  from  below.  The  surface  temperature,  therefore, 
remains  practically  constant  night  and  day.  Daytime  heating  can  cause 
a  thin  surface  layer  of  higher  temperature,  but  wave  action  usually 
rapidly  mixes  this  layer  with  the  cooler  water  below. 

Evaporation  from  the  sea  surface,  on  the  other  hand,  causes  air 
parcels  to  be  carried  aloft  by  buoyancy,  and  thus  gives  rise  to  a  steep 
gradient  of  water  vapor  pressure  (and  hence  of  N)  in  the  first  few  hxm- 
dred  feet  above  the  sea.  While  instantaneous  N-profiles  show  wide  fluc¬ 
tuations,  the  average  over  several  minutes  shows  a  profile  which  is  ex¬ 
ponential  in  shape  when  the  air  and  sea  are  the  same  temperature.  When 
the  air  is  warmer  than  the  sea  the  profile  is  more  nearly  linear;  i.e., 
the  N-gradient  changes  more  slowly  with  height.  In  the  unstable  case 
(sea  warmer  than  air),  the  gradient  decreases  more  rapidly  with  height 
than  in  the  neutral  case  (sea  and  air  temperature  the  same).  An  in¬ 
teresting  consequence  of  the  physical  processes  which  take  place  in  the 
formation  of  an  evaporation  duct,  is  that  in  the  stable  case,  low  wind 
speeds  favor  strong  duct  formation;  while  in  the  unstable  case,  high 
wind  speeds  favor  r-iroiig  ducts.  These  hyi  -'theses  have  been  verified 
in  data  taken  over  the  Irish  Sea  by  Anderson  and  Ccss'»rd  [1].  Criteria 
are  presented  in  section  6,1,3  which  enable  one  to  determine  from 
meteorological  measurements  w,c  bridge  elevation,  whether  or  not  a 
given  wavelength  will  be  trapped  by  the  oceanic  duct. 

(c)  Advection  Ducts 

Another  process  by  which  surface  layers  with  superstandard  N- 
gradients  can  be  formed,  is  the  advection  of  dry  air  over  the  sea.  If 
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the  air  is  colder  than  the  sea .  vertical  convection  takes  place  whii;h 
tends  to  disperse  air  parcels  with  high  refractivity  upward,  thus  pro¬ 
ducing  a  weak  but  thick  surface  layer  which  does  net  usually  affect  prop¬ 
agation  appreciably. 

However,  if  the  air  ic  warmer  than  the  sea,  stability  will  keep  the 
cool,  moist  (high  air  confined  to  a  tliiiuier  layer.  This  gives  rise  to 
surface  layers  in  which  the  N-gradient  can  be .  ong  enough  for  "trap¬ 
ping"  and  at  the  same  time  be  thick  enough  to  affect  longer  wavelength 
energy  than  is  the  case  with  e\  aporation  ducts.  This  situation  is  prev¬ 
alent  over  the  Mediterranean,  Red  Sea,  and  Gulf  of  Arabia. 

5.1,2  Elevated  Layers 

Thus  far,  processes  which  cause  persistent  departures  of  N  in 
atmospheric  layers  cl:)se  to  the  earth  or  sea  surface  have  been  con¬ 
sidered.  These  layers  are  caused  by  processes  taking  place  at  the  sur¬ 
face  itself.  There  is  another  class  of  phenomena  which  gives  rise  to 
strong  N-gradients  at  some  height  above  the  surface.  Refractomeler 
studies  have  shown  that  these  elevated  layers  exist  to  some  extent  over 
wide  geographic  areas.  In  fact,  very  few  refractometer  profiles  have 
been  taken  which  do  not  show  the  presence  of  at  least  one  elevated  layer 
within  10,000  feet  of  the  earth’s  surface. 

Most  of  these  layers  are  only  a  few  hundred  feet  in  thickness  and 
the  change  of  N  from  top  to  bottom  averages  about  1 0  units.  Some  layers 
can  be  traced  for  tens  of  miles  in  horizontal  extent  while  others,  par¬ 
ticularly  the  weaker  ones,  are  less  than  10  miles  in  extent.  Figure  5.1 
shows  an  N-profile  through  an  idealized  elevated  layer.  Several  para¬ 
meters  of  interest  are  shown;  the  intensity,  ANj^  ,  the  thickness.  Ah  , 
and  the  heights  of  the  top  and  bottom  of  the  layer,  h-p  and  hs.  respec¬ 
tively.  It  will  he  noted  that  the  N-gradient  is  shown  as  being  standard 
below  and  above  the  ’>  iy£’. .  This  lo  aJually  tr...  rzse  although  nominal 
departures  from  standard  are  often  found. 

Figure  5.2,  shows  the  effec.  of  such  an  elevated  layer  on  rays 
emanating  from  a  transmitter  located  well  below  the  layer  bottom.  The 
four-thirds  earth’s  radius  modification  has  been  used  in  constructing 
tso  oq,rth  and  layer  boundaries.  Since  (as  shown  in  fig.  5.1)  the  N-gra- 
dien  rl^nd-.rd  outsir'e  In*,  tayer,  the  ray  paths  above  and  below  the 
layer  are  straight  lines.  In  the  layer  itself,  the  N-gradient  exceeds 
standard,  therefore,  the  rays  bend  downward. 


Figure  5.1.  N- Profile- -Elevated  Layer. 


Using  these  geometric  considerations  and  Snell’s  Law,  and  with 
the  help  of  expx'essions  previously  derived  for  ray  bending  in  a  medium 
of  constant  N-gradient,  one  obtains  the  following  expressions  for  the 
distance,  d,  between  the  radar  and  the  point  of  entry  into  the  layer,  and 
Ad,  the  distance  the  ray  travels  in  the  layer. 


(5.1) 


and 


2AnL  +  2  (Ah/ag)  +  2(hg 


- »(5.2) 


Ad  = 
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where  a„  is  the  effective  earth’s  radius,  a  is  the  inclination  ans'le  of 
the  ray  as  it  leaves  the  transmitter,  hg  and  hg  are  the  heights  of  the:  lay¬ 
er  base  and  the  transmitter  respectively,  ^  is  the  layer  thickness,  and 
An^;  i.e.,  (ANg  x  10"®);  is  the  change  in  the  index  of  refraction  through 
the  layer  (see  figs.  2.6  and  5.2).  The  distance  from  thv  ‘^ransmitter  at 
which  the  ray  emerges  is  simply  the  sum  of  these  two  expressions. 
From  this  point,  the  ray  proceeds  in  a  straight  line  as  shown  in  figure 
5.2.  Figure  5.2  also  illustrates  the  behavx^  of  the  rays  as  the  eleva¬ 
tion  angle,  a^,  is  ehaiiged.  For  high  elevation  rays  (above  2°)  the  bend¬ 
ing  is  slight,  but  as  bec<  mes  smaller,  the  rays  are  deflected  more 
and  more  from  their  original  direction. 

Figure  5.3,  illustrates  the  effect  of  raising  the  transmitter  height 
(or  lowering  the  layer  height).  The  high  angle  rays  still  penetrate  the 
layer  as  before,  but  now,  as  is  decreased,  the  rays  are  bent  more 
and  finally  (in  the  tldrd  ray  from  the  left)  emerge  tangent  to  the  layer 
top.  If  is  decreased  ex^en  further,  the  rays  do  not  penetrate  the  layer 
at  all,  but  are  totally  refracted  back  down  toward  the  surface  again. 
The  reason  for  this  will  be  evident  from  a  consideration  of  the  expres¬ 
sion  for  A^.  If  (hg  -  ^r)  and  ffjj  are  sufficiently  small,  the  first  radical 
becomes  imaginary,  signifying  no  emergent  ray.  When  this  radical  is 
just  equal  to  zero,  the  emergent  ray  is  tangent  to  the  layer  top.  How¬ 
ever,  as  becomes  negative,  there  will  again  be  a  value  of  cc^,  which 
is  larger  than  the  remaining  terms  in  the  radical,  and  rays  will  again 
emerge  from  the  layer  top.  This  situation  is  depicted  in  figure  5.3  by 
the  lowest  ray  in  the  bundle.  Between  the  two  emerging  tangent  rays 
is  a  region  where  no  energy  penetrates  the  layer,  A  radar  target  in 
this  region  would  not  be  detected,  nor  would  communications  be  possi¬ 
ble  between  an  airplane  in  the  region  and  a  surface-located  transmitter. 
This  region  xs  called  a  “radio  or  radar  hole,”  A  discussion  of  methods 
for  locating  these  regions,  as  well  as  avoiding  them,  is  given  in  section 
6.3,1. 


In  figure  5.4,  the  transmitter  is  above  ihc  layer.  In  lias  case, 
the  upper  rays  do  not  intersect  the  layer  at  all,  and  are  unaffected  by 
it.  There  are,  however,  two  rays  which  again  define  a  “radio  hole”  in 
this  situation.  The  first  is  the  ray  tangent  to  the  layer  top,  and  the 
second  is  a  ray  going  down  at  a  steeper  angle,  which  comes  back  up 
+hT'ongh  the  layer  farther  ''ut.  The  distance  from  the  transmitter  to  the 
ta  .'gent  poi"v.  I  start  C'  the  radio  hole)  is  simply 

d»  =x J2a^  (hjj  -  h,,,)  , 


(5.3) 


Figure  5.3.  Ra; 
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where  and  hj^  are  the  height*?  of  the  layer  top  and  transmitter  res¬ 
pectively  and^g  is  the  effective  earth’s  radius.  The  width  of  the  radio 
hole  is  more  complicated  to  express  in  this  case,  than  when  the  trans¬ 
mitter  is  below  the  layer.  In  section  6.3.1,  practical  methods  are  pre¬ 
sented  for  operational  use. 

The  main  difference  between  the  transmittpr  being  below  the  layer 
and  being  above  the  layer,  is  that  in  the  former  x.  se  one  can  eliminate 
the  radio  hole  completely  by  lowering  the  transmitter  far  enough  below 
the  layer.  If  hg  -  hj^  is  large  enough,  all  rays  will  emerge  from  the 
layer  regardless  of  £q.  When  the  transmitter  is  above  the  layer,  on  the 
other  hand,  raising  the  transmitter  moves  the  hole  further  away  and 
gradually  narrows  it  but  can  never  eliminate  it. 

Substituting  10  units  for  AN^:(An£_  =  10x10“®)^  and  200  feet  for 
Ah  in  the  expression  for  Ad,  shows  that  if  the  transmitter  is  more  than 
80  feet  below  the  layer  base,  there  will  be  no  radio  hole.  It  can  be  seen 
then  that  a  typical  elevated  layer  will  have  little  effect  until  the  trans¬ 
mitter  is  nearly  at  layer  height.  For  the  same  layer,  with  the  radar 
80  feet  above  the  layer,  a  radio  hole  would  exist  at  radar  height  between 
22  and  60  miles  from  the  transmitter. 

The  ray  tracing  expressions  discussed  above,  are  theoretically 
valid  only  when  the  layer  is  smooth  and  has  a  uniform  N-profile  o\er 
a  large  horizontal  area  surroimding  the  transmitter.  Even  though  actual 
elevated  layers  rarely  are  as  smooth  as  the  mathematical  model  re¬ 
quires,  the  results  of  many  experimental  flights  above  and  below  such 
layers  give  results  as  indicated  by  Doherty  [11],  which  are  in  surpris- 
ingly  good  agreement  with  theory.  Field  strengths  in  radio  holes  are 
found  to  be  considerably  lower  than  outside  them,  and  the  distances  to 
the  boundaries  agree  very  well  with  those  calculated  from  ray  theory 
and  measured  N-prof:lcs. 


In  certain  geographical  areas,  strong  subsidence  of  dry  aix'  from 
aloft  gives  rise  to  intense  elevated  inversions  (layers)  with  AN’s  of  the 
order  of  40-60  units  and  heights  of  less  than  3,000 feet.  These  areas 
are  near  30®N.  (and  S.)  latitude  in  those  belts  where  atmospheric  sub¬ 
sidence  is  strongest.  In  the  Northern  Hemisphere,  the  areas  of  occur- 
renetj  v-f  ‘ih*’  sfrongesi.  elev'.-i.'d  inversion  (layers)  are  located  where 
the  30  ’i\.  iatituae  belt  h.tersects  the  west  coasts  of  large  continents; 
namely,  southern  California,  northern  Mexico,  ’'nd  Morocco.  Little  data 
is  available  for  the  Southex’n  Hemisphere,  although  such  evidence  as 
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exists  indicates  that  similar  areas  will  be  found  on  the  west  coasts  of 
the  continents  in  that  hemisphere  as  well.  Ducts  of  this  intensity  are 
frequently  able  to  trap  energy  from  surface  radar  and  also  at  communi¬ 
cation  wavelengths. 

In  layers  where  ANj.  is  large,  causing  the  total  refraction  of  rays 
with  relatively  large  roughnesi.  the  layer  becomes  a  very 

important  factor  in  feeding  energy  into  the  radio  holes.  In  experimental 
studies  off  southern  Calilo  'nia,  it  has  been  found  that  in  such  situations, 
the  field  strength  in  the  radio  hole  is  weaker  than  it  is  above  or  below 
the  hole,  but  it  is  stronger  than  when  no  layer  is  present.  This  comes 
about  because  of  the  Rayleigh  criterion  of  roughness,  as  follows: 

h  =  X/8  sin  ,  (5.4) 

where  h  is  the  permissible  height  of  random  roughness  elements  for  a 
reflecting  surface  to  appear  smooth  at  wavelength  X  and  with  the  inci¬ 
dent  angle  equal  to  £q. 

Suppose,  for  example,  that  a  layer  has  the  following  character¬ 
istics: 


and 


AN  =  50  , 


Ah  =  200  feet  . 

We  find  that  total  refraction  will  occur  for  Cq  as  large  as  9.3  milli- 
radians.  In  this  case,  the  roughness  elements  must  be  less  than  13 
wavelengths  in  height  for  the  layer  to  appear  smooth.  On  the  other  hand, 
if  AN  is  10  \mits,  as  it  was  in  our  previous  layer,  ^  can  only  be  as 
large  as  2,4  mill:.!  aa’.anc,  u  perinis :  Ible  roughness  element  height 

of  53  wavelengths.  From  this,  it  is  evident  thnt  r.  s+T’ong  elevrited  layer, 
which  is  capable  of  trapping  energy  incident  at  rather  large  angles., 
must  also  be  much  smoother  .han  is  usually  foimd  in  nature  in  order  to 
give  a  well  defined  radio  hole.  Weaker  layers,  which  require  small  ray 
angles  for  radio  holes  to  form,  can  be  relatively  much  rougher  and  si  ii 
produce  well  defined  holes. 

5,2  Ray  Bending 

The  previous  section  considered  the  behavior  of  rays  in  the  pres- 
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ence  of  steep  N-gradients.  Of  primary  concern  were  those  rays  wMch 
were  bent  sufficiently  to  emerge  tangent  to  the  top  of  the  layer.  Now 
we  wish  to  consider  situations  where  the  bending  is  not  so  drastic.  It 
should  be  realized,  of  course,  that  such  bending  occurs  m  any  atmos¬ 
phere  where  N  changes  with  height.  If  N  decreases  with  .  eight,  as  it 
does  over  the  earth,  the  rays  are  bent  downward.  As  was  seen  ir  sec- 
tion  tills  l?0ndinrr»  t,  3.t  s.Jiy  point  rsy  tiith— 

portional  to  the  N-gradient  normal  to  the  ray  ptch. 

In  order  to  determine  t  exactly,  one  must  know  how  N  varies  with 
height.  In  general  this  relationship  esmnot  be  expressed  in  simple  ana¬ 
lytical  form.  Furthermore,  the  N-profile  changes  from  place  to  place 
and  from  time  to  time.  In  practice,  it  is  more  satisfactory  to  approxi¬ 
mate  a  given  N-profile  by  a  series  of  straight  lines  and  then  to  deter¬ 
mine  the  incremental  bending  caused  by  each  segment.  The  total  bend¬ 
ing  up  to  a  given  height  is  simply  the  sum  of  the  individual  increments 
up  to  that  height.  This  procedure  is  particularly  suited  to  N-profiles 
derived  from  radiosonde  data,  since  significant  levels  are  originally 
specified  using  linear  interpolation  between  levels  to  recover  the  orig¬ 
inal  trace. 


Using  this  approach,  one  may  express  the  incremental  bending 
At  caused  by  a  layer  with  constant  N-gradient  as. 


At  (milliradians)  = 


(Ng  -  Niji) 
500  (tan  Pg  +  tan 


(5.5) 


where  the  subscript  B  refers  to  the  bottom  of  the  layer  and  T  refers  to 
the  top  of  the  ?ayer  (see  fig,  2,6).  The  values  of  tan  p  are  determined 
at  each  level  in  terms  of  a^,  Ng  -  Nj,,  and  h,  using  Snell’s  Law.  Figure 
5.5  gives  values  of  500  tan  ^  in  terms  of  these  parameters. 


5,2,1  Elevation  Angie  Errors 


In  tracking  targets  with  radars  capable  of  measuring  elevation 
angles  very  accurately  one  is  more  interested  in  correcting  the  apparent 
elevation  angle,  a^,  to  the  true  angle,  than  in  the  value  of  t  itself.  Fig¬ 
ure  2.6,  shows  this  correction  or  error  angle,  6.  If  one  knows  the  total 
benouig  belweer  radar  h-'ight  and  target  height,  h^,  then, 

7^  tan  4{-^  j-  (Nf  -  Nt) 

6+ (milliradians)  - - .  (5.6) 

T^  +  tan  -  tai:  af^ 


-  79  - 


IIJ-il.lLILlLiMliliiii' 

GRAPHIC  REPRESENTATION  OF 
SNELL'S  LAW  FOR  FINDING 
500  T."'*  fi 

h  In  thousonds  of  feet 


1000 


Figure  5.5.  Graphical  Representation  of  Snell's  Law. 
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For  low  target  heights,  where  the  alxausphere  may  be  regarded  as  a 
single  layer  of  constant  N-gradient,  Q  is  exactly  haK  of  t.  At  great 
heights  ^  approaches  the  value  of  t.  These  two  limits  are  frequently 
useful  in  estimating  ^’s  without  the  necessity  of  carryir'f  out  detailed 
computations. 

Figure  5.6  shows  curves  of  both  rand  ^  "or  various  a^’s,  over  a 
height  range  of  1,000  to  1,000,000  feet.  These  /alues  were  computed 
assuming  Ng  =  350,  and  an  N-profile  as  shown  in  figure  3.5.  It  can  be 
seen  that  a  ray  departing  tangentially  from  the  earth  ( =  0“)  is  bent 
almost  1“  in  traversing  the  entire  troposphere.  At  an  elevation  angle 
of  30“,  the  tropospheric  bending  amounts  of  2’  of  arc.  It  is  obvious  that 
the  bending  falls  off  rapidly  as  the  elevation  angle  increases.  Even  at 
1“,  the  bending  is  only  two- thirds  of  the  tangential  value,  and  at  5®  it  is 
one-fourth  of  the  tangential  value.  It  is  interesting  to  note  that  at  low 
angles,  most  of  the  bending  takes  place  at  low  altitudes,  while  at  higher 
angles,  the  contributions  at  greater  heights  are  relatively  more  impor¬ 
tant.  To  a  first  order  of  approximation,  one  may  regard  total  bending 
as  proportional  to  N^;  this  is  quite  true  at  angles  above  10°,  At  very 
low  angles,  t  varies  more  rapidly  with  respect  to  Ng,  but  the  linear 
approximation  is  still  reasonably  good. 

The  values  shown  in  figure  5.6  refer  to  a  standard  atmosphere, 
with  an  exponential  decrease  of  N  with  height.  When  layers  exist  which 
have  steeper  N-gradients,  additional  bending  will  occur.  Layers  close 
to  the  surface  are  much  more  effective  in  increasing  r  and  ^  than  are 
those  at  heights  of  several  thousand  feet.  The  additional  bending  caused 
by  such  layers  is  confined  largely  to  the  lowest  few  degrees.  Above  5°, 
even  strapping  layer  will  cause  no  additional  bending.  This  is  because 
the  distance  the  ray  travels  in  the  layer  falls  off  rapidly  with  increasing 
elevation  angle.  The  bending  caused  by  any  given  profile  of  refractive 
index  maybe  computed  «  .  '’.fimction  of  the  ray  elevation  angle  at  the 
radar,  using  the  method  described  above.  The  additional  ^  <'aused  by 
several  typical  layers  is  discussed  in  section  6,2.2  and  is  shown  in 
figure  6.11. 

It  should  be  pointed  out  that  the  above  discussion  of  radio  ray 
bending  has  been  confined  to  refraction  caused  by  the  troposphere.  When 
therndaj: ’.o.  get  in  orabc^'c  the  ionosphere  (heights  above  100  miles), 
additional  bending  is  caused  by  ionospheric  layers.  Total  bending  by 
the  ionosphere  is  roughly  equal  to  the  tropospheric  coxitribution  at  100 
megacycles  (  X  -  3  meters).  At  shorter  wavelengths  the  ionospheric 


Figure  5.6.  Total  Bending  and  Elevation  Angle  Error  for  a  Standard  Atmosphere. 
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bending  decreases  as  the  square  of  wavelength. 

5,2.2  Range  Errors 

Because  radar  waves  travel  slower  in  the  atmosj^.  're  than  in  a 
vacuum  and  since  radar  range  calibrations  are  made  using  che  vacuum 
velocity  as  a  reference,  measured  ranges  are  •« lightly  larger  than  actual. 
These  range  errors  are  for  most  purposes  entii  ly  negligible.  A  method 
for  determining  range  errors  is  given,  however,  in  section  6.2.3. 


THE  COMPUTATION  OF  ATMOSPHERIC  REFRACTIVE 
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8.  RADAR  PERFORMANCE 


In  parts  I  and  11  of  this  manual  the  various  mechanisms  were  dis¬ 
cussed  by  which  radar  or  radio  energy  is  propagated  from  a  transmitter 
to  a  distant  point  in  space.  Also  considered  were  the  >.  'veral  ways  in 
which  meteorological  variables  affect  these  propagation  mechanisms. 
It  v/as  seen  that  relatively  simple  ray-trz  -ing  concepts  are  useful  in 
describing  the  behavior  of  radio  energy  in  various  situations.  In  this 
chapter,  these  concepts  are  put  to  practical  use  in  order  to  solve  the 
operational  problem  of  determining  the  variations  in  radar  performance 
caused  by  meteorological  variables. 

6.1  Surface-to-Surface  Situations 


6.1.1  “Normal"  Detection 

A  surface  radar,  such  as  those  on  ships  or  low  shore  installa¬ 
tions,  normally  is  limited  to  detecting  targets  above  the  horizon  line 
(see  fig.  4.1).  Its  antenna  pattern  is  characterized  by  interference  max¬ 
ima  and  minima  as  shown  in  figure  6.1.  Inside  the  fingers  the  fields 
are  above  the  free  space  fields,  because  the  direct  and  reflected  waves 
reinforce  one  another,  and  in  the  regions  between  fingers  the  fields  are 
less  than  free  space  fields,  because  the  two  waves  partially  cancel  one 
another.  For  low  radar  heights,  the  vertical  angles  of  the  center  lines 
of  maximum  strength  are  given  by 

,  ,  ,  1,000  \ 

=  (n  -  i)  — — -  ,  (6.1) 

where  n  is  1,  2,  3,  etc,,  hj^  is  the  radar  height,  and  >  is  the  wavelength, 
both  expressed  in  the  same  units.  Thus,  for  a  radar  height  of  60  meters 
and  a  wavelength  •.n-^ters.  the  interference  maxima  v/ill  occur  at 

2.5,  7.5,  12.5,  etc.,  milliradians.  The  minima  occur  at  5.  10,  15,  etc., 
milliradians.  The  range  at  which  the  tips  of  the  fingers  occur  is  take* 
arbitrarily  as  the  maximum  d<=‘tection  range  of  the  radar;  this  range 
depends,  of  course,  on  transmitted  power,  size  of  target,  and  other 
factors  not  related  to  propagation  itself. 

.u  xS  e’.-idaiit  Lsat  -....lace  targets  will  not  be  detected  unless  they 
ai  e  high  enough  to  protrude  into  the  lowest  maximum  of  the  pattern. 
Tills  will  occur  under  normal  atmospheric  ondition?  at,  or  ju‘:t  beyond. 


Figure  6.1.  Vertical  Coverage  Pattern — Surface  Radar. 
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the  horizon  range  for  most  shipboard  radar  installations  and  most  sur¬ 
face  targets.  Other  things  being  equal,  a  short-wavelength  radar  will 
detect  surface  targets  farther  away  because  its  lowest  interference  lobe 
lies  closer  to  the  horizou  than  that  of  a  longer-wavelength  radar.  In 
the  presence  of  strong  surface  or  elevated  ducts,  the  lov/er  finger  of 
the  pattern  maybe  bent  do  vnward  considerably  as  sho>. '  by  the  dashed 
pattern  in  figure  6.1.  The  lower  side  of  the  bottom  finger  extends  well 
beyond  the  horizon  in  this  case. 


6.1.2  Evaporation  Ducts  (Surface) 


Evaporation  or  “oceanic”  ducts  are  not  particularly  effective  in 
extending  the  coverage  of  most  shipboard  radars  because  they  are 
usually  shallow  (75  ft.  is  about  average).  In  order  for  such  a  duct  to  be 
effective,  the  radar  must  be  within  50  feet  of  the  sea  surface  and  its 
wavelength  must  be  shorter  than  30  centimeters.  Submarine  radars  of 
3-centimeter  (X-band)  and  10-eentime.ter(S-band)  wavelengths  as  well 
as  similar  radars  on  small  craft  are  the  only  ones  which  are  usually 
affected  by  the  evaporation  duct.  The  normal  practice  of  installing 
radars  on  the  highest  point  on  the  ship  places  most  radars  above  evap¬ 
oration  duets  which  may  be  present.  For  low-sited  radars,  however, 
the  duct  is  often  useful  hi  extending  the  normal  limit  of  coverage.  It 
turns  out  that  one  can  determine,  from  simple  meteorological  data, 
whether  such  a  radar  at  a  particular  wavelength  will  have  extended 
coverage  or  not. 

The  Ductogram,  figure  6.2,  was  designed  at  Bureau  of  Aeronautics 
Project  AROWA  [25]  to  facilitate  the  prediction  of  extended  radar  cover¬ 
age  associated  with  an  oceanic  (evaporation)  duct.  Use  of  the  Ducto¬ 
gram  is  limited  to  low-level  oceanic  ducts.  For  input  data,  the  Ducto¬ 
gram  utilizes  surface  observations  of  air,  dew  point,  and  sea  water 
temperatures  taken  onboard  the  ship.  The  quantity  AN*  is  the  approxi¬ 
mate  difference  the  r^^^activity.  N,  at  the  elevation  of  ship¬ 

board  observation  (bridge  height)  and  the  rstVautivity,  Ng,  at  the  water 
surface. 

In  an  air  mass  which  has  stagnated  over  an  ocean  surface  for 
some  time,  the  values  of  N  at  observation  height  (determined  from  a'r 
and  dew  point  temperature)  and  Ng  at  the  ocean  surface  (determined 

*ho  '■-ca  water  ter;:p,^rature)  normally  do  not  differ  by  more  than 
one  or  two  units.  However,  when  the  meteorological  conditions  favor 
extended  surface  rada  r  propagation,  the  4,radient  of  N  may  reach  much 


DUCTOGRAM 
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Figure  6.2.  Dustogram. 
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larger  values  and  is  especially  pronounced  in  the  lowest  fe  w  feet  of  the 
atmosphere.  Thus,  whei.  AN*  is  large,  the  gradient  of  N  for  a  depth  of 
several  himdredsof  feet  will  also  be  large.  For  this  reason,  AN*  is  in¬ 
troduced  as  a  parameter  to  be  used  in  estimating  the  presence  of  an 
oceanic  duct.  The  magnitude  of  AN*  necessary  for  e'^ended  radar  ranges 
varies  with  the  wavelength  of  the  radar,  ranging  frv,  67  units  for  the 
L-band,  to  36  for  the  S-band,  and  16  for  the  X-band„  For  longer  wave¬ 
lengths,  such  as  the  P-band  {100-500  m  .),the  oceanic  duct  will  be  in¬ 
effective,  and  one  must  use  the  methods  pi  esented  in  the  following  sec¬ 
tion  to  determine  whether  trapping  will  occur.  Since  air  turbulence 
affects  the  gradients  of  temperature  and  moisture  throughout  the  layer 
in  which  surface  ducting  occurs,  the  surface  wind  velocity  is  introduced 
as  a  factor  in  the  prediction  of  surface  ducting. 

The  procedure  to  be  used  with  the  Ductogram  is  as  follows:  given 
the  air  temperature,  the  dew  point,  and  the  sea  water  temperature  (all 
in  “F.),  find  the  intersection  of  the  pertinent  air  temperature  and  dew 
point  lines  and  project  horizontally  to  the  left  to  the  given  sea  water 
temperature.  Read  the  resulting  AN* -value  from  the  sloping  lines  at 
the  point.  Then,  refer  to  figures  6,3,  6,4,  and  6,3,  to  obtain  the  mini¬ 
mum  AN'- value  required  for  ducting  as  a  function  of  radar  wavelength, 
surface  windspeed,  and  the  air-water  temperature  difference  (in  ®C.), 

Although  the  vertical  extent  of  an  oceanic  duct  may  be  estimated 
from  local  observations,  it  should  be  pointed  out  that  the  horizontal 
space  distribution  of  refractive  conditions  must  be  analyzed  from  a 
number  of  ship  observations  or  from  a  consideration  of  the  synoptic 
situation  (see  part  IV).  Unless  ducting  conditions  are  expected  to  pre¬ 
vail  over  a  large  horizontal  area,  extended  ranges  should  not  be  pre¬ 
dicted.  In  those  instances  when  ducting  is  indicated,  the  amount  of  range 
extension  cannot  be  accurately  computed.  However,  one  may  safely 
estimate  that  radar  coverage  will  be  extended  to  at  least  double  the 
normal  range.  •\r*icular!y  favorable  conditions  (strong  duct,  low 

antenna,  short  wavelength)  the  range  may  be  extended  as  much  as  five 
times  normal. 

6.1,3  Advection  Ducts  (Surface) 

As  was  discussed  in  section  5.1.1,  ducts  may  occur  when  warm, 
ar^  a:;'  flows  out  ovf  r  '  cool  sea.  They  do  not  extend  more  than  100 
miles  or  so  from  land  unless  there  is  a  strong  advection  wind  carrying 
the  air  offshore.  Usually  these  ducts  a.  e  found  over  seas  lying  close 


-  90  - 


X-SAWD  (3-6  CM.  ) 


Fig-jre  6.3.  X-Band  Trapping  Criteria— Oceat.:::  Duct. 


MINIMUM  DIFFERENCE  IN  REFRACTIVITY  BETWEEN 
SEA  SURFACE  AND  OBSERVATION  LEVEL— AN'mIn. 
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S-BAND  (6-20  CM.  ) 


Figure  6,4.  S-Sand  Trapping  Criieria--Oceanic  Duct. 


MINIMUM  DIFFERENCE  IN  REFRACTIVITY  BETWEEN 
SEA  SURFACE  AND  OBSERVATION  LEVEL -AN 'min. 
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L-oAND  (20  TO  iOO  CM.  ) 


Figiiru  6.5.  L-Band  Trapping  Criteria--Occanic  Duct. 
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to  desert  areas,  such  as  in  the  Mediterranean,  Red  Sea,  Guli'  of  Arabia, 
and  similar  areas.  To  assess  their  effects  on  propagation,  one  needs 
shipboard  radiosonde  or  other  data  from  which  to  plot  an  N-profile  to 
heights  of  500  feet  or  more.  In  order  to  determine  whether  the  duct  is 
capable  of  trapping  energy  of  a  particular  wavelengt'  'and  hence  whether 
extended  radar  coverage  maybe  expected),  theN-proi*ie  is  plotted  and 
the  height  of  the  layer  top,  hip,  and  ANt  for  the  layer  are  determined 
as  shown  in  figure  6.6.  Using  the  criter^w  urves  in  figure  6.7,  one  can 
determine  whether  a  particular  wavelength  will  be  trapped.  If  a  point 
corresponding  to  the  mea  sured  hp  and  AN  lies  to  the  right  of  the  curve 
corresponding  to  the  radar  frequency  in  question,  the  duct  will  trap  the 
emitted  energy  and  give  extended  range  coverage.  The  farther  to  the 
right  of  the  curve  the  point  lies,  the  greater  v/ili  be  the  radar  coverage 
extension. 

The  difference  between  the  trapping  criteria  presented  in  figure 
6.7  and  those  given  in  figures  6.3,  6.4,  and  6.5  lies  in  the  fact  that  fig¬ 
ure  6.7  is  applicable  to  any  duct,  but  one  must  have  an  N-profile  to  de- 


Figure  6,6.  N-Profile--Surface  Layer, 
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termine  whether  or  not  £  given  wavelength  will  be  trapped.  In  figures 
6.3,  6.4,  and  6,5,  on  the  other  hand,  one  needs  only  N-values  at  the  sea 
surface  and  at  bridge  height.  From  these  two  values,  one  can  deduce 
the  N-profile  if  it  is  assumed  that  the  air  has  been  over  the  sea  a  long 
enough  time  for  equilibrium  to  have  been  establisi.  In  general,  the 
criteria  shown  in  figures  6.3,  6.4,  and  6.5  are  applicaoie  far  from  land, 
while  those  in  figure  6.7  are  applicable  anywhere.  No  comprehensive 
comparison  between  the  two  methods  is  j-. .  available, 

6.1.4  Elevated  Ducts 

One  may  use  a  similar  procedure  for  analyzing  the  effect  of  ele¬ 
vated  ducts  on  radar  performance.  In  this  case,  plot  an  N-profile  as 
before;  figure  6.8  illustrates  an  elevated  duct.  Figure  6.7  may  again 
be  used  as  the  trapping  criterion  if  one  plots  the  point  corresponding  to 
the  height  of  the  layer  top  divided  by  two,  (-^)  and  ANt.  for  the  layer. 
If  the  point  lies  to  the  right  of  the  radar  frequency  curve  in  question, 
trapping  will  occur  and  radar  coverage  will  be  extended.  As  in  the  case 


Figure  6,8.  N-Profile--Elevated  Layer. 
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of  evaporation  and  adveotion  ducts,  one  cannot  accurately  predict  the 
extent  of  anomalous  covera^fe,  but  only  whether  or  not  it  will  occur. 
With  experience  one  may  estimate  whether  there  will  be  a  large  exten¬ 
sion  (up  to  5  times  normal  range)  or  a  moderate  one  (twice  normal) 
according  to  how  far  the  plotted  point  lies  to  the  right  of  i.  pertinent 
radar  fre'^uency  curve. 

The  criteria  curves  shown  in  figure  6,7  are  .u  ^sed  on  a  theoretical 
treatment  of  the  guiding  of  radio- radar  energy  by  ducts.  Nurcerous 
attempts  have  been  made  to  treAt  the  problem  theoretically,  but  this 
particular  method,  by  S.  A.  Schelkunoff  of  Bell  Telephone  Laboratories 
[18],  has  been  the  only  one  to  successfully  stand  the  test  of  experimental 
verification. 

6.2  Surface-to-Air  and  Air-to-Surface  Situations 


In  the  previous  discussion  it  was  noted  that  various  types  of  N- 
profiles  with  strong  gradients  could  channel  radar  energy  near  the  sur¬ 
face  and  thus  illuminate  targets  which  were  well  below  the  horizon  line. 
The  propagation  effects  of  refractive  layers  when  one  terminal  is  well 
above  the  ground  are  considerably  different  frOiH  those  when  both  ter¬ 
minals  are  low.  In  principle,  the  coverage  pattern  of  a  surface  radar 
looking  at  a  high-flying  target,  should  be  the  same  as  that  of  an  airborne 
radar  looking  at  surface  targets.  One  significant  difference,  however, 
is  that  in  the  air-to-surface  situation  "sea  clutter"  complicates  detec¬ 
tion,  particularly  of  small  surface  targets. 

6.2.1  Sea  Clutter 

The  myriad  of  echoes  frequently  observed  by  airborne  radars 
looking  over  a  sea  surface  is  called  "sea  clatter"  or  "scope  clutter". 
Its  effect  is  to  incrscis  «?  thp>  background  of  echoes  against  which  one 
must  detect  the  surface  target,- and,  in  the  case  of  tim-al3  targ'^ts  such 
as  submarine  snorkels  and  small  craft,  it  can  seriously  interfere  with 
detection.  Sea  clutter  is  caused  by  reflection  from  individual  waves. 
As  one  would  expect,  it  is  worse  in  heavy  seas  than  in  calm  ones,  and 
the  echoes  are  stronger  for  shorter  wavelengths  than  for  longer  ones. 
Furthermore,  the  echo  intensity  increases  as  the  cube  of  the  angle  of 
Thus,  clutter  echoes  extend  from  the  sea  surface  directly 
undert)  ;  rirc'-af-,  out  to  some  distance,  as  given  in  figure  6.9,  depending 
upon  the  altitude  of  the  aircraft  and  the  sea  st  de.  The  presence  of 
strong  N-gradients  has  little  effect  on  the  extent  of  the  clutter  echoes. 


Figure  6,9.  Graph  for  Kstimating  the  Radius  and  Intensity  of  Sea  Clutter. 
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Figure  6.9  was  adapted  from  a  chart  originally  produced  by  Lock¬ 
heed  Aircraft  Company  [23],  It  was  modified  to  its  present  form,  on 
the  basis  of  in-flight  sea-ciutter  data,  by  personnel  of  the  U.  S.  Fleet 
Weather  Central,  Argentia  [24J.  To  use  the  figure: 

1.  Enter  at  upper  right  with  swell  height  plus  wave  height 
in  feet,  or  with  the  coded  value  of  mean  maximum 
height  of  waves. 

2.  Proceed  hori7,ontally  to  the  appropriate  surface  wind 
speed  line  in  knots. 

3.  Proceed  vertically  down  to  the  appropriate  aircraft  alti¬ 
tude  curve, 

4.  Read  the  sea  clutter  radius  on  the  scale  at  lower  right 
in  nautical  miles. 


For  example,  suppose  the  wave  height  plus  swell  height  is  8  feet  and 
the  surface  wind  is  10  knots: 

If  the  aircraft  is  at  5,000  feet,  R  (the  radius  of  clutter)  =  64  miles. 
If  the  aircraft  is  at  10,000  feet,  R  (the  radius  of  clutter)  =  93  miles. 
Jf  the  aircraft  is  at  20,000  feet,  R  (the  radius  of  clutter)  =134 
miles. 


Because  of  the  effect  of  wave  shape  on  the  reflected  energy,  one  should 
add  15  percent  to  the  clutter  radius  in  the  upwind  area  and  subtract  20 
percent  in  the  downwind  area. 


As  stated  on  the  chart,  the  values  given  are  for  a  2  megawatt, 
lO-centimetex-  radar  ha"’ncf  an  antenna  gain  of  35  decibels.  To  obtain 
clutter  ranges  for  other  wavelengths  multiply  the  indicated  ranges  by 


10 


3/2 


X  (cm. ) 


for  other  radar  powers  multiply  the  indicated  ranges  by 


P  (mw. ) 


n  1/3 
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and  for  other  antenna  gains  multiply  the  indicated  ranges  by  the  antilog 
of  {o.087  [G  (db) -35j|  . 

6.2.2  Elevation  Angle  Error 

An  important  problem  encountered  in  surface,  aii  search  radar 
operation  is  that  of  elevation  angle  error.  The  standard  N-profile  causes 
radar  rays  to  bend  downward  as  they  trave..  vhrough  the  troposphere. 
The  total  ray  bending  along  the  tangent  line  amounts  to  about  3  /4‘’  as 
the  ray  proceeds  from  the  surface  up  through  the  stratosphere.  About 
90  percent  of  this  bending  takes  place  between  the  surface  and  30,000 
feet;  as  a  result  an  air- search  radar  capable  of  measuring  precise  ele¬ 
vation  angles  will  always  indicate  higher  elevation  angles  than  the  true 
ones. 


The  elevation  angle  error,  is  a  function  of  the  measured  eleva¬ 
tion  angle,  a^,  the  N-profile  and  the  target  height,  1^.  Standard  atmos¬ 
phere  corrections  can  be  made  to  take  average  bending  into  account 
(see  fig.  5.6).  In  the  presence  of  layers  where  the  N-gradient  is  stronger 
than  standard,  however,  the  bending  error  can  be  increased  consider¬ 
ably,  If  one  has  a  refractivity  profile,  the  radar  elevation  angle  error 
maybe  determined  quite  accurately,  according  toWeisbrod  and  Ander¬ 
son  [26],  in  the  following  way: 

1.  Subtract  all  N-values  (computed  from  radiosonde  data) 
from  the  refractivity  at  the  surface  (radar  height),  Ng. 
Tabulate  the  resulting  Ng  -  Nj^  values  for  each  height, 
h,  above  the  radar. 

2.  Using  figure  5.5,  determine  500  tan  ^  at  each  height. 

This  is  done  by  entering  the  left  margin  of  the  figure 
witxi  the  pT.«npr  Ng  -  N^.  Next,  move  horizontally  to  the 
curve  correspondiiig  to  the  height,  h.  which  is  expressed 
in  thousands  of  feet.  Now,  proceed  vertically  to  the  curve 
corresponding  to  the  measured  elevation  angle,  Qr^,  which 
is  expressed  in  milliradians.  Finally,  read  500  tan  ^  by 
moving  horizontally  to  the  right-hand  scale.  This  is  a 
graphical  solution  of  Snell’s  Law  (see  section 3.2).  Oc- 
casiorapiv.  -  Nj^  will  be  negative;  in  this  case  enter 
the  left  margin  with  Ng  -  Njj  and  use  the  dashed  height 
cur  /es,  proceeding  as  above  to  determine  500  tan  ,8 . 
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3.  By  dividing  the  atmosphere  into  layers  whose  boundaries 
are  the  radiosonde  heights  tabulated  above,  the  raybend- 
ing  At  caused  by  each  layer  can  now  be  determined  from 
the  expression 

ANt 

At  (milliradians)  = - = -  , 

500  tan  /5j-  +  500  tan  P>j, 

where  500  tan  is  the  value  tabulated  at  the  lower 
boundary  of  a  given  layer  and  500  tan  prp  is  the  value 
tabulated  for  the  next  higher  level;  ANt.  is  equal  to  N  at 
the  bottom  of  the  layer  minus  N  at  the  top  of  the  layer. 

For  each  consecutive  pair  of  heights  one  gets  a  value  of 
At  in  milliradians. 

4.  Add  the  At  s  from  the  radar  height  (surface)  to  the  tar¬ 
get  height,  h^.  The  sum  is  t^  in  milliradians. 

5.  Compute  6  corresponding  to  the  assumed  and  target 
height,  1^,  by 

1000  T^  tan  j3|.  -  (Ng  -  N; )  x  10  +-y 

6  (milliradians)  - - ,  (6.2) 

+  1000  tan  j3^  -  1000  tan 


where  is  the  refractivity  at  target  height  and  1,000 
tan  P  t  is  twice  the  value  of  500  tan  p  at  target  height. 

6.  This  will  give  a  value  of  £  corresponding  to  the  assumed 
Uq  and  target  height.  It  must  be  repeated  from  step  2  oq 
for  other  a q’s  and  target  heights  to  cover  a  useful  range 
of  conditions 

The  procedure  outlined  above  is  straightforward  though  tedious. 
It  is,  however,  quite  exact,  and  is  the  simplest  way  (short  of  using  an 
electronic  computer)  for  determining  6  from  radiosonde  data. 


One  can,  however,  estimate  the  values  of  6  to  be  expected  under 

reasonably  clos-  to  the  standard  stmosphere  by  using  figure 

5. Ill,  Ihs  values  shovn  are  for  an  Ng  of  350.  To  convert  to  other 

values  of  surfaci  refractivity,  multiply  the  indicated  elevation  angle 

errors  by.J!^fi..  The  presence  of  refractive  layers  can  cause  significant 
350 


ELEVATION  ANGLE  ERROR  IN  STANDARD  ATMOSPHERE  (8^,^  IN  DEGREES 
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departures  from  these  values,  particularly  at  low  elevation  ang.’es.  In 
figure  6.11,  the  differences  between  the  elevation  angle  error  under 
layered  and  standard-atmosphere  conditions,  6  layered  “linus  Sgjandard » 
are  shown  for  several  types  of  refractive  layers.  To  estimate  refrac¬ 
tion  effects  then,  one  adds  the  values  obtained  from  figure  6.11  to  those 
found  in  figure  6,10  to  obtain  -values  of  the  total  elevatio.  ngle  error. 

6.2.3  Range  Srror 

Because  radar  waves  in  the  atmosphere  travel  slower  than  they 
do  in  a  vacuum,  a  range  error  results.  Since  radar  reuige  calibrations 
are  made  using  the  vacuum  velocity  as  a  reference,  the  measured  ranges 
are  always  slightly  larger  than  the  actual.  The  resulting  error,  AR,  is 
usually  too  negligible  to  be  considered  in  practical  air  search  opera¬ 
tions.  It  amounts  to  approximately  700  feet  for  the  round  trip  to  a  tar¬ 
get  above  the  stratosphere  with  =  0,  Figure  6,12  shows  the  one-way 
range  error  as  a  fimction  of  measured  range,  R,  for  several  a^^'s  and 
assuming  ^  is  350.  For  other  one  may  assume  that  AR  is  pro¬ 
portional  to  Ng  with  sufficient  accuracy  for  most  practical  purposes. 
The  effect  of  strong  N-gradients  on  range  error  is  at  most  a  few  per¬ 
cent  of  the  values  shown  in  figure  6,12  so  it  can  safely  be  ignored. 

6.2.4  Radar  Holes 

The  occurrence  of  radar  holes  or  blind  spots  in  a  surface,  air- 
search  radar  coverage  pattern  is  infrequent.  It  was  noted  in  section 
5,1.2,  that  intense  elevated  layers  cause  such  holes  but  that  the  fields 
in  them  are  stronger  than  when  the  layer  is  absent,  so  they  do  not  re¬ 
duce  radar  coverage.  Weak  elevated  layers  are  incapable  of  causing 
surface  radar  holes  unless  they  are  just  a  few  hundred  feet  above  the 
radar.  Surface  ducts  caused  by  advection,  however,  can  cause  surface, 
air- search  radar  holes.  The  ray  tracing  results  for  this  case  also  apply 
to  low,  weak  eievateu  lav:”':.  In  both  cases,  the  radar  is  below  the  top 
of  the  layer.  The  neai-  edge  of  the  radar  hole  will  be  located  at  the  top 
of  the  layer  at  the  distance 

d  =  hrp\/ - - -  ,  (6.3) 

T  y  ANl  -  G.  048  h,j.  ' 

wnert:  I.,,  :::  the  height  of  top  of  the  layer  above  the  surface  and 
ANj^  .s  equal  to  N  at  the  surface  minus  N  't  the  top  of  the  layer,  Ng  - 
Nrp.  It  can  be  seer,  that  unless  the  N-gradient  t  tween  the  radar  amenna 
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Fig.’ure  6.12.  Range  Error  in  a  Standard  Atmosphere;  Ng  =  350. 
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height  (surface)  and  the  top  of  the  layer  is  greater  than  4  S  units  per 
1,000  feet,  d  will  be  imaginary  and  no  radar  hole  will  be  present.  The 
same  expression  applies  to  evaporation  ducts,  although  in  this  case 
as  defined  above  will  almost  never  be  large  enough  to  give  a  radar 

hole, 

6,2.5  Summary 

The  general  effects  of  refractive  layers  on  surface-to-air  (and 
air-to- surface)  radar  performance  are  summarized  following: 

1,  In  air-to- surface  situations  the  presence  of  "sea  clut¬ 
ter",  which  is  produced  by  reflections  from  individual 
waves,  may  seriously  interfere  with  the  detection  of  small 
surface  targets, 

2,  Actual  ray  bending  may  be  considerably  in  excess  of 
standard-atmosphere  bending.  This  will  cause  radar 
elevation  angle  errors  on  the  order  of  a  large  fraction 
of  a  degree, 

3,  Range  errors,  which  can  be  as  great  as  700  feet  over  a 
range  of  200  miles  in  a  standard  atmosphere,  are  not 
significantly  changed  by  the  occurrence  of  strong  N-gra- 
dients. 

4,  Radar  holes  or  blind  spots  in  the  radar  coverage  pat¬ 
tern  are  rarely  observed  in  surface-to-air  situations, 
but  may  occur  occasionally  in  the  case  of  advection  ducts, 

6.3  Air-to-Air  Situations 


Since  botr.  the  .aJ  jr  afid  +he  target  are  usually  well  above  the  sea 
surface,  evaporation  and  advection  layers  have  little,  if  anv  effect,  on 
radar  performance  in  this  instance.  Such  layers  are  too  far  below  .-ie 
line  of  sight  to  affect  propagrtion,  and  even  when  the  line  of  sight  grazes 
the  horizon,  they  have  little  effect  for  the  reasons  discussed  in  the 
preceeding  paragraph. 

6.3,1  P'.dav  Hedcs 


Elevated  layers,  on  the  other  hand,  :  an  be  quite  effective  in  causing 
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radar  holes.  It  was  shown  in  chanter  5,  that  radar  rays  can  be  totally 
refracted  when  the  radar  is  below  the  layer,  thereby  causing  a  region 
into  which  little  energy  can  penetrate. 

(a)  Location  of  Holes  (Radar  Below  Layer) 

The  distances,  R-  and  R^,  from  the  radar  to  the  near  edge  and 
to  the  far  edge  of  the  radar  hole,  respectively  '»t  the  top  of  the  layer) 
can  be  expressed  in  terms  of  the  refractivities^j^,  Ng,  and  N,j,  at  the 
radar,  layer  base,  and  layer  top,  respectively,  together  with  the  thi’ee 
corresponding  heights  hg,  hg,  and  h,j,  (see  fig.  6.13).  For  this  purpose 
two  parameters  Dj^  and  D2  are  defined 

0.232  (h,j,  -  hg) 

""  Ng  -  N.p  -  0.048  (h,p  -  hg)  ’ 


and 


0.232  (hg  -  hg) 

^  Ng  -  Ng  -  0.048  (hg  -  hg)  ' 

where  Dj  and  D2  are  in  nautical  miles  and  h^,  1^,  and  hg  are  in  feet. 
It  can  be  shown  that 

B  (naut.  miles)  =  (Dj^  -  Dg)  '-^Ng  -  N^  -  0.048  (h,j,  -  hg) 

±  D2  -  N^  -  0.048  (bp  -  hg)  .  (6.6) 

Subtracting  the  two  terms  gives  the  range  to  the  near  edge  of  the 
hole,  and  adding  them  gives  R”,  the  range  to  the  far  edge  (see  fig.  6.13). 
It  should  be  noted  tha:,  L2.  »*x'-cb  i:  d''terminc'^  by  the  N-gradient  below 
the  layer,  will  usually  be  negative,  since  tb.e  grad^'^ns  will  oe  less  than 
48  units  per  1,000  feet;  on  the  other  hand  will  always  be  positive 
when  there  is  a  radar  hole.  When  ihe  radar  is  at  the  height  of  the  layer 
base,  the  two  radicals  are  equal.  When  the  radar  is  within  the  layer, 
Dj  and  are  equal, and  the  range  to  the  near  edge  of  the  radar  hole  is 
given  by  the  last  term  of  the  range  equation. 

The  width  of  the  radar  hole,  W,  at  the  level  of  the  top  of  the  layer 
is  approximately  given  as  follows: 
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W  (naut.  miles)  2  Dg  -  0.048  (h-j.  .  (6.7) 

The  foregoing  equation  assumes  the  radar  hole  width  is  equal  to  the 
range  to  the  far  edge  minus  the  range  to  the  near  edge;  i.e,,  V/  =  R’*  -  R*. 
In  figure  6,13,  we  see  that  this  is  only  strictly  true  if  theap'»ie  between 
R*  end  K”  equals  zero;  in  all  practical  cases  this  angle  wiix  be  small 
and  the  error  introduced  by  the  assumption  will  be  entirely  insignificant, 

(b)  Avoidance  of  Holes  (Radar  Below  Layer) 

It  was  pointed  out  in  chapter  5,  that  if  the  radar  is  sufficiently 
far  below  the  layer,  there  will  be  no  radar  hole  at  all.  TMs  can  be  seen 
from  equation  6.5;  if  [o.048  (h,j,  -  hgy[  becomes  equal  to  •* 
the  radical  will  be  equal  to  zero,  and  the  radio  hole  will  have  zero  width. 
For  a  still  lower  hp^,  the  hole  width  will  become  imaginary  and  no  hole 
will  exist.  The  radar  height  at  which  a  radar  hole  will  first  appear  is 
therefore. 


hj^  (max.) 


0.04811^  -  Ng  +  Nrj, 
0.048hjj  -  Ng  +  Ng 


(6.8) 


where  ^  is  the  surface  refractivity, 

(c)  Location  of  Holes  (Radar  Above  Layer) 

When  the  radar  is  above  the  elevated  layer,  the  distances  to  the 
near  and  to  the  far  edges  of  the  radar  hole  are  given  at  radar  height 
instead  of  at  the  bright  of  the  layer  top  as  was  done  previously  (see  fig. 
6.14),  The  near  edge  distance,  is  given  by 

0  '165  (h^  -  h^) 

R’  (naut,  miles)  =  ■  ■  - - — -  - ^  (8.3) 

^yo.048  (hj^  -  h,^')  -  (N,j.  -  Ng) 


where  all  heights  are  expressed  in  feet.  The  width  of  the  radar  hole, 
when  the  radar  is  above  the  elevated  layer,  turns  out  to  be  a  complex 
express, which  is  plotted  g^s.'hically  in  figure  6.15.  The  abscissa 
is  C2/C^  v/l.'Ci  s;  aud  C2  are  defined  as 
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Fij^re  6.15,  Width  of  Radar  Hole  Chart--Radar  Above  Layer. 
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-  Ng 

^2  - +  0.048  ,  (6.10) 

h,p  -  hg 

tuid 

Ng  - 

Cj  =  - +  0.04b  .  (6.11) 

Hg  - 

The  ordinate  is  Ahj^/Ah2,  where  is  the  height  of  the  radar  above 
the  layer  top,  (hg  -  h^),  and  Ahy,  is  the  layer  thickness,  (hfji  -  hg). 
The  family  of  curves  give  hole  width,  W,  at  radar  height,  divided  by  j^, 
the  distance  to  the  near  edge  of  the  hole  at  radar  height. 

In  the  expressions  for  C  j  and  above,  it  will  be  seen  that  Cj  is 
always  positive,  but  that  for  strong  layers  can  become  negative. 
Thus,  in  figure  6,15,  negative  values  of  C2/C1  refer  to  strong  layers 
and  positive  values  to  weak  layers.  When  C2/C1  is  zero,  the  N-gradient 
in  the  layer  is  equal  to  the  earth’s  curvature,  and  when  C2/Cj  is  equal 
to  positive  one,  the  layer  becomes  indistinguishable  since  its  N-gradient 
is  equal  to  that  above  the  layer. 

The  curves  of  figure  6,15  were  adapted  from  theoretical  work 
carried  out  at  Cornell  University  [11],  and  the  results  have  been  com¬ 
pared  with  actual  flight  data  taken  by  the  Wright  Air  Development  Cen¬ 
ter.  It  was  found  that  actual  elevated  layers  produce  radar  holes  whose 
location  and  width  agree  quite  well  with  the  theory.  Furthermore,  it 
was  found  that  appreciable  radar  holes  occur  even  when  C2lCj^  is  as 
much  as  0.9  (as  indicated  by  radiosonde  data).  Since  this  value  of  C2/ 
Cl  is  associated  7/i+*’an  N-gradient  in  the  layer  only  30  percent  greater 
than  standard,  il  is  evident  tnat  even  very  weak  elevated  layers  can 
cause  significant  radar  holes. 

(d)  Numerical  Examples 

It  is  likely  that  the  descriptions  given  in  paragraphs  6.3.1(a),  ('<), 
['')  above  for  locatuv?  radar  holes  involve  concepts  unfamiliar  to  the 
•■eadei.  Vo  clarify  these  concepts,  numerical  examples  are  given  in  the 
following  subparagraphs.  Suppose  the  radiosonde  data  given  in  table  6.1 
is  representative  over  the  area  where  -.  adar  holes  are  to  be  located. 
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TABLE  6.1 


Pressure 

Temperature 

Dew  Point 

Height 

N 

mb. 

“C. 

“C. 

ft. 

N- Units 

1,013 

15.0 

12.0 

336 

902 

8.4 

7.0 

3,260 

296 

885 

12.2 

-13.0 

3,780 

850 

9.5 

-18.0 

4,900 

240 

700 

-1,3 

-35.0 

9,900 

201 

Sample  Radiosonde  Data,  Heights,  and  N-Values  for  the  Numerical  Examples. 


The  first  step  is  to  compute  N  at  each  level,  using  chart  I  or  chart  II; 
these  values  are  tabulated  in  the  last  column  above.  Next  plot  N  versus 
height  to  obtain  a  refractivity  profile  as  described  in  section  3.3.2  (see 
fig.  6.16).^  From  this  profile  the  parameters  of  the  trapping  layer, 
which  is  present,  are  determined  to  be 


Ng  =  296 

hg  =  3,260  feet 

N,j.  =  250 

h,j.  =  3,760  feet 

AN  =  Ng  -  N.J.  =  46 

Ah  =  hip  -  hg  =  500  feet 

(i.)  Radar  Below  the  Layer 

To  determine  the  radar  hole  location  when  the  radar  is  below  the 
layer,  substitute  the  above  values  into  equations  6.4,  6.5,  and  6.6  as 
follows: 


0.232  (500) 

Di  =  - 

46  -0,048  (500) 


116 

- -  5.28  nautical  miles  , 

46  -  24 


^"'Jienv  ♦nod  lie'  .ribed  ,.3.4  will  give  somewhat  different  numerical  answers, 

ro  •  tile  purposes  of  illustration  the  simpler  presentation  (of  N-profiie)  is  used. 
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Figure  6.16.  N- Profile- -Numerical  Example. 


0.232  (3,260  -  hj^) 

= - 

Nj^  -  296  -0.048  (3,260  -  hj^) 


and 

R  =  (5.28  -  D2)\/46  -  24  ±  DgV  Nj^  -  250  -0.048  (3,760  -  hj^)  . 

To  determine  the  actual  location  of  ths;  rada”  hole,  assume  a  radar 
height  hj^  and  read  off  1^  at  that  height  from  figure  6.16. 

Example  1 


Suppose  =  3,000  feet;  Np.  (taken  from  fig.  6.1G)  is  then 


301.  D2  now  becomes 


Q.»  23X^2^Q.1=  -8.05  nautical  miles, 
5  -  12. 0 
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R  =  13.3  \J22  ±  8.0  /51-  36.5  , 

R  =  62  ±  3fi  nautical  miles  . 

The  near  edge  of  the  radar  hole,  at  the  layer  top  (3,760 
ft.)  is  thus  at  62  -  30  =  32  nautical  miles  and  the  far  edge, 
Rj^,  is  at  62  +  30  =  9?  nautical  miles. 

Example  2 


Now  lower  the  radar  height  to  2,500_feet;  (from  fig.  6.16) 
then  Np  =  307.  Do  becomes  =  "2.9  nautical 

miles, 

R  =  12.2  \/22  ±  6.9  pi  -  60.5  . 

Since  the  quantity  under  the  second  radical  is  negative,  the 
term  is  imaginary,  and  there  will  be  no  radar  hole. 


From  the  above  examples,  it  is  seen  that  a  radar  hole  is  present 
when  the  radar  height  is  not  far  below  the  layer  base,  but  is  absent  when 
the  radar  is  lowered  an  additional  500  feet.  For  this  particular  layer, 
one  may  find  the  maximum  radar  height  which  will  not  produce  a  radar 
hole.  Substituting  the  parameters  of  the  example  into  equation  6.8  we 
have 


hj^  (max.) 


3,260  (180  -  86) 
156  -  40 


2,640  feet  . 


In  this  example,  a  radar  height  of  2,640  feet  will  give  a  radar  hole  width 
of  zero;  at  greater  heights  a  radar  hole  v;ill  appear.  This  then  is  the 
critical  heigiu  for  avoiding  a  radar  hole  for  the  particular  profile  of  the 
example.  Using  equation  6.8  one  can  readily  compute  h„  (max.)  for  any 
elevated  layer,  using  appropriate  values  of  hg,  h,j„  Ng,  Ng,  and  N.^. 

(ii.)  Radar  Above  the  Layer 

Consider  now  the  problem  of  locating  the  radar  hole  when  the 
radar  is  above  the  same  layer  as  specified  in  figure  6.16. 

Example  3 

For  comparison  with  example  1  above,  assume  (hg  -  hp) 
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in  the  present  case  is  equal  to(hg  -  hjj)  in  the  previous  case, 
or  hj^  =  4j020  feet:  thenfrom  figure  Sje,  Nj^  =  248  and  sub¬ 
stituting  into  equations  6.9,  8.1?,  and  6.11  we  have 

0.465  (260)  121 

R'  =  “puzznZZZZl  = -  =38  nautical  miles  , 

\/0.048  (260)-2  3.2 

C2  =  -^  +  0.048  =  -J.044  , 


and 

Cl  =-^+  0.048  =  0.040  . 


Then  C2/C1  =  -1.1,  and 


Entering  these  values  in  figure  6.15  gives  =  1.8.  Thus, 
the  hole  width  at  radar  height  is  given  as 

W  =  1.8  (38)  =  68  nautical  miles. 

The  distance  to  the  near  edge  of  the  radar  hole  at  radar  height 
is  38  nautical  miles,  and  to  the  far  edge  it  is  68  plus  38  or 
106  nautical  miles. 

E:;ample  4 


In  this  example,  assume  (hp^  -  h^,)  is  equal  to  (hg  -  hg) 
in  exaiiv.  Ic  760  feet.  Then  hg  =  4,520  feet  and 

from  figure  6.16,  Ng  •“  243.  In  this^'^se,  from  equation  6,9, 
we  have 

0.465  (760)  354 

R*  =  —  "  "■= - =  66  nautical  miles  , 

\/0.048  (760)-7  5.4 

and  anu  Cp  will  be  the  same  as  in  example  3  above  (C2/ 
Cp  =  1.1).  The  ratio  Ahi/Ah2  will  now  be --  1.5  and 
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W 

from  figure  6.15,  ^1=  0.94,  Thus, 

W  =  0,94  (66)  =  62  nautical  miles  , 


and 


R*»  =  66  +  62  =  128  nautical  miles  . 

The  results  obtained  in  examples  3  and  4  show  that,  when  the 
radar  is  not  far  above  the  layer  top  (260  ft.),  the  near  edge  of  the  radar 
hole  is  38  nautical  miles  away  and  has  a  width  of  68  nautical  miles. 
Raising  the  radar  height  500  feet,  moves  the  near  edge  of  the  hole  out 
to  66  nautical  miles,  but  does  not  remove  the  hole.  When,  however,  the 
radar  is  an  equal  height  below  the  layer  base  (example  2),  the  hole  dis¬ 
appears  entirely.  This  illustrates  the  statement  made  previously; 
namely,  if  the  rad*^r  is  above  the  layer,  raising  its  height  moves  the 
hole  outward  and  narrows  it  slightly,  but  does  not  eliminate  it.  When 
the  radar  is  below  the  layer,  one  can  eliminate  the  hole  entirely  by 
locating  the  radar  several  hundred  feet  below  the  layer  base. 

(e)  Shape  of  Radar  Holes 

The  previous  paragraphs  have  discussed  the  location  of  the  radar 
hole  (1)  where  it  intersects  the  top  of  the  elevated  layer,  for  the  case 
where  the  radar  is  below  the  layer,  and  (2)  where  it  intersects  the  radar 
altitude,  for  the  case  wh{  .  the  radar  is  above  the  layer.  Since  the 
radar  target  will  rarely  be  at  either  of  these  particular  altitudes,  one 
would  like  to  be  able  to  locate  the  boimdaries  of  the  hole  at  any  target 
height. 

(i.)  Radar  Below  the  Layer 

When  the  raJar  is  tclow  ihc  layer,  is  particularly  simple. 
Since  the  rays  defining  both  edges  of  the  emerge  tangent  to  the 
layer  top,  they  will  remain  the  same  distance  apart  at  all  heights  above 
the  layer  top;  in  other  words,  the  width  of  the  radar  hole  above  the  layer 
is  constant  with  altitude.  This  effect  is  illustrated  in  figure  6.13.  At 
a  height,  1^,  (above  the  top  of  the  layer)  the  radar  hole  width  is  equa 
to  W^;  comparison  will  show  that  W  -  W^.  The  only  problem  then  is  to 
o.'fin''  the  Tange  be''unU  the  layer  top,  versus  altitude  (above  the 
layer  top)  relatic.iship  of  one  edge  of  the  hole  (see  fig.  6,13),  and  the 
problem  is  compL.'tely  solved.  Figure  6.3  7  shows  various  curves  de- 
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R)(-  RAWOi.  BEYCNt>  LAYER  TOP  (NAUTICAL  MILES) 
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200  2S0 
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350 


"^^''^LAYEHTop 


OS  op 


1.  N-GRAOIENT  ABOVE  LAYER  =0. 

2.  N-GRAOIENT  EXPONENTIAL  tBOVE  LAYER 

Ng  ‘  250  (dry  ahnotphere). 

3.  N-GRAOIENT  EXPONENTIAL  ABOVE  LAYER 

N{  :400  (moist  atmosphere). 

RADAR  HEIGHT^  1,000 FT.  ABOVE  SURFACE. 

4.  N-GRADIENT  EXPONENTIAL  ABOVE  LAYER 

N,  :  400  (moist  atmosphere). 

RADAR  HEIGHT  =  5,000  FT.  ABOVE  SURFACE. 


Figure  6.17.  Sliape  of  Radar  Hole  Chart, 

fining  this  relationship  for  several  different  refractive  co.ntiitions.  In 
the  figure,  the  layer  top  is  shown,  together  with  segments  of  curves 
drawn  at  10,000-foot  increments  above  the  layer  top.  The  range  scale, 
is  in  nautical  miles  beyond  the  point  where  the  ray  emerges  from 
the  layer  top.  The  exact  trajectory  of  the  ray  will  depend  on  tlie  char¬ 
acteristics  of  the  N-profile  above  the  layer.  The  upper  straight  line, 
curve  1,  assumes  no  change  of  N  with  height  above  the  layer,  which  of 
course  is  unrealistic.  The  second  curve,  2,  assumes  a  reasonably  dry 
atmosphere  a’covt  the  iayei,  which  co>.  re.‘=ponds  to  an  exponential  pro¬ 
file  with  Ng  equal  to  250.  The  next  twocu»  assume  a  vary  high  sur¬ 
face  N  of  400.  In  these  cases  the  amount  of  atmosphere  below  the  source 
is  of  importance;  thus,  the  third  curve,  3,  is  for  a  radar  height  of  1,000 
feet,  and  the  fourth  curve,  4,  is  for  a  radar  at  5,000  feet.  An  exponential 
decrease  of  N  with  height  above  the  layer,  as  portrayed  in  figur*  3.5, 
lias  been  assumed  i~  all  cases  except  the  top  line.  The  proper  curve 
to  use  in  most  ^ases  will  lie  between  2  and  3,  and  probably  closer  to 
curve  2. 
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(ii.)  Radar  Above  the  Layer 

In  this  situation,  the  near  edge  of  the  raaar  hole  is  again  deiined 
by  figure  6,17,  since  this  edge  is  tangent  to  the  layer  top.  The  far  edge 
of  the  hole  is  defined  by  the  next  ray  emerging  from  the  layer  top,  which 
emerges  at  a  small  upv/ard  angle  rather  than  being  exactly  tangent  (see 
fig,  6,14),  The  hole,  therefore,  gets  narrower'  (width  u  '•reases)  with 
increasing  altitude  above  the  layer,  at  a  rate  .vhieh  depends  on  the  layer 
characteristics.  The  parameter  given  in  figui  e  6,15  is  the  ratio  of  hole 
width,  W,  to  the  distance  to  the  near  edge,  at  the  same  altitude  as 
the  radar  antenna.  It  turns  out  that  above  this  altitude  the  hole  narrows 
very  slowly  with  altitude;  thus,  if  one  predicts  a  hole  of  constant  width, 
he  will  be  a  little  pessimistic  in  that  the  hole  at  greater  heights  will 
actually  be  somewhat  narrower  than  predicted, 

(f)  Practical  Limitations 

In  applying  ray  tracing  techniques  to  pr  .ctical  propagation  prob¬ 
lems,  several  limitations  should  be  kept  in  mind.  First,  elevated  layers 
as  they  occur  in  nature  are  rarely  uniform  over  the  100-200  mile  ranges 
called  for  in  the  foregoing  ray  treatment.  Even  if  one  had  an  accurate 
refractive  index  profile  at  one  point  along  the  propagation  path,  this 
would  not  necessarily  be  representative  of  the  layer  characteristics 
100  miles,  or  even  50  miles,  away.  Usually  elevated  layers  can  be  re¬ 
garded  as  being  reasonably  uniform,  but  one  must  remember  that  a 
single  radiosonde  ascent  may  not  be  representative  of  the  layer  in  ques¬ 
tion,  It  is  always  necessary  to  consider  carefully  all  available  synoptic 
information. 

Another  characteristic  of  layers  is  that  they  have  small-scale 
irregularities  which  make  them  more  like  a  ground  glass  mirror  than 
a  smooth  mirror  in  the  region  of  total  refraction.  This  has  the  effect 
of  scattering  energy  in'" the  radar  hole  from  adjacent  rays.  This  "fill¬ 
ing  in"  of  the  theoretical  well-defined  radar  hole  is  ino*~c  apparent  at 
the  edges  of  the  hole,  particularly  at  shorter  wavelengths,  since  these 
wavelengths  are  more  affected  by  small- sized  roughness  elements. 
For  example,  the  southern  California  subsidence  inversion  layer  looks 
quite  smooth  at  a  3-meter  wavelength,  but  it  looks  quite  rough  to  0.5- 
meter  energy.  As  pointed  out  in  chapter  5,  weaker  layers  appear 
sri'outl'cr  than  intense  Iry-.-.  s;  for  this  reason  when  a  hole  is  present 
dui  TO  a  Weak  layer  u  is  less  likely  to  be  filled  in  by  scattering;  i.e.. 
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Vr'ill  have  more  well  defined  boundaries^ 


Another  complicaiing  factor,  which  is  particularly  common  over 
the  sea,  is  that  energy  is  reflected  from  the  sea  surface.  This  energy 
will  be  incident  upon  the  layer  at  much  steeper  angles  than  the  direct 
ray  and  thus  will  not  be  affected  by  the  layer  to  any  appreciable  extent, 
Energy  reflected  from \he  sea  will,  therefore,  tend  to  fill  in  radar  holes, 
particularly  at  longer  wavelengths  (due  to  more  e*  'oient  reflection). 
The  fact  that  deep  radar  holes  are  obs  'rved,  however,  even  over  sea 
paths,  attests  to  the  fact  that  one  should  not  rely  on  sea  reflections  to 
eliminate  radar  holes  hut  rather  only  to  reduce  their  intensity. 

The  above  remarks  are  not  intended  to  convey  the  impression  that 
atmospheric  structure  is  so  erratic  that  useful  predictions  of  its  effect 
on  radio  or  radar  propagation  cannot  be  made.  It  is  wise,  however,  to 
realize  that  various  complicating  factors  are  operating  to  modify  the 
simple  assumptions  on  which  the  computations  are  based.  Observa¬ 
tional  evidence  shows  that  radiosonde  and  refractometer  data  are  usually 
adequate  for  practical  performance  calculations, 

6,3,2  Extended  Radar  Ranges 

In  the  case  of  air-to-air  layered  situations,  it  is  found  that  the 
energy  in  the  region  just  beyond  the  far  edge  of  the  radar  hole  is  stronger 
than  when  no  layer  is  present.  Theory  [11]  predicts  that  this  should  be 
the  case,  and  experimental  evidence  supports  the  theory.  As  the  tar¬ 
get  closes  range  it  will  first  pass  through  this  region,  which  is  usually 
about  the  same  width  as  the  radar  hole.  As  it  moves  in,  the  target 
passes  through  numerous  narrow  fingers  of  alternately  high  and  low 
illumination.  The  field  in  the  maxima  increases  as  the  far  edge  of  the 
radar  itOle  is  approached,  and  usually  detection  probability  increases 
markedly  over  that  expected  if  no  layer  were  present.  Finally,  when 
the  far  edge  of  the  radar  hole  is  passed,  the  target  will  be  lost  until  it 
emerges  fron.  th.- .iC-'r  edg*'  of  theradfir  hole  at  relatively  close  range. 
The  layer,  therefore,  serves  to  extract  radar  eneigy  from  the  region 
of  the  radar  hole  and  to  concentrate  it  in  the  region  beyond.  As  .cn- 
tioned  above,  the  enhancer'.ent  is  strongest  just  beyond  the  radar  hole; 
it  decreases  gradually  through  a  distance  about  equal  to  the  hole  width, 
and  merges  with  the  free  space  field  beyond. 

'’Vhen  a  l?ye:  >6  present  which  causes  a  radar  hole,  the  best  oper¬ 
ational  procedure  is  to  fly  the  radar  above  the  layer  at  a  height  such 


-  120  - 


that  the  far  edge  of  the  radar  hole  coincides  with  the  radar’s  normal 
detecticn  range*  In  t£ii&  wajr*  the  enhancement  heyond  the  radar  hole 
will  serve  to  increase  the  detection  range  by  at  least  half  the  width  of 
the  radar  hole. 

6.4  Method  for  Dcterminii.,^  Refractive  Characterist'^^s  from  Radar 
Performance  Bata 


Under  certain  circumstances,  such  as  in  aircraft  carrier  opera¬ 
tions  during  radio  silence,  ont !  may  net  be  able  to  obtain  radiosonde  data 
fromwliich  to  make  predictions,  and  at  the  same  time  airborne  refrac- 
tometers  may  not  be  available.  It  may  be  feasible,  hov/ever,  to  obtain 
initial  detection  ranges  on  incoming  aircraft  at  several  altitudes  with 
surface  search  radar.  Such  information  will  serve  to  define  the  shape 
of  the  radar  horizon  line  and  thus  indicate  the  presence  of  an  elevated 
layer. 


It  was  stated  in  chapter  5,  that  the  range  to  the  radio  horizon  in 
an  atmosphere  with  a  linear  N~gradient  is 

R  =  /2ae  hj^  +  , 

where  ag  is  the  effective  earth’s  radius  for  the  N-gradient  present, 
is  the  radar  height,  and  ht  the  target  height.  If  one  plots  ^  versus  h^, 
for  various  N-gradients,  the  result  is  shown  in  figure  6.18.  The  dis¬ 
tance  is  made  up  of  two  parts,  corresponding  to  the  two  terms  expressed 
above.  The  first  t  srrn  is  the  distance  from  the  radar  to  the  radio  hori¬ 
zon  ( hj^)  and  is  plotted  to  the  left  of  zero.  Since  this  term  is 
usually  quite  small  compared  to  the  second  term,  the  effect  of  varying 
ag  is  insigniticant,  and  the  standard  N-gradient  of  12  tinits  per  1,000 
feet,  has  been  assiimed.  If  the  radar  is  at  a  height  of  100  feet,  the  hori¬ 
zon  distance  is  12  miles.  If  the  N-gradient  is  constant,  the  total  hori¬ 
zon  distance  for  a  ;;ivtVi  tu.  get  is  the;'.  1 2  miles  plus  the  horizon 

distance  at  which  the  cur\e  corresponding  to  the  sting  gradient  inter¬ 
sects  the  target  height,  1^. 

6.4.1  Layer  Characteristics 

Tf  one  has  initial  detection  range  data  at  several  target  heights 
ai;  I  plcLit.h;--m  on  fig  ire  6.18  (after  subtracting  the  radar  height  term), 
the  points  v/ill  lie  on  a  single  curve  if  no  elevated  layers  are  present. 
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Figuri  6.18  Height  of  Radar  Horizon  versus  Range  and  Refractivity  Gradient.  (Chart  VI) 


o  1  OO  - 


In  this  case  the  position  of  the  plotted  points,  as  shown  by  the  circles, 
in  relation  to  the  curves  already  plotted  on  figure  6,18,  enables  one  to 
determine  the  N-gradient  directly. 


When  an  elevated  layer  is  present,  the  points  obs^  -ed  for  the 
targets  (aircraft)  below  the  layer  height  will  Me  on  a  -’ngle  curve  as 
before.  When  the  targets  are  above  the  layer,  n.  wever,  ti,^  points  will 
lie  further  out,  as  shown  by  the  crosses  in  figure  6,18,  The  reason  for 
this  is  that,  as  the  rays  penetrt.te  the  layer  they  are  bent  dov-mv/ard  (as 
illustrated  in  fig,  5.2),  thus  extending  the  horizon  somewhat. 


Using  points  plotted  as  in  figure  6.18,  one  can  determine  the  layer 
characteristics  in  a  particular  situation  in  the  following  way: 


1,  Observe  that  the  lower  points  (crosses)  lie  on  the  curve 
corresponding  to  an  N-gradient  of  15  units  per  1,000  feet. 
Thus,  one  can  say  that  this  is  the  N-gradient  below  the 
layer. 

2.  Next,  place  the  overlay  (fig,  6,19)  on  top  of  figure  8.18. 
This  overlay  is  simply  a  duplicate  of  the  curves  plotted 
in  figure  6.18  to  the  right  of  the  origin  (h  =  0,  R  =  0). 
Now  move  the  overlay  up  and  to  the  right,  keeping  the 
axis  at  the  bottom  parallel  to  the  horizontal  lines  in  fig¬ 
ure  6.18,  A  position  of  the  overlay  will  be  found  such 
that  'vhe  upper  points  (crosses)  fall  on  one  of  the  overlay 
curves.  When  this  position  is  found,  the  overlay  curve 
corresponding  to  the  plotted  points  indicates  the  N-gra¬ 
dient  abovf^  'he-  (4?  ~  12  in  case). 

3 .  Read  off  the  height  h  j ,  and  range  Rj^ ,  where  the  two  curves 
intersect  (•^'  =  15  of  tigure  6,18  and  ^  =  12  of  the 
overlay).  Inthiscase  hj  =  3,500feet  and  Rj^  =  76  miles. 

4,  Head  on  xht  a2»  and  distance  R2,  of  the  overlay 

origin  (cross).  In  this  case  h2  =  1,100  feet  and  R2  =  16 
miles. 
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5. 


The  intersection  height  of  the  two  curves  h2,  is  th>  layer 
height.  The  total  change  of  N  through  the  layer  is  deter¬ 
mined  by 


where  0.0541  =  Nj  and  0..  j41 

In  thi.s  case 


Rji  -  R2^ 


(5.12) 


hj  3,500 

—  = - =  46.0  feet /mile  , 

Rj  76 


and 

hj  -  hg  3.500  -  1,100  2,400 

- -  -  = -  =  40.0  feet/mile  , 

Rj  -  Rg  76  -  16  60 


6,  Figure  6,20  is  a  plot  of  (or  Ng)  versus-^  .  Nj,  cor¬ 
responding  to  (-j^)  =  46,0,  is  115,  and  N^,  corresponding 
to  r|)  =  40.0,  is  87.  Thus,  AN^  =  Nj  -  Ng  = 

115  -  87  =  28  iinits. 


In  summary,  the  points  plotted  on  figure  6,18  indicate  that  theN-profile 
has  agradient  of  15  units  per  1,000  feetup  to  3,500  feet,  then  a  decrease 
of  28  units,  followed  by  an  N-gradient  of  12  luiits  per  1,000  feet  above 
this  height.  It  is  not  feasible  to  determine  the  layer  thickness  from 
detection  range  data.  One  would  need  a  large  number  of  very  precisely 
measured  points  and  would  have  to  assume  a  smooth,  absolutely  uniform 
layer  in  order  to  derive  the  layer  thickness.  Instead,  one  has  obtained 
the  effective  wmch  the  would  have  if  it  had  zero  thickness; 

i.e.,  if  it  were  a  very  sharp  (tliin)  layer. 

In  addition  to  figures  j,18  and  6,19,  full-size  reproductions  (with¬ 
out  the  plotted  data  show'n)  have  been  printed  on  transparent  paper  and 
are  provided  as  loose-leaf  inclusions  to  this  manual  (charts  VI  and  '•’’IJ). 
.'i+hnugh  charts  "VI  andVII  maybe  used  in  their  present  form,  it  is  sug¬ 
gested  that  working  copies  be  made  on  transparent  plastic  stock  using 
the  Ozalid  'or  otherl  duplicating  process.  Routine  use  of  these  charts 
in  their  present  form  will  soon  result  ..i)  their  destruction. 


Figure  6.20.  Refractivity  versus  h/R. 


6,4.2  Prediction  of  Radar  Holes 
(a)  Radeir  Below  the  Layer 


When  an  airborne  radar  is  feet  below  a  sharp  layer,  the  dis¬ 
tances,  ^  and  R”,  to  the  boundaries  of  the  I'adar  hole  are  given  by 


ti  (ixaut.  miles)  r 


1  -0.0209(^ 


- 


ANl  -  0.048h'[l  -0.0209(^)j]  I  .  (6.13) 


]  .  (6.3 


where  is  the  N-gradient  below  the  layer  in  units  per  1,00C  feet 

aiid  AN.  ij  N  U.e  Aayer  itself.  Using  the  minus  sign 

before  tue  second  radical  gives  the  distance,  to  the  near  edge  of  the 
hole  at  layer  height,  and  using  the  plus  sign  give  a  the  distance,  ]^,  to 
the  far  edge. 
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Using  the  figures  given  in  the  preceding  example  we  obtain 


R  = 


4.8? 


\/28  ±  \/28  -0.048[1  -0.0209(15)]hn  .  (6,13) 


or 


1  -0.0209(15) 

R  =  7.1  jj.S  ±  ^28  -0.033h« J  , 
R  =  37.5  Ijl  ±  ''0118h»]  . 


If  the  radar  heigh-*  is  850  feet  or  more  below  the  layer,  the  radical  will 
be  imaginary  and  no  hole  will  be  found.  If  it  is  closer  to  the  bottom  of 
the  layer  than  this,  say  500  feet  below  the  layer,  there  will  be  a  hole 
between  13.5  miles  and  61,5  miles. 


(b)  Radar  Above  the  Layer 


When  the  radar  is  above  a  sharp  layer,  the  distance  to  the  near 
edge  of  the  radar  hole  at  radar  height  is  similar  to  that  given  in  section 
6,3,  namely. 


R’  (naut.  miles)  * 


(6.14) 


where  ^  is  the  r^dar  height  in  feet  above  the  layer  and  (^■)2  is  the 
N-gradient  in  N-units  per  1,000  feet  above  the  layer.  The  distance  R” 
to  the  far  edge  of  the  hole  (again  at  radar  height)  is  obtained  from  fig¬ 
ure  6,21,  In  this  figure,  it  is  assumed  that 


~  units/1,000  feet  , 


which  will  be  ably  well  satisfied  most  of  the  time.  This  figure 

is  similar  to  figure  6,15,  which  applies  to  layers  of  fiidte  thickness,  and 
gives  the  ratio  of  hole  width  W  to  ^  in  terms  of  and  radar  heit'‘  I, 
h**,  above  the  layer. 

(c)  Summary 

To  cuniinarize  use  of  the  radar  method  for  obtaining  refrac¬ 
tive  characteristics,  one  finds  that  if  initial  detection  ranges  on  incom¬ 
ing  aircraft  are  available  for  at  least  f  ur  different  aircraft  altitudes. 


Figure  6.21.  Width  of  Ruaar  Hole  Chart--Radaf  Above  a  Sharu  Discontinuiij  in 
N-Gradient. 


it  is  possible  to  deduce  the  effective  characteristics  of  the  N-profile, 
These  characteristics  are:  (1)  theN-gradient,  (■^•>2,  belowthe  layer, 
(2)  the  layer  height.  (3)  the  effective  change  of  N  through  the  layer, 
AN^,  anu  \4;  Lhe  N-gi-adJerit,  Ig*  above  the  layer.  The  effective  N 
change  hi  the  layer  is  based  on  a  thin  layer,  and  it  is  not  practica)  fo 
attempt  to  determine  layer  thiclmess.  Using  the  de.  i  ved  characteristics. 
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one  can  predict  the  location  and  extent  of  radar  holes  for  airborne 
radars.  In  order  to  do  this  without  layer  thickness  information,  one 
must  modify  the  methods  of  section  8.3  as  has  been  shown. 

It  should  be  emphasised  that  theN-profile  characteristics  derived 
in  this  way  are  dependent  upon  the  accuracy  withwhic'  *he  horizon  line 
of  ship-borne  radar  can  be  defined.  As  'an  be  seen  xn  figure  8.1s, 
errors  of  a  mile  or  so  in  defining  this  line  ..bove  the  layer  will  cause 
appreciable  errors  in  the  derived  layer  characteristics.  It  is,  there= 
fore,  much  better  to  use  refractometer  or  radiosonde  data  if  it  is  avail¬ 
able. 
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7.  RADIO  PERFORMANCE 


The  physical  concepts  of  radio  wave  propagation,  which  were  de¬ 
veloped  in  parts  I  and  II  of  this  manual,  apply  equally  well  whether  one 
is  concerned  with  radar  detection  or  with  radio  commur’ nation.  The 
application  of  these  concepts,  however,  is  somev'hat  differeui,  in  the  two 
cases.  Since  radar  detection  involves  reflects.  .1  of  energy  from  the 
Largct  and  a  return  trip  to  the  radar  antenna,  a  great  cieajL  moi’e  power 
is  required  for  the  detection  cf  a  target  at  a  given  distance  than  for 
radio  communication  over  the  same  distance.  Furthermore,  the  mete¬ 
orological  effects  needed  to  increase  radar  detection  ranges  are  much 
stronger  than  those  needed  to  extend  communications  ranges  by  the 
same  amoimt.  In  this  section,  the  principles  presented  in  parts  I  and 
II  are  applied  to  the  problem  of  radio  communications  over  the  ocean. 
Before  proceeding,  however,  it  is  necessary  to  become  familiar  with 
the  term,  “decibel”. 


Since  radio  power  or  field  strength  may  vary  over  a  very  large 
range,  covering  many  orders  of  magnitude,  it  is  convenient  to  express 
these  quantities  on  a  logarithmic  scale.  In  most  radio  propagation  work, 
the  power  or  field  strength  is  expressed  in  “decibels  below  free  space”. 
To  understand  this  and  similar  terms,  one  must  bear  in  mind  that  deci¬ 
bels  always  refer  to  a  ratio.  The  definition  of  decibel  (db)  is 

P  F 

db  =  10  log-^  =  20  log-^  ,  (7.1) 

where  P  is  power  (or  intensity)  and  ^  is  field  strength  (or  amplitude). 
The  zero  subscript  refers  to  a  reference  power  or  field.  If  the  refer¬ 
ence  is  the  free  space  power  or  field  at  a  given  distance  from  the  trans¬ 
mitter,  then  expressing  the  power  at  “x”  decibels  below  free  space, 
means  that  ten  times  the  g  '•?  the  rr.tio  of ''fr‘'e  space”  to  "measured 
power”  at  a  given  distance  is  equal  to  For  <»vample,  suppose  the 

measured  power  received  at  20  miles  from  a  transmitter  is  one-fifth 
of  the  computed  free  space  power  rt  20  miles.  Then-^  =  5,  but  log  5 
=  0.70;  thus,  on  a  decibel  scale  the  measured  power  iSpiO  (0.7)  =  7 
decibels  below  free  space.  In  the  same  example,  since  =  5,  then 
=■  xfF  =  2,24, but  log  2,24  =  0,35,  and  20(0.35)  =  7 decibels.  Thus, 
the  c  =cib..'  is  the  same  v.hether  the  power  ratio  or  the  amplitude 

ratio  is  used.  The  following  sections  discuss  the  effects  meteoro¬ 
logical  variables  on  radio  communications. 
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7,1  Surface-to-Surfacfe  Situations 


7,1.1  Linear  N-Profile 

Consider  now  the  case  of  two  ships  proceeding  away  from  each 
other  with  omnidirectional  antennas  radiating  at  a  wavelength,  X,  and 
located  H  feet  above  the  sea  surface.  How  does  the  rad*  "ignal  strength, 
as  measured  at  one  of  the  ships,  vary  wr  h  distance?  The  result  is 
shov/n  schematically  in  figure  7.1, 

At  very  close  range,  there  is  interference  between  the  direct  and 
surface-reflected  rays.  This  will  cause  the  field  to  pass  through  max¬ 
ima  and  minima  as  the  ships  move  apart.  At  distances  where  the  two 
waves  arrive  just  out  of  phase  the  received  signal  will  be  very  low. 
These  distances  are  d  =  W  where  n  is  an  integer  (1,  2,  3,  etc.).  The 
greatest  distance  where  a  minimum  will  occur  corresponds  to  n  =  1  or 
d  =  -2^  ■ ,  Suppose  for  example  that  H  =  100  feet  and  X  =  1  meter  (300 
me.);  then  d  =  1.0  nautical  mile.  At  twice  this  distance  the  signal  will 


z 


Figure  7.1.  Radio  Field  Strength  versus  Distance  lor  Vanous  Propagation  Conditions. 
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again  pass  through  a  maximum.  Since  the  two  waves  are  in  phase  the 
signal  amplitude  will  be  double  the  free  space  field  at  this  point  (pow<jr 
will  be  4  times  free  space).  Thus,  the  field  is  10  log  4  =  6  decibels 
above  free  space.  From  this  point  on,  the  signal  will  decrease  steadi.^,' 
with  increasing  distance.  When  the  distance  is  doubled  again  (d 
or  4.0  miles  in  the  example  above),  the  field  will  be  equ.  to  the  free 
space  field.  At  this  point  the  two  waves  arri  'e  90°  out  of  phase  with 
each  oxher ,  sc  the  measured  fielo  is  simply  that  of  either  wave  alone. 


As  the  ships  continue  to  separate,  they  will  reach  the  horizon  dis¬ 
tance.  Since  both  antennas  are  assumed  at  equal  height  this  distance  is 


dh=  2  -  2.2  / - 2——  ,  (7.2) 

V  1  -0.0209(-^ 

where  H  is  the  antenna  height  in  feet  and-^is  the  N-gradient  in  the 
layer  below  the  antenna.  For  H  =  100  feet  and-^^  =  12  units  per  1,000 
feet,  d  =  25.5  nautical  miles.  The  signal  strength  at  the  horizon  dis¬ 
tance  turns  out  to  be  very  nearly  constant  (when  expressed  as  db  below 
free  space)  over  a  wide  range  of  wavelengths  and  N-gradients,  Figure 
7.2  shows  the  horizon  signal  strength  versus  X  for  a  substandard  N-gra¬ 
dient  (8  units/1,000  ft.)  and  for  a  superrefractive  gradient  (30  units/ 
1,000  ft.).  The  solid  curves  represent  100-foot  antenna  heights,  and  the 
dashed  curves  represent  50-foot  antenna  heights.  It  can  be  seen  that  a 
value  of  24  ±  1  decibels  below  free  space  represents  the  horizon  signal 
(field)  strength  for  all  but  the  shortest  wavelengths  quite  closely.  The 
main  effect  of  changing  the  N-gradient  is  to  increase  the  horizon  dis¬ 
tance  according  to  the  above  expression. 


Beyond  the  horizon  distance,  the  radio  waves  entei  the  diffraction 
region  in  which  the  ettenuation  rate  (decrease  of  field  strengtli  with 
distance)  in  decibels  pei  nautical  mile  is  cojistai.t.  The  attenuation 
rate.  A,  in  this  region  is 


A  (db/naut.  mile)  =  0.72 


AM  \  2/3 

-0.0209-^j 


(7.3) 


whci  j  X  is  :hc  wa\elengvh  iu centimeters  and  is  the  N-gradient  be¬ 
low  antenna  heigh^.  Using  the  previous  exar'.ple  (  X=  iOO  cm.,-^^  = 
12units/l,000  ft.)  gives  arate  of  1.3  decibels  per  nautical  mile;  smaller 
X’s  and  smaller  N-gradients  will  increase  the  rate.  Table  7.1  shows 
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Figure  7.2.  Horizon  Field  Strength  versus  Wavelength. 


the  attenuation  rate  for  various  \’s  and  N-gradients.  This  constant 
attenuation  rate  will  control  the  fields  out  to  a  distance  of  about  50  miles 
(for  shipboard  antenna  heights). 

From  this  point  on,  the  radio  waves  enter  the  "scatter"  region, 
in  which  the  dominant  propagation  mechanism  is  scattering  from  atmos¬ 
pheric  inhomogcr.cit'''s.  The  attenuation  rate  in  this  region  has  been 
found  experimentally  to  be  about  l>.i5  of  aoecibel  per  nautical  mile  and 
is  reasonably  independent  of  wavelength.  Ine  scatter  rej^lon  is  not 
usually  of  particular  interest  in  shipboard  communications  because  tne 
fields  are  so  weak  that  very  high-power  transmitters  and  large,  highly- 
directional  antennas  are  required  to  obtain  a  usable  signal.  Shore  in¬ 
stallations,  between  islands  for  example,  often  use  such  measures,  d 
i  Uiu  ’^het  .scatter  comm  .'^'-cations  are  very  reliable.  The  fields  in  tem¬ 
perate  latitudes  are  usually  somewhat  stronger  in  summer  than  in  win¬ 
ter,  because  .he  presence  of  more  water  ''apor  results  in  more  intense 
scattering  parcels. 
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TABLE  7.1 


X 

(units/ 1.000  ft.) 

Ah 

cm. 

r\ 

O 

10 

12 

15 

20 

25 

800 

0.73 

0,71 

0.68 

0.65 

58 

0.50 

0.43 

300 

0.93 

0.90 

0.87 

0.84 

0.74 

0.64 

0.55 

lail 

1.17 

1.13 

1.09 

1.05 

0.93 

0,81 

0.69 

60 

1.58 

1.53 

1.47 

1.42 

1,25 

1.09 

0.93 

30 

1.99 

1.92 

1.86 

1.79 

1.57 

1.37 

1.17 

1-5 

2.51 

2.43 

2.34 

2.25 

1.98 

1.73 

1.48 

Attenuation  Rate  in  the  Diffraction  Region  (db/naut.  mile)  for  Linear  N-Gradients. 


7,1.2  The  Effect  of  Layers 

Thus  far,  the  effects  of  constant  N-gradients  on  communications 
have  been  considered.  Under  these  circumstances  diffraction  theory 
may  be  applied  quite  easily,  since  a  linear  N-profile  may  be  closely 
approximated  by  simply  changing  the  effective  earth’s  curvature.  In 
the  case  of  nonlinear  N-profiles,  with  surface  and/or  elevated  layers, 
theory  has  not  been  so  well  verified  by  experimental  evidence.  It  is 
necessary,  therefore,  to  resort  to  semiempirical  methods,  as  derived 
byGossard  and  Anderson  [14],  to  obtain  useful  results.  In  practice,  one 
is  not  interested  in  the  signal  strength  as  such,  but  in  how  much  addi¬ 
tional  communicating  distance  will  result  from  given  layer  character¬ 
istics.  This  cdditioral  distance  is  given  by 


AD  (naut.  miles)  =  (Dq-2.55\/h) 


105 


B  log  ANj^  +  C  log  hrp 


where  Dq  is  the  normal  communicating  distance  in  nautical  miles,  H  is 
the  antenna  height  in  feet,  _AN,  is  the  cliange  in  N  through  the  layer,  h,ji 
is  th--  l?y-3i  -tf-f/  height  i.i  fest,  and  A.  B,  and  C  are  given  in  figure  7.3 
as  functions  of  wa-''elength. 
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Numerical  Example 

Suppose  for  example  X=  100  centimeters,  H  =  75  feet, 
and  Dq  36  xiii  !he  layp'.  clraracteristics  are 

=  2,000  feet  and  ANj^  =  10.  Fto..*  fig\’re7.3;  A  = 
-116,  B  =  -62,  C  =  +63 

105 

AD  =  (35  -  22)  -  1)  =  32.5  nautical  miles  . 

Increasing  the  Jayer-top  height  to  3,000  feet  makes 


AU  =  20  nautical  miles 
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Equation  7.4  applies  only  over  a  limited  range  of  elevated  layer 
characteristics.  It  would  not  bo  used  under  strong  ducting  conditions; 
i.e„  when,  A  +  B  log  +  C  log  h^,  is  less  than  about  10  decibels. 
Under  these  conditions,  the  assumption  of  a  constant  signal  at  the  hor,. 
zon  distance  becomes  invalid.  For  conditions  between  standard  and 
moderate  ducting,  and  for  extended  ranges  within  100  miles  from  the 
transmitter,  equation  7.4  will  be  reasonably  reliable. 

When  layers  are  present  which  give  very  large  (or  negative)  AD's 
using  the  above  expression,  there  is  no  method  (empirical  or  otherwise) 
for  predicting  the  expected  extension  in  communicating  distance.  When 
trapping  layers  are  encountered,  the  fields  at  the  horizon  and  beyond 
arefoimd  to  be  nearly  at  the  free  space  level  for  great  distances.  When 
the  h,j,  and  ANj^  of  a  layer  exceed  the  trapping  criteria  shown  in  figure 
8.7,  it  is  safe  to  conclude  that  communication  signals  will  be  very  much 
stronger  than  usual,  and  that  communications  may  be  carried  out  be¬ 
tween  ships  at  much  greater  distances  than  normal. 

7.2  Surface-to-Air  and  Air-to-Surface  Situations 


The  meteorological  effects  on  surface-to-air  communications  will 
be  considerably  less  pronoimced  than  those  described  in  chapter  6,  The 
region  of  maximum  signal  enhancement  caused  by  atmospheric  layers 
extends  from  the  sea  surface  up  to  the  top  of  the  layer.  Since  strong 
effects  are  rarely  caused  when  the  layer  top  is  above  a  few  thousand 
feet,  the  aircraft  v/ill  usually  be  well  above  the  region  where  strong 
effects  occur. 

In  cases  where  appreciable  surface-to-surface  effects  occur,  the 
signals  above  the  layer  will  also  be  enhanced,  but  to  a  degree  which 
falls  off  rapidly  with  height.  Experimental  evidence  shows  that  incases 
of  strong  ducting,  the  field  for  the  first  few  thousand  feet  above  the  layer 
top  is  appreciably  stioagr-  than  when  no  duc^-  is  present.  Also,  when 
the  aircraft  is  at  high  altitudes  (30,000  -  40,000  ft.),  the  racio  horizon 
distance  is  increased  by  as  much  as  5  percent  of  its  standard  value  of 
1.23  (\^S^  +  where  and  are  the  two  terminal  heights. 

No  systematic  study  has  been  made  for  intermediate  altitudes. 
It  is  likely,  however,  that  the  degree  of  communication  range  extension 
will  n,.'.  a  fewpcrcen*  '.’''•enever  the  aircraft  is  more  than  a  few 

thousa,  .d  feel  above  the  layer  top. 
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In  the  absence  of  layers  capable  of  causing  appreciable  surface  ■ 
to- surface  signal  enhancsrnent,  the  surface-to-air  enhancement  will 
be  negligible.  One  can  conclude,  therefore,  that  the  meteorological 
effects  on  surface-to-air  communications  will  generally  be  very  much 
less  than  over  surface-to-surface  paths,  and  that  they  will  be  appre¬ 
ciable  only  when  the  aircraft  is  at  low  altitude  (such  that  surface-to- 
surface  criteria  apply). 

7,3  Air-to-Air  Situations 


When  both  terminals  of  a  communications  link  are  well  above  the 
surface,  the  effects  of  •weak  ele'vated  layers  become  more  important 
than  they  are  when  one  terminal  is  near  the  surface.  In  this  situation, 
radio  holes  can  plcy  a  significant  role  in  redistributing  the  transmitted 
energy. 

In  section  6.3  (air-to-air  radar),  it  was  noted  that  radio  holes, 
caused  by  relatively  weak  elevated  layers,  can  cause  radio  energy  to 
be  very  weak  in  the  radio  hole,  and  that  beyond  the  hole,  for  a  distance 
roughly  equal  to  the  hole  width,  the  energy  is  stronger  than  usual.  Thus, 
as  two  aircraft  move  apart,  the  signal  will  decrease  at  the  free  space 
rate  until  the  near  edge  of  the  hole  is  reached.  Then  it  will  decrease 
markedly  for  a  distance  of  many  miles,  until  the  far  edge  of  the  radio 
hole  is  passed.  At  this  distance,  the  signal  will  increase  markedly  in 
intensity,  pass  through  a  maximum,  and  then  will  gradually  fall  off  to 
the  free  space  values  at  a  distance  roughly  equal  to  the  width  of  the 
radio  hole.  Beyond  this,  the  layer  will  have  little  effect,  and  the  signal 
strength  will  again  decrease  at  the  free  space  rate  (6  db  per  octave  of 
distance).  The  methods  for  locating  the  boundaries  of  the  radar  hole 
which  were  developed  in  section  6.3,  are  of  course  equally  valid  for 
the  radio  hole,  since  the  phenomena  are  identical. 

It  should  bf-fe*  .^hasized  that  all  of  the  above  effects  occur  within 
the  radio  horizon  as  determined  by  d  =  1.33  .  This  dis¬ 

tance  for  high  altitude  aircraft  is  rather  large  (353  miles  at  20,000  f+ 
for  example),  so  free  space  attenuation  will  usually  limit  the  communi¬ 
cation  range  before  the  earth’s  shadow  becomes  effective. 


PART  IV 

METEOROLOGICAL  STUDIES 
OF  THE 

ATMOSPHERIC  REFRACTIVE  INDEX 
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8.  REFRACTIVE  INDEX  CUMATOLOGY 


In  part  II  of  this  manual  we  have  discussed  methods  for  describing 
the  refractive  properties  of  an  atmosphere  in  terms  of  various  refrac- 
tivity  profiles,  and  in  part  III  we  have  outlined  procedures  for  predicting 
radio- radar  performanct  in  terms  of  specific  refractivity  conditions. 

In  the  absence  of  specific  information,  and  particularly  for  plan¬ 
ning  purposes,  it  is  useful  to  develop  a  knowledge  of  the  climatology  of 
refractivity.  An  under  star  ding  of  the  refractive  properties  of  air  masses 
will  also  be  helpful  in  answering  questions  as  to  the  horizontal  extent 
and  time  persistence  of  refractive  layers.  For  these  r  easons  such  a 
climatology  is  presented  in  the  following  paragraphs, 

8,1  Seasonal  Distribution  of  Refractivity  Near  the  Surface 

If  we  examine  the  first  term,  the  “dry”  term,  on  the  right-hand 
side  of  equations  3.3  or  3,4,  we  see  that  it  is  proportional  to  air  density. 
From  typical  atmospheric  data,  it  can  also  be  shown  that  the  “dry” 
term  constitutes  at  least  60  percent  of  the  value  of  N,  For  these  rea¬ 
sons  the  refractivity  is  strongly  affected  by  pressure  changes. 

In  order  to  study  the  distribution  of  refractivity  at  the  surface,  it 
is  desirable  to  suppress  those  variations  produced  solely  by  surface 
topography.  For  tliis  purpose  we  use  a  technique  analogous  to  that  of 
reducing  station  pressures  to  sea  level.  Solving  the  expression  for  the 
standard  refractivity  profile,  equation  3.5,  for  Ng  we  have 

■KT  TVT  Ch 

Nq  =  Ng  exp. 

where  Ng  is  the  refractivity  at  the  e'lrth’s  surface,  Njj  is  the  refrac¬ 
tivity  at  height  h  above  the  surface  and  £  is  the  exponential  coefficient 
of  a  particular  Si.uicirird  i^rcfiie.  We  then  define  a  datum  refractivity, 
Nq,  as  follows 

ch 

Nq  =  Ng  exp.  s  ,  (8.1) 

where  Ng  is  the  surface  refractivity,  hg  is  the  height  of  the  surface,  and 
f'  is  again  the  exponenti'^l  coefficient  of  a  particular  standard  profile. 
■'hpc>.t'>Tr  refracti'  ityis  then  the  value  obtained  if  one  extrapolates  the 
surface  value  along  a  standard  refractivi^'y  profile  to  zero  height,. 
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In  order  to  facilitate  the  reduction  of  surface  refractivity  values 
to  the  zero-height  datum  figure  8,1  has  been  prejiared,  using  the  value 
of  c  =  0.1057  as  arbitrarily  adopted  by  Bean  and  Riggs  [3]  as  being  the 
most  representative  standard  profile.  Equation  8.1  then  becomes 


N  =  w 
^'O  *^s 


0.1057  hg 


where  Nq,  Ng.  and  h^  (in kilometers)  are  defii.  as  in  equation  8.1,  in 
use,  one  locates  the  point  on  figure  8.1  corresponding  to  the  surface  re¬ 
fractivity,  Ng,  and  the  surface  elevation,  hg.  By  laying  a  straight  edge 
through  the  point  and  parallel  to  the  set  of  sloping,  dashed  lines,  the 
datum  refractivity,  Nq,  may  be  read  from  the  ordiiiate  scale.  Other 
values  of  the  exponential  coefficient,  c,  could  be  represented  on  figure 
8.1  by  additional  sets  of  slopdng,  dashed  lines  of  different  slopes. 


Using  the  reduction  method  described  above  Horn,  Bean,  and  Riggs 
[15]  have  prepared  world- wide  maps  of  datum  refractivity  for  various 
months  and  seasons.  Figures  8.2  and  8.3  show  these  results  for  the 
months  of  February  and  August.  It  can  be  seen  that  Nq  varies  from 
390  in  humid  tropical  areas  to  290  in  desert  and  plateau  regions.  Con¬ 
tinental  interiors  and  mountain  chains  in  middle  latitudes  are  reflected 
by  low  values  of  Nq  as  compared  to  coastal  areas.  Such  major  climatic 
details  as  the  Indian  monsoon  and  orographic  rain  shadows  are  also 
indicated  by  the  Nq  contours. 

Figure  8,4  illustrates  the  climatic  response  of  surface  refrac¬ 
tivity  and  shows  a  world- wide  analysis  of  the  annual  range  of  monthly 
mean  values  of  Ng,  Climatic  controls,  such  as  the  prevailing  transport 
of  maritime  air  inland  over  the  west  coasts  of  North  America  and  Europe, 
are  indicated  by  relatively  small  annual  ranges  of  Ng  (20  to  30  units). 
On  the  other  hand,  a  range  of  40  to  50  units  or  more  along  the  east 
coast  of  the  United  States  reflects  the  frequent  invasion  into  this  area 
of  such  diverse  air  nr,  ^  u'  cont:n''  ''tal  ar^tv?  and  maritime  tropical. 
The  largest  annual  ranges  of  Ng  (90  units)  are  obc?rved  in  the  Sudan  of 
Africa  and  in  regions  affected  by  the  Indian  monsoon. 


8 . 2  Refractive  Index  Characteristics  of  Air  Mass  Types 

Air  masses,  having  approximately  homogeneous  horizontal  dis- 
trib.’iiofij  of  tPi'^peraiure,  La.aidity,  and  lapse  rate,  have  been  defined 
[6J  and  are  commonly  identified  by  the  practical  meteorologist  in  his 
study  of  daily  weather  development.  Nearly  uniform  characteristics 
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Figure  8.1.  Chart  for  the  Reductior  of  Ng-Values  to  No-Values. 
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Figure  8.2.  World-Wide  Distribution  of  Datum  Refractivity.  Nq,  for  the  Month  of 
February. 
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Figure  8.4.  World-Wide  Distribution  of  the  Annual  Range  o;  Monthly  Mean  Refrac- 
tivity  at  the  Surface,  Ng, 
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are  acquired  when  air  remains  over  a  large  source  region,  with  uniform 
surface  temperature  and  humidity  conditions,  until  equilibrium  is  reach¬ 
ed  through  the  action  of  convection  and  radiative  processes.  Because 
refractive  properties  are  also  determined  by  the  distribution  of  tem¬ 
perature  and  humidity,  ei  ch  air  mass  may  be  associated  with  charac¬ 
teristic  refractivity  values  and  profiles.  These  air  mass  characteris¬ 
tics  will  now  be  described. 

Because  it  is  desire  i  to  suppress  the  pressure-height  variations 
in  refractivity,  the  A-unit  presentation  of  refractive  index  (see  section 
4,4,3)  will  be  used  in  constructing  all  profiles.  The  profiles  thus  be¬ 
come  more  sensitive  indicators  of  changes  in  temperature  and  humid¬ 
ity. 


8.2.1  Continenta} -Arctic  and  -Polar  Air  Masses 

Continental-arctic  air  masses  are  produced  over  the  arctic  (and 
antarctic)  fields  of  snow  and  ice  and  over  the  snow- covered  portions 
of  the  continents,  where  the  wind  systems  are  predominantly  anticy- 
clonic  and  large-scale  subsidence  is  typical.  Physical  modification  of 
continental  arctic  air  as  it  moves  southward  over  the  continents  changes 
it,  by  warming,  into  a  continental-polar  air  mass  in  which  the  tempera¬ 
ture  regime  is  less  severe. 

The  refractivity  profiles  of  the  two  air  masses,  continental  arctic 
and  continental  polar,  are  similar  as  indicated  in  figure  8,5,  which 
shows  representative  profiles  for  both  summer  and  winter.  Continental 
polar  air  is  characterized  summer  and  winter  by  a  profile  showing  a 
very  slight  increase  in  refractivity  with  height.  Continental  arctic  air 
shows  the  same  slight  increase  with  height  above  one-half  kilometer, 
but  exhibits  a  sharp  increase  in  the  lowest  levels  above  the  ground. 
Large  surface  v  coji*’n?mal  arctic  air  indicate  extremely  low 
ground  temperatures  characteristic  of  n  pronounced  arctic  inversion. 

When  continental  polar  air  is  found  away  from  its  arctic  source 
region  (in  the  United  States  for  example),  it  has  usually  been  modified, 
and  the  strong  surface  inversion  has  been  erased,  A  nearly  linear  p’'c- 
file,  resembling  the  standard  refractivity  profile,  results  and  supf.r- 
refrcctlvj  effecis  arc-  .  -.duced.  On  the  other  hand,  ducting  occurs  re¬ 
latively  frequently  near  the  ground  in  the  arctic  source  region. 


8.2.2  Maritime- Polar  Air  Masses 


The  maritime  polar  source  region  is  in  reality  an  oceanic  zone 
of  transition  separating  the  tropical  maritime  source  from  the  polar 
and  arctic  continenUii  scarce.  Durmg  the  winter-,  maritime  polar  air 
over  land  exhibits  a  linear  profile  above  about  tne  1 -Kilometer  level 
with  a  decrease  in  refractivity  in  the  lower  layers  (see  fig.  8.6),  due 
to  the  drying-out  effects  of  an  over- land  trajectory.  Summertime  pro¬ 
files  over  land  and  oceanic  profiles  show  an  increase  in  refractivity 
near  the  ground  due  to  the  addition  of  water  vapor  (at  warmer  tempera- 
Pmes). 

As  a  conscouence,  v'inter  profiles  ove'-  land  generally  produce 
standard  or  even  substandard  refractive  conditions,  while  summer  and 
oceanic  profiles  tend  to  become  superrefr active. 
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Figure  8,6.  Representative  Profiles  for  Maritime-Pular  ^ir  Masses. 


8.2.3  Continental-Tropical  Air  Masses 

Source  recfions  for  these  air  masses  are  found  in  the  Northern 
Hemisphere  over  tne  iarge  land  .aasses  North  Africa,  Asia,  south¬ 
western  United  States,  and  southern  Europe  u,  rummer:  and  over  North 
Africa  in  winter, 

Refractivity  profiles  for  continental  tropical  air  show  a  decided 
increase  in  refractivity  with  height  and  typically  have  surface  value  , 

300  units,  as  shovm  in  figure  8,7.  A  positive  gradient  of  about 
1;  1  A'  Vii'iil'-  per  kilo!'ieiei-  exists  aloft.  The  extremely  low  surface  re¬ 
fractivity  values  reflect  the  prevailing  high  temperatures,  extensive 
vertical  mixing,  and  coctreme  dryness  typ^Lcal  of  the  desert  source  ra- 
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Figure  8.7,  Representative  Profiles  for  Tropical  Air  Masses, 


gions  of  this  air  liiass.  Continuous  subrefractive  effects  can  be  expected 
in  continental  tropical  air  in  view  of  the  large  positive  A-gradient, 

8.2,4  Maritime-Tropical  Air  Masses 

The  maritime  tropical  source  region  comprises  the  oceanic  areas 
on  the  equatorial  side  of  the  subtropical  high-pressure  belt  in  both  hem¬ 
ispheres, 

vbnwr.  in  figive  8,7.  'profiles  for  this  warm,  humid  air  mass 
showc  i;iS./  Aed  decrease  in  refractivity  with  height,  often  amounting  to 
40A-units  or  more  it  the  3-}dlometer  level.  Sc  dace  refractivity  values 
are  about  10  A-units  higher  in  summer  due  to  the  increased  insolation 
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over  Northern  Hemisphere  source  regions,  Superrefractive  effects  can 
be  expected  within  this  air  mass  both  summer  and  W’intsr,  but  thej'  can 
be  expected  to  be  more  pronounced  in  summertime  due  to  the  increased 
water  vapor  gradient. 

8.2.5  Superior  Air  Masse> 

Superior  air  is  essentially  a  high-levr'  air  mass.  .  '■hough  it  may 
appear  at  the  surface.  It  appears  most  frequei  cly  over  the  southwestern 
part  of  North  America  and  is  believed  to  be  created  largely  by  large- 
scale  subsidence  aloft. 

This  air  mass  is  characterized  by  a  near -linear  profile,  as  shown 
in  figure  8.8.  A  slight  decrease  of  refractivity  with  height  is  evidenced 


Figure  8.8.  Representative  Profiles  for  Superior  Air  Masses. 


-  150  - 


in  most  summer  profiles.  Winter  profiles  typically  show  about  the  s£ime 
rate  of  decrease  to  the  2 -kilometer  level,  above  which  a  slight  increase 
in  refractivity  is  observed.  Superior  air  mass  profiles  approximate 
the  standard  refractivity  profile,  with  only  a  slight  tendency  toward 
superrefraction.  Normal  propagation  at  microwave  frequencies  can  be 
expected  within  a  superior  air  mass;  ducting  is  most  unlikely  due  to  the 
dryness  of  superior  air. 

8.2.6  Air  Mass  Summary 

In  order  to  summarize  the  salient  refractive  characteristics  of 
each  of  the  several  air  masses  discussed  in  the  preceding  paragraphs 
table  8.1  has  been  prepared  for  convenient  reference. 


TABLJS  8.1 


Air  Mass 

A*  Profile 
Characteristics 

Summer 

Winter 

■ra 

A  at  3  km. 

m 

Continental 

Arctic 

IBB 

330 

340 

330 

Continental 

Polar 

Near -Linear 

330 

330 

320 

325 

Maritime 

Polar 

Near-Linear 
(above  1  km.) 

340 

330 

320 

330 

Continental 

Tropical 

Strong  Positive 
A-Gradient 

300 

340 

290 

330 

Maritime 

Tropical 

Negative 

A-Gradient 

380 

340 

360 

330 

Superior 

Near- Linear 

330 

320 

310 

310 

Typical  Air  Mass  Refractivity  Characteristics  Between  the  Surface  and  3  Kilometers. 


8.3  Mechanisms  of  Azr  f.'nss  Modification 


We  have  described,  in  very  general  terms,  the  refractive  prop¬ 
erties  of  the  six  most  commonly  ideT>tified  air  masses.  These  proper¬ 
ties  are  typical  for  each  air  mass  in  or  near  its  source  region.  As  an 
air  mass  moves  from  its  source  region  it  is  modified  by  various  phy¬ 
sical  mechanisms  which  come  into  play.  In  order  to  predict  the  prob- 
abl?.  ch-'."=‘ctpr)stics  of  such  c  '  -arisient  air  mass  we  shall  now  discuss 
each  ri  =riia:;isni  in  turn. 
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8.3.1  Radiation 

An  important  physical  process,  producing  nocturnal  atmospheric 
stratification,  is  the  continuous  loss  of  energy  by  long-wave  radiation 
from  the  surface.  For  an  average  earth  temperature  of  287®  K.  this 
radiaticin  has  its  maximum  intensity  at  10/u  in  the  infrared  portion  of 
the  electromagnetic  spectrum.  Most  of  this  outgoing  radiation  is  cap¬ 
tured  by  the  carbon  dioxide  and  v/ater  ■^npor  absoj. .  ’  ''n  bands  in  the 
atmosphere,  so  that  the  groimd  cools  after  c  mdown  at  a  relatively  slow 
rate.  Through  conduction  and  convection  processes  the  air  layer  next 
to  the  ground  also  loses  heat.  An  inversion  layer  several  hundred  meters 
thick  may  be  produced  in  this  fashion.  Conditions  favoring  the  formation 
of  strong  nocturnal  radiation  inversions  are  clear  skies,  low  humidity, 
and  light  surface  winds.  Large  humidity  gradients  are  also  common 
in  radiation  inversions,  so  that  superrefractive  conditions  and,  occa¬ 
sionally  ducting,  result  from  the  large  negative  refractivity  gradients 
produced. 

Over  ocean  surfaces  the  effects  just  described  do  not  occur  to 
any  appreciable  extent  since  as  surface  water  cools,  it  becomes  heavier, 
sinks,  and  is  replaced  by  warmer  water  from  below  ♦he  surface.  Thus, 
ocean  surface  temperatures  do  not  vary  appreciably  from  day  to  night. 

8.3.2  Turbulence  and  Convection 

During  the  day  surface  heating  is  produced  by  insolation.  Surface 
heating  is  at  a  maximum  on  hot,  clear,  summer  afternoons  when  the 
ground  may  reach  temperatures  of  70°C,  and  higher,  as  shown  by  Sin¬ 
clair  [21].  The  atmosphere  above  this  superheated  ground  surface  is 
warmed  by  conduction  and  small-scale  convection.  Eventually,  super- 
adiabatic  gradients  are  produced  in  the  air  layer  above  the  ground,  and 
overturning  occurs.  This  process,  together  with  the  turbulence  created 
by  winds  movii!i*u\r'' rough  land  surfaces,  produces  large- scale  atmos¬ 
pheric  mixing  u'hich  tends  to  destroy  vertical  stratifications.  Accord¬ 
ingly,  standard  or  subrefractive  profiles  are  the  rule  during  daytiwo 
hours  over  land  areas  with  clear  skies  and/or  where  the  terrain  is 
rough  and  moderate  to  strong  winds  exist. 

As  in  the  case  of  radiational  cooling,  the  effects  described  ab->ve 
ai  •'  iijt  ob^ervou  over  o  ean  surfaces.  Insolation  striking  a  deep  water 
surface  is  not  absorbed  in  a  shallow  surface  layer  but  rather,  pene¬ 
trates  to  moderate  depths.  Thus,  the  surface  temperature,  trioether 
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with  that  of  a  thick  layer  of  water,  is  raised  only  slightly.  Also,  the 
roughness  coefficient  of  even  heavy  seas  is  small  compared  to  that  of 
most  land  surfaces. 

8.3.3  Subsidence 

Subsidence  is  defined  as  a  large-scale  downward  motion  super¬ 
imposed  on  the  horizontal  air  motions  in  the  atmosphere.  Its  occur¬ 
rence  is  common  in  the  eastern  portions  of  the  oemiperm.*-.  c  anticy¬ 
clones  of  the  Northern  Hemisphere.  The  lower  .  .mit  of  subsidence  is 
well  defined,  and  the  base  of  the  subsiding  layer  is  usually  found  be¬ 
tween  800  and  900  millibars.  T1  e  subsiding  air  is  warmed  by  compres- 
sional  heating  and  is  characterized  by  a  temperature  inversion  and  a 
large  dew  point  depression,  as  shown  by  Petterson,  Sheppard,  and 
Priestly[17].  Sincethe  air  below  the  subsidence  layer  is  usually  moist, 
excellent  conditions  for  the  formation  of  superrefractive  layers  aloft 
are  present.  The  California  coast  of  the  XJnitBu  States,  under  the  eas¬ 
tern  end  of  the  Pacific  anticyclone,  is  noted  for  the  persistence  of  super- 
refractive,  elevated  layers  during  the  summer  months. 

8.3.4  Evaporation  and  Condensation 

The  addition  of  water  vapor  from  a  warm  water  surface  to  over- 
lying  cold  air,  takes  place  rapidly  because  of  the  higher  vapor  pressure 
of  the  water  surface.  At  the  outset  of  this  process,  shallow  superre¬ 
fractive  or  ducting  layers  may  be  formed  close  to  the  surface.  How¬ 
ever,  as  the  evaporation  process  continues  and  the  moisture  has  been 
mixed  through  a  fairly  thick  layer  (by  convective  processes)  the  uni¬ 
form  distribution  of  moisture  precludes  the  formation  (or  maintenance) 
of  superrefractive  layers. 

Such  evaporative  effects  may  also  occur  inland  when  a  cold,  dry 
air  mass  moves  over  moist  soil  or  heavily  vegetated  areas.  The  Florida 
Everglades,  for  example,  have  about  the  same  evaporation  rate  as  the 
ocean, 

8.3.5  Advection 

An  air  mass  moving  over  a  surface  of  dissimiliar  physical  char¬ 
acteristics  (warmer,  cooler,  wetter,  dryer)  will  be  modified  through 
the  action  of  surface  heating  O’'  cooling,  convection,  evaporation,  etc. 
The  rest  >s  iz  oe  particular  case  may  be  deduced  from 
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the  discussion  given  in  the  preceding  paragraphs  and  are  discussed  in 
considerable  detail  in  NAVAJER  50-1P-527  [8], 

8.4  Climatology  of  Ground-Based  Ducts 


For  the  purposes  of  this  section  atmospheric  ducting  (trapping) 
is  defined  to  occur  whenever  geometrical  optics  indicates  that  a  radio 
ray  leaving  the  earth’s  surface  is  sufficiently  refracts  ‘"at  it  eventu¬ 
ally  travels  back  toward  or  parallel  to  the  Cu  th’s  surface.  This  cri- 
terion  has  been  applied  by  Bean  [2]  to  several  years  of  radiosonde  ob¬ 
servations  from  stations  ty  pical  of  polar,  temperate,  and  tropical  cli¬ 
mates  in  order  to  derive  estimates  of  the  variation  of  the  occurrence 
of  radio  ducts  with  climatic  conditions. 

Approximately  3  years  of  radiosonde  data  typical  of  a  polar  cli¬ 
mate  (Fairbanks,  Alaska),  a  temperate  climate  (Washington,  D,  C.),  and 
a  tropical  climate  (3wan  Island,  W.  I.)  were  examined  by  means  of  a 
digital  computer  for  the  occurrence  of  ducts  during  the  months  of  Feb¬ 
ruary,  May,  August,  and  November.  The  percentage  occurrence  of  ducts 
is  shown  in  figure  8.9  for  these  three  locations.  The  maximum  occur¬ 
rences  of  13,8  percent  for  August  at  Swan  Island  and  9.2  percent  for 
Fairbanks  in  February  are  significantly  greater  than  the  values  observed 
at  other  locations  and  times  of  the  year.  The  Washington  data  display 
a  summertime  maximum  of  4,6  percent.  These  data  indicate  that  the 
tropical  zone,  maximum  incidence  occurs  in  late  summer  and  is  nearly 
three  times  the  temperate  zone,  summertime  maximum;  while  the  arctic 
zone  maximum  occurs  in  midwinter  and  is  about  twice  as  large  as  the 
temperate  maximum. 

The  results  of  an  analysis  of  refractivity  gradients  observed  during 
ducting  are  given  in  figure  8,10.  The  maximum  gradient  of  420  N-units 
per  kilometer  was  observed  during  February  at  Fairbanks,  Alaska; 
with  a  value  of  396  N-units  per  kilometer  observed  during  August  at 
Washington,  D.  C.  'i  ne  ^jrauient  oostrved  dux  ing  ducting  at  Swan  Island, 
W.  I.  reached  a  maximum  (285N-imits  per  km.)  ootn  uuring  August  and 
November.  The  annual  mean  values  of  N-gradient  appear  to  show  a 
small  latitudinal  trend  from  a  high  value  of  230  N-units  per  kilometer 
at  Fairbanks  to  a  value  of  190  N-units  per  kilometer  at  Swan  Island. 

Another  propsT-ty  of  radio  ducts  is  their  thickness;  typical  thick- 
ne'  ■;  ’Clu-'s  .>rc  shov.  i  in  figure  8.11.  Again  a  small  latitudinal  trend 
is  observed  as  the  median  thickness  increas''  s  from  66  meters  <:t  Fair- 
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Pigiirp  8.9.  Frequency  of  Occurrence  of  Ground  Based  Ducts. 
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Figure  8.10.  Typical  Refractivity  Gradients  of  Ground-Based  Ducts. 


Percentage  of  Observotions  That  Equal  or  Exceed  the  Ordinate  Value 
Figure  8,1 1.  Observed  Ground-B^.sed  Duct  Thickness. 


banks  to  106  meters  at  Swan  Island, 


We  may  summarize  the  data  of  figures  8,10  and  8,11  by  saying 
that  polar  climates  (Fairbanks)  are  characterized  by  shallow  layers, 
with  relatively  intense  gradients,  tropical  dims 'es  (Swan  Island)  are 
characterized  by  thick  layers  with  relatively  weak  gradients,  and  tem¬ 
perate  climates  (Washington,  D.  C,)  lie  in  between. 

It  is  interesting  to  consider  the  temperatui .  '.nd  humidity  distri¬ 
bution  that  accompsinies  a  radio  duct^  An  example  typical  of  each  sta¬ 
tion  is  given  in  figure  8,12,  The  polar  duct  illustrated  was  accompanied 
by  a  surface  temperature  of  -25®C,  with  a  strong  temperature  inversion 
and  a  slight  humidity  lapse,  indicating  an  inversion  associated  withwin- 
tertinie  cooling  of  the  air  next  to  the  ground.  The  temperate  example 
appears  to  be  typical  of  a  temperate  zone,  temperature  inversion.  The 
tropical  profile,  however,  shows  a  moderate  lapse  of  both  temperature 
and  humidity.  This  apparent  inconsistency  is  explained  by  the  strong 
lapse  in  vapor  pressure  which  is  associated  with  a  moderate  lapse  in 
temperature  and  humidity  when  the  surface  temperature  is  near  30®C, 
This  strong,  vapor  pressure  lapse  presumably  arises  from  evaporation 
from  the  sea  surface. 

The  structure  of  ducts  was  further  examined  by  studying  the  per¬ 
centage  of  the  total  N-gradient  in  each  duct,  contributed  by  the  gradient 
of  the  "dry”  and  "wet”  terms.  The  median  contribution  of  the  "dry” 
term  gradient,  summarized  in  table  8,2,  displays  strong  seasonal  and 


TABLE  8.2 


Month 

Fairbanks 

(Polar) 

W  ashington,  D,C. 
(Temperate) 

Swan  Island 
(Tropical) 

February 

102?: 

73.0% 

9.5% 

May 

40.5% 

33.5% 

2.0% 

August 

37.0% 

26.5% 

4.5% 

November 

62.0% 

-  -  -  —  --  -  1 

55.0% 

6.0% 

Median  Contribution  to-^^  for  Ducting  Conditions, 

-n  ah 
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Figure  3.12,  Temperature  and  Humidity  Profiles  for  Typical  Ground-Based  Ducts. 
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geographic  differences.  The  “dry"  term  contribution  decreases  from 
winter  to  summer  and  from  arctic  to  tropical  climates.  The  tropical 
ducting  gradients  are  then  at  least  90  percent  due  to  humidity  lapse, 
while  the  polar  wintertime  maximum  is  due  to  the  strong  temperatui*e 
inversion  which  is  associated  with  very  low  surface  temperatures.  In 
fact,  under  these  conditions  vapor  pressure  increases  with  height,  and 
the  dry  term  contributes  more  than  100  percent  of  the  ducting  gradient, 

8,5  Climatology  of  Superrsfractive  Elevated  ^vers 

Although  no  systematic  study  of  elevated  layers  has  been  com¬ 
pleted  on  a  scale  large  enough  to  support  statistically  valid  conclusions, 
a  few  general  statements  may  be  made.  The  great  majority  of  all  strong, 
superrefractive,  elevated  layers  are  associated  with  subsidence  in  the 
atmosphere.  Since  subsidence  involves  an  adiabatic  process,  a  given 
layer  may,  therefore,  be  associated  with  a  certain,  constant  potential 
temperature.  In  other  words,  the  layer  is  in  an  isentropic  surface. 
This  concept  may  be  used  operationally  to  estimate  the  horizontal  ex¬ 
tent  of  strong  superrefractive  layers  by  indicating  the  particular  heights 
most  likely  for  layer  occurrence  at  points  removed  horizontally  from 
the  known  sounding.  It  must  be  remembered,  however,  that  the  exis¬ 
tence  of  the  “same"  potential  temperature  alone  is  not  sufficient  evi¬ 
dence  for  extending  an  elevated  layer  horizontally  into  ‘  'xmknown' '  areas; 
other  considerations  must  also  indicate  a  “good"  probability  for  its 
existence. 

Inmaritime-polar  air  masses  over  the  western  Atlantic  and  North 
Pacific,  the  boundary  between  the  moist,  lower,  mixed  layer  and  the 
unmodified,  dry,  upper  layer  can  frequently  become  an  elevated  super¬ 
refractive  layer.  When  the  existence  of  such  a  layer  has  been  estab¬ 
lished  by  measu’'ement,  the  tops  of  the  convective  clouds  in  the  mixed 
lower  layer  can  be  assumed  to  be  at  the  height  of  the  layer  base. 
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9.  SYNOPTIC  STUDIES  OF  REFRACTr/E  PROPERTIES 


As  indicated  in  the  introduction^  no  practical  techniques  exist 
which  will  pei’mit  one  to  make  useful  predictions  of  the  refractive  state 
of  the  atmosphere  beyond  a  6'-  to  12-hour  range.  For  such  "short" 
periods  one  may  properly,  either  assume  no  change  in  measured  con¬ 
ditions,  or  predict  trends  in  the  modificat-  '**  of  obsei  ved  refractivity 
profiles  expected  as  a  result  of  forecast  meteorological  changes. 

One  approach  to  the  problem  of  Jonger-range  refractivity  fore¬ 
casting  lies  in  the  direction  of  predicting  the  detailed  pressure,  tem¬ 
perature,  and  humidity  structure  at  several  points  within  the  area  of 
interest,  and  from  these  predictions  calculate  the  required  refractivity 
profiles.  The  precision  realized  in  12-  to  36-hour  forecasts  of  P, 
and  RH  with  present-day  methods  is  not  sufficient  to  produce  useable 
refractivity  predictions  for  these  forecast  intervals.  This  approach, 
for  periods  in  excess  of  18  hours,  can  be  expected  to  lie  beyond  the 
"state  of  the  art"  for  many  years  to  come. 

Another  approach  would  be  to  attempt  to  develop  3 -dimensional 
models  of  refractivity  which  could  be  directly  associated  with  the  class¬ 
ical  synoptic  features  (fronts,  highs,  lows,  etc.)  and  v/hich  would  be 
sufficiently  representative  of  the  real  refractive  situation  to  permit 
useful  forecasts  to  be  made  from  a  knowledge  of  expected  synoptic  de¬ 
velopments.  In  order  to  explore  this  approach  to  the  refractivity  fore¬ 
casting  problem,  a  detailed  case  study  has  been  completed  by  Horn, 
Bean,  and  Riggs  [15]  and  is  described  in  the  following  paragraphs. 

9,1  Analysis  of  the  Storm  of  18-21  February  1952 

A  wintertime  storm,  displaying  strong  contrasts  of  continental 
polar  and  maritim,  U  aii ,  "vs  choser  fo^  study  of  the  refractivity 

model  problem.  This  particular  storm  comi.icr.cs'’  with  an  outbreak 
of  polar  air  across  the  Great  Plains  from  the  18th  to  the  21st  of  Feb¬ 
ruary  1952,  during  which  time  a  pronounced  cold  front  developed  and 
moved  rapidly  across  the  United  States. 

9-1  Description  of  the  Synoptic  Development 

At  the  start  of  the  period  a  frontal  system  had  come  in  from  the 
Pacific  Northwest  and.  was  composed  of  a  cold  front  extending  from 
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northern  Utah  southward  into  Arizona,  with  a  stationary  front  extending 
northeastward  into  Wyoming^  This  frontal  system  had  moved  very 
slowly  for  several  days  prior  to  the  start  of  the  sequence.  With  a  fresh 
outbreak  of  polar  continental  air  east  of  the  Rocky  Moimtains  this  old 
front  became  more  active  and  {moving  ahead  of  the  fast-moving,  polar 
continental  front  sweeping  across  the  Great  Plains)  was  reported  as  a 
squall  line  by  the  time  it  crot  sed  the  Mississippi  River  early  on  the 
morning  of  the  20th  of  February,  During  the  latter  part  of  the  period 
the  polar-maritime  cold  front/ squall  line  wasl  '•ated  in  tl  „oveloping 
warm  sector  of  the  polar-front  wave.  The  entire  frontal  ensemble  of 
cold  front,  polar-front  wave,  and  squall  line  had  moved  rapidly  to  the 
east  coast  by  the  morning  of  the  21st,  which  is  the  last  day  of  the  period 
studied. 

This  storm  sequence  is  a  rather  well  developed  one,  with  strong 
temperature  and  humidity  contrasts  between  the  air  masses  within  the 
warm  sector  and  behind  the  cold  front.  The  existence  of  this  strong 
contrast  was  one  of  the  main  reasons  for  studying  this  particular  period, 
because  if  refractivity  is  not  a  good  synoptic  indicator  in  this  case  there 
is  little  hope  of  its  being  one  in  more  subtle  synoptic  disturbances. 

9.1.2  Datum  Refractivity  (Nq)  Analysis 

In  order  to  supress  the  effects  of  topography,  the  datum  refrac¬ 
tivity,  Nq,  (see  section  8,1)  has  been  used  as  the  refractive  parameter 
plotted  in  all  horlzontsd  analyses.  Charts  of  Nq  were  prepared  from 
surface  observations  taken  every  12  hours  from  0130  EST,  February  18 
until  1330  EST,  February  20,  1952,  or,  in  other  words,  every  12  hours 
during  the  period  studied.  These  Nq  charts  are  presented  in  figures 
9,1  through  9.6,  Each  chart  ihows  the  Nq  field  with  the  U,  S,  Weather 
Bureau  frontal  ^nalysis  (surface)  superimposed;  conventional  frontal 
symbols  have  been  used  in  thesefigures.  The  cold  front  which  originally 
extended  from  Utah  southward  appears  to  be  rather  weak,  since  it  is 
not  associated  with  an  change  ‘across  the  front*  In  the 

early  stages  of  our  sequence  the  lack  of  air  mass  contrasts  is  evidenced 
by  the  slight  change  in  the  position  of  the  290  NQ-contour  (encircling 
the  west  Texas  and  New  Mexico  rrea)  as  the  frontal  system  moves 
through  that  area.  By  comparison,  the  cold  front  which  later  sweeps 
down  across  the  Great  Plains  has  a  rather  marked  Ng-gradient  across 
the  front;  due  in  large  measure  to  the  northward  flow  of  warm,  moist 
air  inio  u-'t  r/arm  seet  jr,  vi-,\h  had  formed  by  1330  EST  on  the  19th. 
It  IS  perhaps  significant  that  the  Nq- contours  show  the  various  frontal 
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systems  as  transition  zones  rather  than  as  sharply  defined  disconti¬ 
nuities. 

Figures  8,1  through  9.6  appear  to  indicate  that  Nq  does  reflect 
frontal  system  changes  and  that,  with  extensive  studies  of  many  cases, 
a  system  of  direct  objective  forecasting  from  this  parameter  might  be 
developed.  Although  this  conclusion  is  based  on  a  study  of  a  single 
storm,  it  appears  to  be  borne  out  by  prelimi  a^'y  analy-  ^  of  other 
frontal  systems. 

9.1.3  Vertical  Cross  Section  A-Unit  Analysis 

Another  aspect  of  this  case  study  is  illustrated  in  a  series  of 
charts  displaying  vertical  cross  sections  of  the  storm  of  18-21  P’ebru- 
ary  1952  in  terms  of  A-units,  which  are  largely  free  from  a  pressure- 
height  dependency  (see  section  4,4.3).  The  height  distributions  of  re- 
fracti\dty  along  the  base  line  of  the  particular  cross  section  shown  in 
figure  9.7  were  computed  for  every  12  hours  during  the  4-day  period. 
Examples  representing  the  early,  mature,  and  late  stages  of  the  out¬ 
break  are  reproduced  as  figures  9,8,  9.9,  and  9,10,  At  the  outset  of  the 
period  of  observation,  the  A-unit  gradient  was  quite  flat  around  the  polar 
front.  By  1000  EST  on  the  18th  of  February,  figure  9,8,  the  polar  front 
was  located  midway  between  Glasgow,  Mom,  and  Dodge  City,  Kans. 
The  contrast  between  the  southerly  extension  of  polar  air  and  the  nor¬ 
therly  advection  of  tropical  maritime  air,  from  the  Gulf  of  Mexico  into 
the  developing  warm  sector  of  the  polar-front  wave,  is  evidenced  by 
the  relatively  large  A-gradients  in  the  vicinity  of  Dodge  City.  The  core 
of  tropical  maritime  air  has  evidently  not  progressed  far  enough  north¬ 
ward  to  displace  the  warm  but  dry  air  that  had  been  over  the  Great 
Plains  prior  to  ‘he  outbreak.  As  a  result,  a  region  of  low  A-values  is 
found  between  the  front  £ind  the  tropical  maritime  air.  Twenty-four 
hours  later,  figure  9,9,  the  core  of  tropical  maritime  air  has  become 
more  extensive  and  now  I'e  --  her.  to  an  .’».Hitude  cf  ^  kilometers.  A  second 
cold  front  was  now  reported  on  the  daily  weather  ’>’ap,  and  the  area  of 
low  A-values  is  confined  between  it  and  the  main  cold  front.  Finally, 
by  the  morning  of  February  20,  figure  9.10,  the  front  has  passed  to  the 
south  of  Lake  Charles,  La,  and  the  polar  air  just  behind  the  front  is 
characterized  by  relatively  low  A-values. 

9.1.4  Lpj';  Cha»’t  Aiiaiysis 

The  use  of  space  cross  sections  does  not  always  yield  measure- 
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Figure  9.7,  Location  of  Cross-Section, 


Figure  9.8.  Space  Cross-Section  in  A-Units,  1000  EST,  18  February  1952. 
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ments  at  the  most  desirable  points  along  a  frontal  interface.  Another 
method  of  determining  the  refractive  index  structure  in  a  frontal  zone 
is  to  construct  a  time  cross  section  from  the  series  of  radiosonde  ob¬ 
servations  taken  at  a  single  station.  Thus,  as  the  frontal  system  ad¬ 
vances  and  passes  over  the  station,  one  obtains  a  time-history  of  the 
structure  of  the  frontal  zone.  Such  cross  sections  may  be  drawn  using 
time,  or  the  distance  from  the  station  to  the  front  (at  the  surface),  as 
abscissa.  In  the  latter  case,  the  presentat-  m  is  referr''^’  to  as  an  epic 
chart  since  the  observations  are  normalizeo  ^ith  respect  to  the  frontal 
passage.  Epic  charts  are  illustrated  in  figures  9,11  and  9,12.  Figure 
9.11  represents  atypical  continental  station  located  in  the  polar  conti¬ 
nental  air  throughout  the  storm  sequence  under  study.  The  essential 


Figure  9.11.  Fpic  Chart  in  A-Units,  R&..id  City,  South  Dakota. 
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feature  here  is  the  absence  of  detail  in  the  A- structure  due  to  the  presi- 
ence  of  a  uniform  air  mass  over  the  station.  Compare  this  with  the 
epic  chart  for  Oklahoma  City,  Okla.,  figure  9,12,  which  clearly  shows 
a  strong  prefrontal  A-imit  high,  a  strong  gradient  across  the  frontal 
interface,  and  an  A-unit  low  behind  the  front.  Figure  9.12  was  cor  • 
structed  for  a  station  located  it  the  point  of  confluence  of  virtually  un¬ 
modified  polar  continental  and  tropical  maritime  air, 

9.2  Conclusion 


Although  no  general  conclusions  should  be  drawn  from  this  single 
case  study,  it  appears  probable  that  3-dimensional  refractivity  models 


Figure  9.12.  Epic  Chari  in  A-Units,  Oklahoma  City,  Oklahoma. 
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could  be  developed  which  could  be  directly  associated  with  synoptic  fea¬ 
tures  and  which  were  generally  representative  of  the  refractive  prop¬ 
erties  of  the  real  atmosphere.  With  this  object  in  mind  preliminary 
models  for  an  idealized,  continental,  fast-moving  cold  front  and  for  an 
idealized,  continental  warm  front  have  been  constructed  by  Bean,  Riggs, 
and  Horn  [4]  and  are  shown  in  figures  9.13  and  9.14.  The  figures  are 
drawn  in  the  form  of  cross  sections  normal  to  the  front  and  show  the 
frontal  surfaces,  the  directions  of  motion,  'nd  the  refra^’^ivity  fields  in 
the  form  of  A- unit  contours.  Great  care  musi  beexerc  oed  in  the  appli¬ 
cation  of  these  models,  until  such  time  as  the  details  of  their  construc¬ 
tion  can  be  verified  by  tht  analysis  of  a  statistically  valid  number  of 
cases. 


R?f;’acti  e  imit-x  Structure  Through  an  Idealized  Continental,  Fast- 
Moving  Cold  Front  in  A-Units. 
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Figure  9.14.  Refractive  Index  Structure  Through  an  Ideali:^ed  Continental  Warm 
Front  in  A- Units. 
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